Microscopic Models for Ultrarelativistic Heavy Ion Collisions 



S. A. Bass, M. Belkacem, M. Bleicher, M. Brandstetter, L. Bravina, C. Ernst, 
L. Gerland, M. Hofmann, S. Hofmann, J. Konopka, G. Mao, L. Neise, S. Soff, 
C. Spieles, H. Weber, L. A. Winckelmann, H. Stocker and W. Greiner 

Institut fur Theoretische Physik 
Robert-Mayer-Strasse 10 
Johann Wolfgang Goethe - Universitat 
D-60054 Frankfurt am Main 

Ch. Hartnack and J. Aichelin 

SUBATECH, Ecole des Mines 
4, rue Alfred Kastler 
La Chantrerie 
F-44070 Nantes, Cedex 03, France 

N. Amelin 

Joint Institute for Nuclear Research (JINR) Dubna 
141980 Moscow region 
Russia 



February 9, 2008 



Contents 



I Introduction 



2 Transport theory of nuclear collisions 


2.1 


Non-relativistic transport theory 


2.2 


Rclativistic transport theory 


2.3 


The Quantum Molecular Dynamics approach 




2.3.1 Formal derivation of the transport equation 




2.3.2 Inclusion of collisions 




2.3.3 Pauli blocking due to Fermi statistics 




2.3.4 The Relativistic Quantum Molecular Dynamics approach 



3 The UrQMD-Model 




3.1 Initialization 


3.2 Equations of motion .... 




3.3 The collision term 


3.3.1 The collision criterion 




3.3.2 Cross sections . . . . 




3.3.3 Angular distributions 






3.3.4 Resonance lifetimes 

3.3.5 String-excitation and -fragmentation 

3.3.6 Color fluctuations, color opacity and transparency 



4 Probing hot and dense nuclear matter with relativistic heavy ion collisions 

4.1 Baryon stopping 

4.2 Particle production 

4.2.1 Subthreshold particle production from resonance matter 

4.2.2 'temperatures and single particle spectra 

4.2.3 Hadrochemistry: multiplicities and ratios 

4.2.4 Equilibrium properties of hadronic matter 

4.2.5 Resonance matter 

4.2.6 Antimatter and strange matter 

4.2.7 Particle freeze-out" 

4.2.8 Dilcpton production 

4.2.9 Charmonium production and suppression 

4.3 Collective flow 

4.3.1 Bounce-ofl: collective flow in the reaction plane 



4.3.2 Squeeze-out: flow perpendicular to the reaction plane] 
5 Summary and conclusions 



105 
113 



2 



Chapter 1 
Introduction 



The investigation of nuclear matter under extreme conditions is one of the major research topics of nuclear 
and high energy physics (a short and necessarily incomplete overview may be obtained in refs. |L[ [2j, || |], [|, 
H; H; i; H HH HH ) • The motivation for this is the unique opportunity to investigate the equation of state 



of nuclear matter and the search for phase transitions (such as the liquid gas or the quark gluon plasma 
(QGP) phase transition) and the possible restoration of chiral symmetry. Also the general understanding 
of the dynamics of heavy ion collisions over a vast energy range from the Coulomb barrier (several MeV per 
nucleon) to the highest energies currently available or planned for the future is interesting in itself. Here 
one can check the current understanding of the theory of strong interactions (QCD) and different effective 
theories based on hadronic/quark degrees of freedom. The intriguing role of color coherence phenomena 
(transparency and opacity) for fluctuations, stopping and charmonium production can be studied best at 
these energies. 

The measurements above several hundreds of AMeV are done at experimental heavy-ion facilities in 
three energy regimes: i) at about 1 AGeV at BEVALAC in Berkeley or SIS (SchwerlonenSynchrotron) at 
GSI-Darmstadt; ii) the AGS (Alternating Gradient Synchrotron, in Brookhaven) energy regime at about 
2-15 AGeV; Hi) and the SPS (Super Proton Synchrotron, at CERN) energy regime of 40-200 AGeV. 
Much higher energies will be available in the future with the Relativistic Heavy Ion Collider (RHIC) in 
Brookhaven (^i « 200 AGeV) and the Large Hadron Collider (LHC) at CERN (yG « 6 ATeV). 

Unfortunately, there exists presently no theoretical model that provides a consistent understanding 
of the reaction dynamics of heavy ion collisions over the whole energy range. One has to deal with 
quite different reaction mechanisms from compound nucleus formation and deep inelastic scattering at the 
Coulomb barrier over particle and resonance production at intermediate energies up to string-excitation 
and -fragmentation or parton scattering at (ultra-) relativistic energies. While at low and intermediate 
energies, descriptions in terms of hadrons (resonances) are appropriate, at high energies the quark and 
gluon degrees of freedom enter the game. 

One aim of this paper is the description of a theoretical model that incorporates these different reaction 
mechanisms and that is able to yield observables. The model is dubbed UrQMD. Such a microscopic model 
is based on a phase space description of the reaction. It contains a lot of unknown parameters, which 
will have to be checked and fixed by experimental data or by further model assumptions. This theoretical 
approach allows to pin down physical ingredients that determine the values of certain observables like 
particle abundances, collective flow of hadrons, rapidity distributions, cluster formation, etc. Thus certain 
experimental observables can either be uniquely traced to a physical parameter or if they are described 
equally well by different physical assumptions. This is of utmost importance if one wants to find evidence for 
new physical phenomena like the phase transition to the quark gluon plasma. It must be clarified whether 
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the measured data could also be understood in terms of a purely hadronic scenario or by assuming only 
a partial dynamical deconfinement, which is not accompanied by complete equilibration. One example is 
the vividly discussed suppression of J/^ or ty' production in Pb+Pb collisions as compared to p+p or 
p+A collisions, which can be explained in different physical scenarios. 

It is one of the main tasks of theoretical heavy ion physics to link experimental observables to the 
different phases and manifestations of nuclear or - in more general terms - hadronic matter. For this, a 
detailed understanding of the dynamics of heavy ion reactions is essential. Transport theory has played 
an important role in the interpretation of experimental results and in predicting new interesting effects 
in relativistic heavy ion reactions. It is particularly well suited for the non-equilibrium situation, rapid 
time-dependence of the system (even the use of the term "state" seems questionable), finite size effects, 
non-homogeneity, iV-body phase space, particle/resonance production and freeze-out as well as collective 
dynamics. Microscopic and macroscopic (hydrodynamical) transport models attempt to describe the full 
time-evolution from the initial state of the heavy ion reaction (i.e. the two colliding nuclei) up to the 
freeze-out of all initial and produced particles after the reaction. Simplified thermal equilibrium models 
neglect most of these dynamical effects, but make physical assumptions on the initial part of the reaction, 
e.g. thermalization or plasma creation. 

It should be kept in mind that notions like 'phase transition' are in principle restricted to equilibrated 
systems in the thermodynamic limit (N, V — ► oo, N/V =const.). However, in heavy ion collisions we have 
to deal with finite, dynamical systems, which makes the applicability of thermodynamic models doubtful, 
even though they might successfully describe some global observables like particle abundances. Never- 
theless, thermodynamic concepts are quite useful for a global characterization of the hot and compressed 
hadronic matter. In figure |1.1| we therefore present a conjectural phase diagram of nuclear matter where 
only the point at T=0, po = 0.16 fm~ 3 , realized in nuclei, is well known. For low excitation energies (low 
densities p < po and low temperatures T < Tc ~ 15 MeV), one expects a phase transition from a hadronic 
gas to liquefied nuclei. This possibility results from the resemblance of the nuclear interaction to that of 
atoms in gases. Both exhibit a short range repulsion and a long range attraction. There are indications 
that this phase transition may have been recently observed ]13| . 

At higher densities (2-3 p ) the possible existence of density isomers was proposed |T^, Q. At even 
higher densities and/or temperatures, one expects a phase transition or even a smooth cross over to 
the QGP. However, neither the exact critical density and temperature are known, nor the order of the 
transition. The latter is important, since a coexistence region of plasma (strongly interacting quarks and 
gluons not confined to well separated hadrons) and hadron gas might lead to observable effects in the flow 
excitation function [|l^, [7], [16|, [17|, [18| and other observables. 

Figure LI shows the generalization of the phase diagram to net-strangeness and net-antibaryon. Hy- 
permatter and antimatter can only be produced in energetic heavy ion collisions. The production of 
(multi-) strange nuclei containing one or more A-particles fl9], |2(| |2l|, [22H , as well as the production of 
anti-nuclei (anti- matter cluster) have been measured p3| , |24]j . Furthermore, the possibility of metastable 
exotic multi strange objects (MEMOS) |2l| |6|] and strangelets [^7], ^8], ^0], [31], |32], |33| might serve as 



observables for deconfinement during heavy ion collisions. Ultra-relativistic heavy ion collisions offer the 
unique opportunity to probe these different states of highly excited hot and dense nuclear matter under 
controlled laboratory conditions. The dashed lines in figure |1.1| visualize how different regions in the phase 
diagram may be probed simultaneously in the course of one heavy ion collision event, due to distillation 
and local fluctuations . 

Note that in a hot and compressed hadron gas at high excitation energies we deal not only with nucleons 
but with many different baryon (> 50) and meson (> 30) species. With increasing mass, these particles are 
often unstable - the width of the resonances increases. Many of the presently available microscopic models 
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Temperature 




Figure 1.1: Phase diagram of hadronic matter: Only the point at ground state density is well known; the 
dashed lines show areas probed in the course of heavy ion reactions. 



treat resonances as quasi-particles which propagate and scatter during their lifetime and then decay into 
other resonances or stable hadrons (with respect to strong interactions). This concept becomes doubtful 
when the width of the resonances is so large that no well defined quasi particle state (i.e. width/mass 
~ 1) exists. However, the resonances turn out to be important for the microscopic description of heavy ion 
collisions at bombarding energies above 1 GeV per nucleon. On the other hand, the concept of resonances 
cannot be used for internal excitation energies above 2 GeV. Here the quark degrees of freedom must be 
explicitly taken into account. Therefore some of the models presently available, including the UrQMD 
model presented in detail in this review, replace the resonances by continuous string excitations for higher 
energies. 

A microscopic dynamics description of heavy ion collisions is usually based on transport theory. Here 
a sequence of propagations of particles is simulated numerically (representing baryons (resonances) and 
mesons with or without interaction with a potential and subsequent scattering processes or decays). The 
main ingredients in this description of heavy-ion collisions are the cross sections, the two-body potentials 
and decay widths. It is obvious that - at least in principal - these models should be obtained consistently 
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from an underlying theory. Since the particles propagate in a hot and dense and even not equilibrated 
medium of highly excited hadrons, the properties of the particles might change significantly. Consequently, 
properties like effective masses, effective momenta, in-medium cross sections and decay widths should be 
calculated for the actual local situation in which the particle propagates. Unfortunately this is very 
complicated and from the numerical point of view also very time consuming, so that in most of the models 
drastic approximations are made. 

Furthermore, the basic requirements for the applicability of transport theoretical models might not be 
fulfilled in certain realistic situations. First of all the quasi particle limit, i.e. on shell propagation between 
subsequent scatterings, can be a bad approximation in strongly interacting, dense matter. Between two 
subsequent scatterings of a particle with time difference At, the energy of the state of the particle can only 
be fixed up to AE = h/At. For short collision times this uncertainty can become rather large (At = 1 
fm/c — > AE = 200 MeV). Even without collisional broadening the quasi particles acquire a finite life time, 
which requires the introduction of nontrivial spectral functions. Consequently one needs to introduce 
the energy of the quasi particles as an additional independent variable, i.e. the full 8-dimensional phase 
space description of transport theory is required. Secondly, transport theory relies on the description of 
the dynamics of the system as a Markovian process. This means the particles completely 'forget' about 
their critical state before each scattering. At least in the energy regime of several GeV per nucleon the so 



called memory effect does not seem to be negligible ||35| . Moreover, pure transport models describe the 
system by the one body phase space distribution function, which does not contain any information about 
correlations of the particles. These are, however, the most important for the description of fragment or 
cluster formation (hypernuclei, antinuclei, etc.). Such clusters can be obtained from the one body phase 
space distributions only with the help of a statistical coalescence model. In the UrQMD model described 
below, this shortcoming is avoided, since the particles can interact by individual two body forces at least 
at low energies. 

Unfortunately, the generalization of such two body forces to the relativistic region is not a simple task 
and is not incorporated in the present UrQMD model. In principle, the interaction must be mediated by 
fields, which are propagated according to wave equations or one must make use of the so called constraint 
Hamiltonian dynamics |3?|. Many other models avoid the propagation of fields using local density 
approximation. 

Finally, we would like to mention an additional problem arising from the practical realization of trans- 
port models: Even though the underlying transport equations are relativistically invariant, the actual 
realization by a sequence of propagations, scatterings and decays is not (see section |3.3.1| ) . This leads to 
a frame dependence of the observables which contributes to the systematic errors of the respective model 
predictions. 

This review is organized as follows: Chapter |2] contains a brief introduction into the many-body theory 
of nuclear collisions. In chapter ||] we discuss in great detail one specific microscopic transport model, the 
Ultra-relativistic Quantum Molecular Dynamics (UrQMD) model. Chapter |4] focuses on the exploration of 
hot and dense nuclear matter with microscopic transport approaches. Different effects and observables are 
discussed and dependencies to model parameters are pointed out. The main focus lies on baryon stopping, 
collective flow and particle production. A brief summary and conclusions are found in chapter |5|. 
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Chapter 2 

Transport theory of nuclear collisions 



2.1 Non-relativist ic transport theory 

The ultimate goal of a theoretical description of heavy ion collisions is a derivation and solution of equations 
of motion for the elementary degrees of freedom, the quarks and gluons. Since non perturbative dynamical 
solutions of QCD are not possible yet, one is forced to separate the description into two parts: The low- 
energy part (y/s < several GeV), where a description in terms of QCD-quasi particles (hadrons) is possible 
and a high energy part where the quark and gluon degrees of freedom are treated explicitly. For the 
low energy part one can derive solvable equations of motion from phenomenological effective Lagrangians, 
if strong but controllable approximations are made. This means that at least in principle one knows 
what has been neglected in the approximation scheme as long as one accepts an effective Lagrangian with 
baryons and mesons as degrees of freedom as a starting point. One example is the a — u— model p8| . |3£| 



where the interaction of nucleons is described by the exchange of a scalar and a vector meson. While the 
cj-meson has a clear physical significance, it is not quite clear if the cr-meson can be interpreted as a real 
physical particle or rather must be taken as an effective description of the interaction. This model must 
be extended to include other baryonic and mesonic degrees of freedom (e.g. A, iV*, ir, K, p etc.) since 
already at moderate energies of several hundred MeV per nucleon a significant production of new particles 
sets in (mainly pions but also other (pseudo-) scalar and (pseudo-) vector mesons as well as a lot of baryon 
resonances and even anti-baryons). The general theoretical framework for the derivation of practically 
solvable equations of motion based on effective Lagrangians is transport theory. The basic approximation 
which is made here is the description of the dynamics in terms of a semi-classical single particle phase 
space distribution instead of iV-body non-equilibrium Green's functions. 

In order to sketch the corresponding approximation scheme, we start with a simple non-relativistic 
potential model for nucleons and then later on simply write down the proper relativistic generalizations. 
The Hamilton operator shall be given by 

H = Jd 3 x ^(aji , t) (-^Ai^ i> H (xi ,t) + \j d^i d 3 x 2 d 3 x[ dV 2 

(xi,x 2 \V\ x[ , x' 2 ) ^h(xi , t)ip j H (x 2 , t)ip H (x' 2 , t)tp H (x[ , t) , (2.1) 

where the nucleon field operator ip\j{x\,t) and ipn(xi ,t) are given in the Heisenberg representation. For 
simplicity we neglect here the complications introduced by a spin or isospin dependent interaction as 
well as three body forces, since we only want to review the principle steps of the derivation. For details 
we will refer to the corresponding literature as well as to the derivations in the relativistic case. The 



7 



non-relativistic non-equilibrium Green's functions (g c ,g a ) and correlation functions (g > ,g < ) are defined 
by: 

/(MO = ^<TV*(1)A(1')>. 

^(MO = 4(^(1)^(10), 

^(M0 = -^(^(10^(1)), (2.2) 

with 1 a short hand notation for (xx,tx) and/or spin or isospin indices and T c ,T a are the chronological 
and anti-chronological time ordering operators. In the relativistic case ip* must be replaced by ■0. The 
four two point functions can be written in a more compact form as a path ordered Green's function (2x2 
matrix denoted by the underline) 

G(1,1') = 4(P^(1)^(10)- (2-3) 

Here, additionally to the times ti,t±', one has to specify if these times are on the upper or lower branch 
of the integration contour, which yields the four combinations given in Eq. ( |2.2j ). Originally, the theory 
was formulated by Schwinger [40| and the first application to transport problems is due to Kadanoff and 



Baym For an introduction to the formalism of path ordered non equilibrium Green's functions we 

refer the reader to the works of Danielewicz [42] and Botermans and Malfliet [43||. To some extent we will 
follow their work and notations. From the equations of motion for the Heisenberg field operators one can 
derive equations of motion for the Green's functions. However these equations for the one body Green's 
functions couple to the two body Green's functions and so on. This set of coupled equations of motion is 
usually called the Martin Schwinger hierarchy j44|. In order to solve the problem completely one has to 
solve the whole set of coupled equations for all n-body Green's functions. However, one can replace in the 
equation for the one body Green's function the unknown two body Green's function by the also unknown 
self-energy. This leads to the Dyson equation in differential form: 

s(l) G(l, 1') = 5(1 - 10 + / dl"E(l, 1") G(l", 1') (2.4) 
here s(l) is the differential operator 



c 



2m 

The 5-function is generalized to the integration contour 



h 2 

s(l) = ihd tl + —A X1 . (2.5) 



c 

such that /dl'£(l - l')Z(l') = £(1), namely: 

c 



CXJ CXJ 

Jdl = J dl (upper) - J dl (lower) (2.6) 



5{t\ — t 2 ) both ti and t 2 on the upper branch 
5(ti — t 2 ) = { —5(ti — t 2 ) both ti and t 2 on the lower branch . (2.7) 
otherwise 
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Of course the proper self-energy 1') is also a two by two matrix defined on the contour. In more 



detail, equation (2.4) has four components: 



^(1)^(1,1') 



5(1 -!') + / dl" P c (l , 1") g c (l", 1') - E<(1, 1") g>(l", 1')] , 
-5(1 -!') + / dl" P > (1 , 1") g^l", 1') - S a (l, 1") «f (1", 10] 
| dl" [£>(1 , 1") s c (l", 10 - £ a (l, 1") <7>(1", 10] , 

J dl" p c (i , l'O g^i", 10 - s<(i, l") g a (i", l')] , 



(2.8) 



which are known as Kadanoff-Baym equations j}fl ]. In integral form equation ( |2.4| ) is the Dyson equation 
on the contour 

10 = Go(l, 10 + / dl" J dl'" Go(l, 1"') £(1"', 1")G(1", 1) (2.9) 



c c 

where the zeroth order Green's function inverts s(l) 



s(i)G o (l,i0 = 5(1-10 



(2.10) 



The equations of motion ( |2.4|) or ( |2.8| ), respectively, contain in principle the unknown self energies which 
have to be calculated by an appropriate approximation scheme, e.g. in perturbation theory or in T- or G- 
matrix approximation. However, the corresponding expressions will in turn contain the unknown Green's 
functions, such that one obtains a coupled set of equations for the G's and S's. E.g. in the T- or G-matrix 
approximation one calculates the self energies from 



S(l, 1') = -ih J d2 J d2' ( (1 2|T|1' 2') - (1 2|T|2" 1') ) G(2", 2 + ) 



(2.11) 



c o 

where the T-matrix is defined by the integral equation: 



(12|T|1'2') = (12|Z|l / 2 / ) 

+ ih J dl"d2"dl'"d2'"(12|V|l ,,, 2"')G(l",l'")G(2",2'") (1"'2"'|T|1'2') (2.12) 

c 

For the effective interaction in nuclear matter (G-matrix), one has additionally to incorporate the Pauli 
blocking of the intermediate states where the two particles or holes propagate from 2"' to 2" and from 



1'" to 1", respectively. For details, we refer the reader to [43, 5S|. In the dynamical situation of heavy-ion 
collisions, the problem is more difficult. At the beginning of the reaction, one does not have a single 
Fermi sphere at zero temperature as for nuclear matter in the ground state, but rather one has to deal 
with two Fermi spheres separated by the beam momentum. Calculations of the G-matrix for colliding 
nuclear matter have so far not been carried out. However, the effects of the dynamical situation on the 



mean-fields (real part of the self-energy) have been estimated by the Munich group [46]. Only a few 
attempts have been made to include realistic G-matrices into a transport description, see e.g. the work of 



the Tubingen group |47|]. In many other cases, free or in medium corrected cross sections are used together 
with parameterizations of the mean-field (e.g. different versions of the QMD approach) or the latter is 
calculated in local density approximation (see below). In order to avoid confusion with the propagator 
G(l, 1'), we use the term T-matrix or T instead of G-matrix or G. 
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Equation ( |2.12| ) corresponds to a resummation of all the ladder diagrams for the effective interaction 
(T-matrix) in the medium, which is an appropriate approximation for interactions with hard core, like the 
nucleon-nucleon interaction. However, polarization insertions are neglected in this approximation. This 
will be of some importance for the relativistic generalization, where pions (in general mesons), nucleons 
and resonances appear (see below). Unfortunately the coupled equations (|2.9p, Q2.11D and fl2.12| ) form a 
very complicated set of integral equations in an 8-dimensional coordinate space. They might be visualized 



in a comprehensive way by Feynman diagrams in figure \L1 



N — 

r 



A 



V 



+ 



A 

i) 1 " 



V 



V 2' 



+ 




V 2' 



\ 



1 2, 



/ 



V 2' 





+ 1" 



2" 



T 2' N 



Figure 2.1: Feynman diagrams for the T-matrix approximation 

Here the double (single) line denotes the interacting (non-interacting) Green's function in the T-matrix 
approximation, while the wiggled line is the bare nucleon-nucleon interaction. X and T represent the proper 
self-energy (two point function) and effective interaction or T-matrix (four point function), respectively. 
Even in equilibrium nuclear matter calculations one usually has to make further approximations to come 
to solvable equations. E.g. one approximates the two interacting propagators in the T-matrix equation 
neglecting the hole contributions and/or one treats the self-consistency only on the average etc.. Then 
one obtains e.g. the Brueckner-Hartree-Fock equations (in equilibrium!) or in the relativistic version the 
Dirac-Brueckner equation (see e.g. [flop . 



One important requirement for all kinds of approximations for the dynamical problem is the preserva- 
tion of the basic conservation laws of the underlying theory. It has been shown that the above mentioned 
T-matrix approximation leads to energy, momentum and particle number conserving transport equations 
43]. Furthermore it has also been shown loc. cit. that the T-matrix approximation follows from boundary 



conditions for the two body Green's function which correspond to Boltzmann's assumption of molecular 
chaos. This means that the two body Green's function should factorize for times in the infinite past, 
representing two non-interacting free particles. In order to come to solvable approximations one can e.g. 
approximate the T-matrix by the bare interaction alone (T w V). This leads to the time dependent 
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Hartree-Fock approach which will be discussed below. In this approach the particles interact only via a 
time dependent mean field generated by all particles. Direct interactions (correlations) are neglected (see 
graphs in figure |2T2| ) . 

1 



(D = i — p + 2 i 




2' 



r 



Figure 2.2: Feynman diagrams for the Hartree-Fock approximation 

The self-consistency is still preserved, since the full Green's function in Hartree-Fock approximation 
appears. One can give up the self-consistency to some extent but include higher orders in the interactions. 
The Born approximation is obtained via the graphs in figure |2.3|, which leads to the self-energy as sketched 



in figure 2.4. 




+ 



Figure 2.3: Feynman diagrams for the Born approximation 





Figure 2.4: Feynman diagrams for the self-energy in the Born approximation 

However, in most cases the self consistent propagator has to be replaced by some further approximations 
for practical reasons. One can use e.g. a propagator calculated for a local equilibrium situation, i.e. 
replacing essentially the bare mass of the particle (here nucleons) by an effective mass or effective dispersion 
relation, respectively, which is calculated for an equilibrium system. Usually this is called dressing the 
propagator. In this way one can take also polarization insertions into account, which are important for 
the relativistic generalization. An analogous approximation can also be applied to the self-energy itself. 
I.e. in the dynamical theory one uses the self-energy calculated from equilibrium theory (in T-matrix 
approximation) j|7|]. In this approach self-consistency is still preserved in the equilibrium calculation of 
the effective interaction, but then this interaction is also used for the non equilibrium situation. 

Only a few attempts have been made to take the non equilibrium situation at the beginning of a high 
energy nucleus nucleus collision at least in a stationary approximation into account, e.g. for two streams of 
nucleons inter-penetrating each other f£8 |. It should be kept in mind, however, that in most works the 
approximations mentioned above are combined with a gradient and quasi particle approximation which 
then leads to transport equations for the dynamical problem. In order to summarize the necessary steps 
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to derive a transport theory, we do not need to specify a certain approximation scheme for the self-energy, 
but rather take this quantity as calculable and therefore known. 
We introduce the retarded and advanced Green's functions by 



g 

9 



9^-9 
9 C -9' 



9^~9 U 
9 < -9° 



(2.13) 

and analogous equations for the retarded and advanced self energies. The advantage of these quantities 
is that the real part of the retarded self-energy can be interpreted as the (in general non local) mean 
field while the imaginary part describes the absorption or finite life time of the quasi particles (dressed 
nucleons) . 

If the (anti) chronological quantities are eliminated from equation (|2.8|) one obtains in obvious operator 
notation for the two point functions 



± 



1 



> 



(2.14) 
(2.15) 



If one adds or subtracts the corresponding hermitian conjugate equations one can separate the real and 
imaginary parts of these equations, e.g. 



(8(1) 



g 



ReY> + , g< 
+ ^P > ,0 < ]. 



1 



E < , Re g 
P< , g>] + 



(2.16) 



{s(l) + s*(l'))g 



9 



+ 



> 



, Reg' 



(2.17) 



where the brackets, [A, B]±, denote (anti) commutators of the operators. 

Those equations for two point functions A(l, 1') are now transformed into phase space by the Wigner 
transform defined by 

l l \ 

y) . (2.18) 



A w (R,p) = J<p yf *>-y/ h A(R+±y,R-± 



The main advantage of this transformation is that functions sharply peaked with respect to the relative 
variable y and smooth with respect to R translate into smooth phase space functions with respect to both 
arguments, R and p. In this way local operators have momentum independent Wigner transforms etc.. 
To perform the Wigner transformations, equations (|2.16|) and (|2.17| ) are rewritten in terms of the new 



variables R and y and the operation / d 4 j/exp(ip ■ y/h) is then applied to the equations. Derivatives with 
respect to y can be expressed by multiplications with p by partial integration. 

d 



ol m 



2m 4 2m R 



> 










9w = 


Re S + , g< 


+ 


£< , Reg + 






-,w 




-,w 



,w 



2 1 



(2.19) 



> 










9w = 




+ 


S< , Re g + 






+,w 




+,w 



9 K }-w 



-,w 



(2.20) 
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Here the subscript W at the (anti-) commutators denote the Wigner transforms of the whole quantity. 
In order to express the right hand side in terms of Wigner transforms, we need the Wigner transform of 
Operator products, which is in general given by 



(AB) w (R,p) = A w (R,p)ex P r-2\B w (R,p) 



(2.21) 



The differential operator, A, in the exponent is given by d P • Or — Or • d P in four vector notation and 
where the arrows indicate the direction into which the partial derivatives act. 

The essential simplification is now to neglect higher than first order derivatives. More precisely: we 
require quantities of the form 



\% d p A w ■ d R B w 
2 



A W B W 



< 1 



and 



ifi d R A w ■ d p B w 
2 



AwBw 



<C 1 



(2.22) 



to be small every where in phase space i.e. that the operators are sufficiently diagonal (local) in space-time 
representation and slowly varying with R = (x + x')/2, and that higher order derivatives can be omitted. 
It should be mentioned that the smoothness of the Wigner transformed quantities in phase space can be 
enforced to some extent by averaging them over a sufficiently large region in phase space (of volume larger 
than h), what usually is called coarse graining. Of course in doing this some of the information contained 
in the exact Wigner transforms is lost, but one gains in simplicity of the dynamical description and in the 
case of heavy ion collisions this is the only theoretical approximation which leads to (numerically) solvable 
dynamical equations. Other approaches like the Quantum Molecular Dynamics model are motivated by 
the form of these equations but up to now these models are lacking a strict theoretical derivation. Now 
in a gradient expansion up to first order derivatives the Wigner transforms of the (anti) commutators, 
[A , B]± t w, can be expressed by simple products or Poisson brackets of the Wigner transformed operators 
by 

[A , B} +w = 2 A w B w (2.23) 
[A , B}_ w = ih{A w , B w } = ih(d p A w ■ d R B w - d R A w ■ d p B w ) (2.24) 

Furthermore the left hand side of Eq. fl2.19| ) can also be expressed by a Poisson bracket namely ih{u — 

> 

i 9w }' where u denotes the zero component of the momentum four vector, while the derivative on the 



2m 



left hand side of Eq. Q2.20 ) is of second order and consequently should be neglected in a gradient expansion 
up to first order. We arrive at 



OJ 



El 

2m 



ReZ + , g 



ih 



> 

E< 



Reg- 



>„< 



Z>g 



^2.25) 



— Re S" 

2m 



g < - S < Reg- 



4 1 



g K } - 



(2.26) 



where we dropped the index W for the Wigner transformed quantities. These are of course just c-number 
functions in phase space and the above equations of motion are therefore local, in contrast to the exact 
equations. The analogous equations for the retarded Green's functions are 




^ ± ,9 ± 



0. 



^2.27) 
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p2 

U) 

2m 



g ± = l. (2.26 



The right hand side of Eq. ( p. 25 ) is the collision term, with the rate of loss minus the rate of gain (note 
an over all minus sign implicitly contained in the Green's functions). 

Especially the self energies S > ' < can be interpreted as the rate of scattering out (in) per occupied 
(empty) state respectively. The correlation functions g > ' < are connected to the local occupation number, 
f(R,p), and the local spectral function, a(R,p), via the definitions of a and /, 

9 > = ^<l-f), g <=lf. (2.29) 

in in 

such that a = ih(g > — g<) = ih(g + — g~). The interpretation of a and / is motivated by the solutions for 

2 

g > ' < in a non interacting equilibrium system of fermions, where a(R,p) ~ 6(u — (independent of R) 
and f(R,p) is the Fermi distribution. In this case the energy of the particles is fixed to their mass shell 
value. On the left hand side of (|2.25| ), the first Poisson bracket describes the motion of the phase space 
points propagated by a semi-classical Hamiltonian given by p 2 /2m + 3ft S + . The real part of the retarded 
self-energy can therefore be interpreted as a space and momentum dependent potential. The difference to 
a completely classical motion is of course again the possibility of energies off the mass shell. The second 
Poisson bracket on the left hand side of ( 2.25 ) reflects the kinetic effects of the collisions, e.g. the effects 



of the energy and momentum change in collisions on the effective mass etc. In contrast to the transport 
equations (|2.25[) , equations ( |2.26| ) which should be simultaneously fulfilled are not studied in great detail. 

2 

The reason is that the possibility of off shell propagation {u ^ j- + 3ft£ + ) in between collisions makes 
practical solutions rather difficult. Therefore one usually neglects all terms except the first one in equation 
Q2.26|) , which is called the quasi particle approximation. In this approximation the energy of the particles 
is fixed to their mass shell value and therefore the energy is no independent variable anymore. In order to 
see when this is justified we solve equation ( |2.28|) for the retarded and advanced Green's function 



0* = g ^ . ■ • ( 2 -30) 

2m 



io - - £± ± ie 



The ie prescription is only needed if the imaginary part of the self energy vanishes or is neglected. If 
this solution is put into fl2.27| ) the latter equation is automatically fulfilled. Subtracting the two equations 
( f2.27| ) from each other and doing the same with ( |2.28| ) we obtain two equations for the spectral function: 

{u-E,a} = {T,^g + }, (2.31) 

(oj-E)a = T$tg + . (2.32) 

Here we introduced 

T = m(£> - £<) = ift(£+ - IT) = -2h% £+ (2.33) 

2 

in analogy to the spectral function and the abbreviation, E, for E(R,p) = + 3ftS + (i?, p). Observing 



that from Eq. ( 2.30|) it follows that 



^9 + = ^— 2 = 1 ? , (2-34) 

^-EY + {0 (-"£) 2 + 7 2 ' 



with 7 = T/2h, one can solve the algebraic equation ( |2.32| ) for a with the result 



r 7 

^ = 2h- -. (2.35) 

{u-Ef + {0 (-"£) 2 + 7 2 1 ) 
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The spectral function is therefore (formally) a Lorentzian with the width 7 = -3E + around the mass 

2 

shell E = ^— + 9ft S + . It can easily be shown that the solution ( |2.35| ) of ( |2.32| ) is also a solution of ( |2.31| ). 
The quasi particle approximation is valid if the Lorentzian becomes a 5-function in the limit Y/E <C 1. 
The imaginary part of the retarded self-energy varies from —10 MeV at E = 125 MeV to —25 MeV at 
E = 325 MeV so that the fraction Y/E seems to be sufficiently small for nonrelativistic energies. However, 
it should be kept in mind that off shell propagation can cause changes in the order of more than 10% in 
the scattering rates since in the collision term foldings of four spectral functions occur, which increases 
the effective width by a factor of 2 — 4. However, if we adopt the quasi particle approximation and if 
we neglect the momentum dependence of the real part of the retarded self-energy in the collision term, 
one immediately obtains the so called Vlasov- Uehling-Uhlenbeck equation (VUU, sometimes also called 
Boltzmann-Uehling-Uhlenbeck BUU) : 



dt 



m 



2g 



V ri - V n 9ft£ + ■ V 



pi 



fi(ri,Pi,t) 



J d 3 p 2 J d 3 p[ J dV 2 5\ Pl +p 2 - pi - p' 2 ) 



m 2 (2nhf 

x mm - /i)(i - h) - hh{\ - - m 



da 

dn 



(2.36) 



which especially in its relativistic version (see below) is a starting point for many practical simulations 
of heavy ion collisions (g = 4 denotes the nucleon degeneracy factor). If the r.h.s of Eq. ( |2.36| ), the 
collision integral, is neglected one obtains the Vlasov equation |4i|. On the other hand, if the real part 
of the retarded self-energy 9ft£ + is neglected and if the so called Pauli blocking factors (1 — /) in the 
collision integral are approximated by 1, one obtains the famous Boltzmann equation. The modification 
of Boltzmann's collision integral by the Pauli blocking factors, which take into account the Pauli principle 
for the final state of the two-body collisions, was first introduced by Nordheim |5(J, and latter worked 
out by Uehling and Uhlenbeck [51, 52]. Examples for the application of Eq. ( j2.36j ) to heavy ion collisions 
can be found e.g. in [Q |54], ^5], |56|, 0, |5"8| . It should be noted that the momentum dependence of the 
self-energy causes some practical problems in the collision term, since scatterings change the momenta 
of the particles and therefore also the potential energy. In principle this is correctly incorporated in 
equation ( |2.25|) where the energy E = p 2 /2m + 9ft£ + occurs in the spectral function (5-function in quasi 
particle limit) and from this one can derive more general transport equations than ( |2.36[ ). However, then 
additional factors (Z-factors) appear in the collision integral which make the practical solution more difficult 
¥1, Such collision terms are standard in other fields [50, 51]. In equation ( |2.36 ), only the free on shell 
energy p° = p 2 /2m appears in the energy conserving 5- function. Therefore in some practical computer 
simulations the momentum or the energy respectively must be rescaled to guarantee the conservation laws 
if the potential is momentum dependent. 

The basic ingredients of transport equations like (|2.36|) are the real part of the retarded self-energy (the 
potential energy) and the cross sections which both should be calculated consistently from the same basic 
interaction e.g. in G-matrix approximation, etc. Then the cross section and the G-matrix are connected 
by 



(p lP2 \G\ P ' lP ' 2 ) a - ; 



±9 



(2tt) 5 ^ 7 da 
2 V 3 dn 



ni 



(Pi -£>i)5 3 (Pi +p 2 -P\ 



P2) 



(2.37) 



where the plane waves are normalized to the volume V |47| . 

In practice one solves equations like ( ]2.36|) by Monte-Carlo methods. First practical solutions of the 
Boltzmann equation (no self-energies, no Pauli blocking) were done in the late 70's for hadron induced and 
heavy-ion reactions p 



3]. |64|. |65]. |66]. |67|] . Later, the BUU simulations with test-particle methods were the 
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first to include mean potentials and Pauli blocking [^3], |55], [56], ^ . The phase space distributions can 
be represented by the density of test particles. The latter are propagated classically with the Hamiltonian 
E(R,p) = p 2 /2m + 3ft£ + (i£, p) which yields a solution of the Vlasov part (l.h.s. = 0) of the transport 
equation. At the point of closest approach of two particles one performs a random scattering process, if 
the distance at that point is smaller than \Jatot/^- The angular distribution for the scattering is obtained 
from da/dQ. 

Up to now there is no stringent proof that this procedure really yields a solution of the full transport 
equation, but it can be shown numerically that this prescription yields the correct asymptotic approach to 
equilibrium, which can be deduced from the Boltzmann equation, i.e. equation Q2.36Q , without potential 
and without the (1 — /) Pauli-blocking factors |68|| . 

However, a deficiency of all transport equations for the single particle phase space distribution is that 
clusters or fragments in heavy ion collisions can only be constructed in a statistical way, e.g. with the 
help of a coalescence model. The reason is that in transport theory bound states of many nucleons cannot 
be described directly. This can be avoided if one combines the fully classical equations of motion with 
the random scatterings of transport theory. The classical many body theory with two body forces allows 
for bound states of the particles but completely ignores any quantum effects. Especially the hard core 
of the nucleon-nucleon interaction would lead to a completely different classical behavior as compared to 
quantum mechanical scattering. Therefore in models like QMD, RQMD etc. one translates the soft (long 
range) part of the effective interaction in the medium (ReS + ) into a classical but density and possibly 
momentum dependent two body interaction. However this must be supplemented by random scatterings 
to account for the short range part of the interaction and the other quantum effects like Pauli blocking. 
Unfortunately this procedure can only in part be justified from first principles f69fl , but it can be shown 
that this model yields quite similar results as the transport model as long as one body observables are 
concerned 



70, ITT], [72] . They are different though for many body observables like fragments. 



2.2 Relativist ic transport theory 



We now proceed to generalize the approximation scheme outlined above to the relativistic case. Here one 
should start from an effective Lagrangian with baryons and mesons as elementary degrees of freedom. 
Since the well known o — uo— model has been proven to describe static properties of finite nuclei very well 



73] , |74| it seems to be reasonable to start with this model also for the description of the dynamics of heavy 



ion collisions. It is known that the linear energy dependence of the real part of the self-energy is too strong 
as compared to the measured logarithmic behavior. However, this can be remedied by an effective (energy 
and density dependent) coupling strength [fTER . 

We start from a Lagrangian density for nucleons, deltas and a nucleon resonance like the A^* 440 plus 
the a and to and it meson, for details cf. (76 1. The free part reads: 



+<Mi7/^ - M A ]r A + \d^u - U(a) 



U{uo) + Ud^d^-ml^ 



(2.38) 



and U(a), U(c<j) are the self-interaction part of the scalar field [77] and vector field [78 

U{a) = \rn\o 2 + lb(g° NN a) 3 + \c{g° NN a)\ 



(2.39) 
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U(u) 



-m>^(l + 



2 Z 2 ' 



(2.40) 



respectively. Cj is the interaction Lagrangian density 



£>i — £nn + £n*n* + £aa + £nn* + £atv + £an* 

= g a NN^( X )^( X ) a ( X ) ~ &N$( X )llSP( X ) U ' l ( x ) + gNN^( X hvl5T • ^(x)^7r(x) 

+g^^*(x)^*(:r)<7(:r) - &. N 4\x)lvrt\x)u*{x) + gW^^h/^r ■ i/;*(x)d^(x) 
+g AA ipAu{x)ip A {x)a{x) - glA^^{xh^ A {x)uj^{x) + gl A $ Av (x)^ 5 T ■ ip A {x)d^ir{x) 
+[g' 7 NN 4*{x)i)(x)a{x) - g^frix^ix^ix) + g* NN 4*(x) W r ■ tP(x)d»iz(x) 
-gl N $ Atl (x)d>*<K(x) ■ S + i>(x) - gl N 4 Afl (x)d^(x) • S + if>*{x) + H.c] 
= g* N 4>{x)T N A ^{x)<S> A {x) + g^ N 4*{x)Tfr{x)^ A {x) + g A AA i, Au {x)T^l{x)^ A {x) 

+[ g * N 4*(x)rT tP(x)$ a (x) -gi N i, Alx {x)d^{x) ■ $+4,{x) 

-gl N 4 Afl (x)d»<K(x) ■ S + r(x) + h. c] (2.41) 

Here, r/>, ip* are the Dirac spinors of the nucleon and iV*(1440) , and ip All is the Rarita-Schwinger spinor 
of the A-baryon. r is the isospin operator of the nucleon and iV* (1440), T is the isospin operator of 
the A, and S + is the isospin transition operator between the isospin 1/2 and 3/2 fields. gJf N = fn/fn-*, 
g AN = /*/m 7r ; = T^* = 7a t a , T a = j A T A , A=cr, lu, ir, the symbols and notation are defined in 
table [21]. 

In the closed time path formalism one obtains quite analogous Dyson equations for the different particle 
species from which one obtains in turn transport equations for nucleons, N* uo and A 12 32 |75|, [76], |75[| . 
However, in the relativistic case, even more difficulties arise as for the nonrelativistic tranport equations. 
T- or G-matrix calculations starting from a Lagrangian like ( |2.41| ) are very cumbersome even for the static 
equilibrium situation. Also the non-locality of the Fock term in the mean-field approximation, leads to 
a non-trivial momentum dependence of the self-energy. The same is true for the Born approximation, 
therefore one usually takes only the Hartree or the mean-field approximation for the Vlasov part in the 
transport equation. Effective masses and momenta enter both the Vlasov part and the collision term 
where the cross section has to be calculated in some approximation. As a qualitative guideline, one can 
calculate cross sections from Lagrangians as (|2.41| ) in the Born approximation. However, in medium dressed 
propagators should be used for the internal lines, taking into account effective mass and momentum in 
Hartree approximation. For a more rigorous derivation, we refer to the work of Cassing et al. 
How the cross sections are determined in the UrQMD model can be found in section ^.3.2. 



10, 8 



One then usually uses the propagators of non-interacting systems of particles with effective mass, e.g. 



for the nucleon and the iV* 440 



9 c,a (x,p) 
9 > (x,p) 
g < (x,p) 



($ + m* N : 



±1 



+ 



7T1 



p2 _ m *2^ -4- j e E{j>) 



■6{p - E(p))f N *(x,p) 



7T1 



E(p) 



7T1 



5(j>o-E(p))[l-f N .(x,p))y + m* N : 



E(p) 



6(po - E(p))f N *(x,p)(j!+m* N *] 



(2.42) 
(2.43) 
(2.44) 



where E(p) = Jp 2 + M 2 ^*. The effective mass and momentum of e.g. 



the iV* 440 are defined as 



p^(x) 



(i)=M w , + 4»(x) 



(2.45) 
(2.46) 
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A 


m A 


Snn 




Saa 




1A 




T A 


$4(2) 


U A 


U A 


a 


m a 


Snn 


cr a 

&N*N* 


Saa 


cr " 

SNN* 


1 


1 


1 




1 


1 


oj 


m w 


— rr w 

SAW 


rr^ 

&N*N* 


Saa 


&NN* 




1 


1 






1 


71 


m^ 


SiViV 


&N*N* 


Saa 


&NN* 


^75 


T 


T 


7T(x) 


1 





Table 2.1: Symbol and notation definitions. 

The self-energies are taken here in Hartree Fock approximation only. Herewith one assumes that the 
contribution of the Born diagrams to the real part of the self-energy is small. Furthermore one often also 
neglects the Fock part of the self-energy for the same reason. The self energies in Hartree approximation 
for the iV*^ for example are given by 



Z S n*(x) = -^f[g NN Ps(N)+g%, N ,p s (N*)+gl A Ps(A)], (2.47) 



E^(rr) = ^[g NN ff v (N) + g N * N *f%(N*) + SaaP^(A)]. (2.48) 



m 2 



where ps{i) and Py{i) are the scalar and vector densities of particle species i and analogous expressions for 
iV and Ai232- Furthermore one usually makes use of the so called local density approximation. This means 
that the mean a and uj fields are determined from the local scalar and vector densities alone, neglecting 
retarded contributions from other space time regions, i.e. in the Klein-Gordon equation for the a (and 
analogous for the u): 

(□ + ml )a + b(g% N ) 3 a 2 + c{g NN f a 3 = g% N p s (N) + g N * N *p s (N*) + <7 AA p S (A) (2.49) 

One neglects the dynamical contributions of the D-operator. This should be a good approximation as 
long as there are no rapid field oscillations on the scale of the (small) inverse mass of a and u meson. 
However, the retardation of the interaction should not be neglected for all energies from non-relativistic up 
to the ultra-relativistic energy regime. It is not clear up to now if one should take related phenomena like 
a, oj Bremsstrahlung seriously |HJ] . Moreover, it has been shown that for intermediate energies of several 
hundred MeV per nucleon the above approximation is acceptable |82 . 



As an example, the RBUU equation e.g. for the N* M0 distribution function now reads: 

0v[<9£ - d^ N ,{x)dl + d-Bu^m + m^Eu^m ^^ 

h N * (p) 

= C N *(x,p). (2.50) 

and analogous equations for the nucleon and the Ai232- A practical solution of the relativistic Vlasov 
equation (l.h.s. of equation (|2.50| ) = 0) can be found in HS3| . Attempts have been made also to solve the 



corresponding quantum mechanical problem in terms of the time-dependent Dirac equation with mean 



fields [84]. Realizations of the full relativistic Boltzmann Uehling Uhlenbeck model (RBUU) can be found 
in 



The collision term on the r.h.s. can be written as 



1 f d 3 p 2 f d 3 p 3 r d 3 p4 
(2^)3 J J^f 1 (>) 
xW Nt (p,p 2 ,p 3 ,p 4 )(F 2 -F 1 ), (2.51) 
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where F 2 , F\ are the Uehling-Uhlenbeck factors 



F 2 = [1 - f N * (x, P, t)] [1 - Jb 2 (x, p 2 , t)]f Ba (x, p 3 , t)/ fl4 (x, p 4 , t), (2.52) 

F a = /jv*(x,p,t)/fl 2 (x,p 2 ,t)[l - /fl 3 (x,p3,t)][l - /fl 4 (x,p 4 ,t)], (2.53) 

B 2 , B 3 , B4 can be N, A and iV*(1440) and possibly many other resonances. W N * (p, p 2 , pz, Pi) is the 
transition probability of different channels, which has the form (for details see ref. |76| ) 



w "' iP - P3 ■ n) = i6^.( rt ^ fe )^fe)^fa) ¥J TDtD - TEtB) +P3 ^ Pa - (2 ' 54) 

Here To, Te are the isospin matrices and $d, &e are the spin matrices, respectively. D denotes the 
contribution of the direct diagrams and E is that of the exchange diagrams. A, B = a, u, n represent the 
contributions of different mesons. The exchange of p% and p± is only for the case of identical particles in 
the final state. The two-body scattering reactions relevant to the iV*(1440) in the N, A and iV*(1440) 
system are as follows: 

(1) Elastic reactions: 

NN* — ► NN*, AN* — ► AN*, N*N* — ► N*N* . 

(2) Inelastic reactions: 

NN < — ► NN*, NA < — > NN*, AA < — > NN*, 

NN* < — > AN*, NN* < — > N*N*, NN < — > AN*, 

NA < — > AN*, AA < — ► AA^*, N*N* < — ► AN* , 

NN < — ► N*N*, NA < — > N*N*, AA < — ► N*N*. 
In principle one also obtains resonance formation and decay terms which are not given here. From 
this it becomes clear that for higher energies (above several GeV per nucleon), where dozens of baryon 
and meson resonances come into play, the number of reaction channels increases tremendously. Even for 
comparatively low energies higher resonances contribute significantly to (subthreshold) particle production 
(see section 4.2. 1| ). Most of the corresponding coupling constants are not known and instead one directly 



parameterizes the cross sections and potentials and decay widths. Furthermore the resonances have a 
spectral function (mass distribution) of finite width which is neglected in the quasi particle approximation. 
These spectral functions might change (width and the position of the pole) in the hot and dense nuclear 



medium, which has caused much interest during the last years |86|, [87], ^8], ^9], 0, [H] . 

However, we would like to stress that a shift of the resonance mass distribution to lower masses is not 
necessarily connected to chiral symmetry restoration (m* — > 0), since such a shift could also be obtained 
in a model without chiral symmetry in the high energy limit. 

In a pragmatic way one can incorporate the mass distributions of resonances by folding the collision 



terms with appropriate mass distributions of the involved resonances |92j (see also section |3.3.2|) . In the 
actual numerical implementations this is usually taken into account by choosing the actual resonance 
mass (mass of the test-particles) according to a probability distribution, which might depend on the local 
density and temperature. 

However, the basic concept of well established quasi particles becomes meaningless for energies beyond 
several GeV (baryons), where there is no pronounced resonance structure on the corresponding cross section 
anymore. For these energies one must incorporate a mechanism for continuous excitations of mesons and 
baryons. Since for such high energies also the internal quark and gluon degrees of freedom become more 
and more important one usually supplements the actual resonances by continuous string excitations above 
a (model dependent) energy threshold. These excited strings can then fragment into pre-hadrons by 
string fragmentation algorithms, which have been proven to be successful as long as rescatterings of newly 
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produced hadrons can be neglected ||93|| . A well known example is the FRITIOF model f94| , |95| , which 
however does not incorporate a complete phase space picture but rather works only in momentum space so 
that rescatterings are difficult to be included. Perturbative QCD effects like multiple minijets have been 
included in the PYTHIA model 
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In the UrQMD model, to be discussed in more detail in the following chapter, a string excitation and 
fragmentation scheme is combined with the transport theoretical approach. Other models working along 
the same line are the RQMD (Relativistic Quantum Molecular Dynamics) model [|37]] or the HSD (Hadron 
String Dynamics) |97|]. Some other models are based mainly on the dynamics of strings like VENUS 
and QGSM [|17] or partons like the parton cascade model PCM ||99| , [L00|1 or the dual parton model 
XD2J. 
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Since parton cascades are mainly based on the cross sections from perturbative QCD, a basic require- 
ment is that the typical momentum transfer ^[Cp in parton scatterings is larger than ~ 10 GeV so that 
running coupling constants become sufficiently small. Even if this is fulfilled in first collisions, secondary 
interactions are much softer and require a non-perturbative treatment. Therefore these models have to be 
supplemented by hadronization prescriptions like the parton-hadron conversion model HIJING ||103|| and 



prescriptions for rescattering of newly produced hadrons like the hadron cascade HIJET |104 



A common problem for all models involving a transition from the string or parton picture to hadrons 
is the so called formation time. This is the time after which newly produced quarks and anti-quarks can 
be considered to be in a well established hadronic state (meson or baryon) that can interact again with 
other particles. Not much is known from the theoretical side how such a parameter can be calculated or 
how an appropriate model prescription could be derived. For the prescription used in the UrQMD model, 



we refer to section 3.3.5 



Even though models like RQMD or VENUS are impressively fine tuned for CERN SPS data the price 
to be paid are additional mechanisms or model assumptions like string droplets in VENUS [105] or color 
ropes in RQMD |TUB| 
range (e.g 
QMD/IQMD 



Furthermore, many of these models are not designed to work over a broad energy 
the relativistic cascade ARC 

H |ID§, |109l, TO fTTTH or BUU | 



107|| is only designed for AGS energies and models such as 
% |m[ may only be applied to the BEVALAC/SIS 
energy domain). On the other hand the search for the quark gluon plasma in high energetic heavy ion 
collisions will require a systematic study of the excitation functions of observables over the broad energy 
range from SIS (~ 1 AGeV) up to SPS (~ 200 AGeV) or even higher energies in colliders like RHIC 
(Vs = 20 AGeV) or LHC (y/s > 1000 AGeV). 

A general problem of hadron or parton cascades or transport models is the geometric interpretation 
of the cross section in most of the actual model implementations. Even though the underlying theory or 
model is usually Lorentz invariant the practical realization might not be [113]. This is due to the fact 



that the sequence of collisions according to a specific collision criterion can be reference frame dependent. 
The reason is that points of closest approach of several particles in space and time might be interchanged 
in the time sequence after a Lorentz transformation, such that in different reference frames one obtains 
different results for the sequence of scatterings. In the UrQMD model this problem has been minimized by 
the choice of the collision criterion (see section |3.3.1| ) such that the frame dependence of the observables 
(going from the CM to the target rest frame) is only in the order of ~ 5% at CERN/SPS energies. 



2.3 The Quantum Molecular Dynamics approach 



In the following we discuss an approach which goes beyond the one-body descriptions as discussed above, 



which is the Quantum Molecular Dynamics (QMD) model [5S, 108 



110]. The QMD model is a 
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iV-body theory which simulates heavy ion reactions at intermediate energies on an event by event basis. 
Taking into account all fluctuations and correlations has basically two advantages: i) many-body processes, 
in particular the formation of complex fragments are explicitly treated and ii) the model allows for an 
event-by-event analysis of heavy ion reactions similar to the methods which are used for the analysis of 
exclusive high acceptance data. 

The major aspects of the formulation of QMD will now be discussed briefly. For a more detailed 
description we refer to ref. [ 6lj . The particular realization of the UrQMD model will be discussed later. 



2.3.1 Formal derivation of the transport equation 

In QMD each nucleon is represented by a coherent state of the form (we set H, c = 1) 

/ 2 \ 3 / 4 f 2 1 1 

M*;q,P,t)=(-) exp {-j(x - qi (t)) 2 + -i Pi (t)xj (2.55) 

which is characterized by 6 time-dependent parameters, q { and p {1 respectively. The parameter L, which 
is related to the extension of the wave packet in phase space, is fixed. The total n-body wave function is 
assumed to be the direct product of coherent states (|2.55|) 



^ = l[ ( f) l (x,q i ,p i ,t) (2.56) 

i 

Note that we do not use a Slater determinant and thus neglect antisymmetrization. The computational 
time scales like (A p + A t Y in the case of a Slater determinant, while in QMD it is like (A p + A t ) 2 . First 
successful attempts to simulate heavy ion reactions with antisymmetrized states have been performed for 



small systems [|114| , |1 1 5|| , for a recent review see Ref. [ |116|| . A consistent derivation of the QMD equations 
of motion for the wave function under the influence of both, the real and the imaginary part of the G- 
matrix, however, has not yet been achieved. Therefore we will add the imaginary part as a cross section 
and treat them as in the cascade approach. How to incorporate cross sections into an antisymmetrized 
molecular dynamics is not yet known. This limits its applicability to very low beam energies. 

The initial values of the parameters are chosen in a way that the ensemble of At + Ap nucleons gives 
a proper density distribution as well as a proper momentum distribution of projectile and target nuclei. 

The equations of motion of the many-body system is calculated by means of a generalized variational 
principle: we start out from the action 

ti 

S = f £[$,$*]cft (2.57) 



tl 








dt 



with the Lagrange functional C 

$\ (2.58) 

where the total time derivative includes the derivation with respect to the parameters. The Hamiltonian 
H contains a kinetic term and mutual interactions Vy, which can be interpreted as the real part of the 
Brueckner G-matrix supplemented by the Coulomb interaction. We will later on describe the components 
of H in detail. The time evolution of the parameters is obtained by the requirement that the action is 
stationary under the allowed variation of the wave function. This yields an Euler-Lagrange equation for 
each parameter. 

If the true solution of the Schrodinger equation is contained in the restricted set of wave functions 
(pi(x; tjfj, Pi , t) this variation of the action gives the exact solution of the Schrodinger equation. If the 
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parameter space is too restricted we obtain that wave function in the restricted parameter space which 
comes closest to the solution of the Schrodinger equation. Note that the set of wave functions which 
can be covered with special parameterizations is not necessarily a subspace of Hilbert-space, thus the 
superposition principle does not hold. 

For the coherent states and a Hamiltonian of the form H = J2%Ti + \ Vij (Ti= kinetic energy, 
= potential energy) the Lagrangian and the variation can easily be calculated and we obtain: 



2 £M*> 



2Lm 



Pi 



m 



Pi 



(H) 



(2.59) 

(2.60) 
(2.61) 



with (Vij) = J d 3 Xi d 3 x 2 (fi*(fijV(xi,X2)4>i(frj- These are the time evolution equations which are solved 
numerically. Thus the variational principle reduces the time evolution of the n-body Schrodinger equation 
to the time evolution equations of 6 ■ (Ap + At) parameters to which a physical meaning can be attributed. 
The equations of motion for the parameters p { and q { read 



Pi 



d(H) 
9q { 



and q { 



d(H) 
dPi 



(2.62) 



and show the same structure as the classical Hamilton equations. The numerical solution can be treated 
in a similar manner as it is done in classical molecular dynamics El |117|, |118|, |119|, |120||. Using trial 



wave functions other than Gaussians in Eq. ( 2.55| ) yields more complex equations of motion and hence the 
analogy to classical molecular dynamics is lost. If (H) has no explicit time dependence, QMD conserves 
energy and momentum by construction. 



2.3.2 Inclusion of collisions 

As stated above the imaginary part of the G-matrix acts like a collision term. In the QMD simulation 
we restrict ourselves to binary collisions (two-body level). The collisions are performed in a point-particle 
sense in a similar way as in VUU or cascade |5"5| : Two particles collide if their minimum distance d, i.e. 
the minimum relative distance of the centroids of the Gaussians during their motion, in their CM frame 
fulfills the requirement: 

d<d = yj^-, otot = <r(Vs, type). (2.63) 

where the cross section is assumed to be the free cross section of the regarded collision type (N — N, 
N-A,...). 

Beside the parameters describing the N-N potential, the cross sections constitute another major part of 
the model. In principle, both quantities are connected and can be deduced from Brueckner theory. QMD- 
calculations using consistently derived cross-sections and potentials from the local phase space distributions 
have been discussed e.g. in [0. Such simulations are time-consuming since the cross-sections and potentials 
do explicitly depend on the local phase space population. 

Within the framework of using free cross section one may parameterize the cross section of the processes 
to fit to the experimental data if available. For unknown cross sections isospin symmetry and time 
reversibility is assumed. 
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2.3.3 Pauli blocking due to Fermi statistics 



The cross section is reduced to an effective cross section by the Pauli-blocking. For each collision the phase 
space densities in the final states are checked in order to assure that the final distribution in phase space is 
in agreement with the Pauli principle (/ < 1). Phase space in QMD is not discretized into elementary cells 
as in one-body models like VUU. In order to obtain smooth distribution functions the following procedure 
is applied: The phase space density // at the final states 1' and 2' is measured and interpreted as a blocking 
probability. Thus, the collision is only allowed with a probability of (1 — — fy)- If the collision is not 
allowed the particles remain at their original momenta. 



2.3.4 The Relativistic Quantum Molecular Dynamics approach 

The Relativistic Quantum Molecular Dynamics (RQMD) approach has been developed to extend the QMD 
model up to relativistic energies (AGS and CERN/SPS domain) [j37| . Its main improvements compared 
to the standard QMD model j69j [721 |lQ9l TO [TT0| , [TTTl |!2l1 are 

1. covariant dynamics 

2. an improved and extended collision term containing heavy baryon-resonances, strange particles and 
string-excitation for high energy hadron-hadron interactions 

In this section we will focus on the first item - the covariant dynamics - since the description of the 
UrQMD collision term in the following chapter will cover in great detail all the techniques which have also 
been employed in the RQMD model. 

The Relativistic Quantum Molecular Dynamics model describes the time-evolution of a many-body 
system using classical covariant equations of motion. The system propagates in a 8iV-dimensional phase 
space with QN degrees of freedom representing the classical configuration- and momentum-space. The 
remaining 2N degrees of freedom contain the eigentime and energy of each particle. 

The necessity of employing an 8iV-dimensional phase-space is based on the no-interaction theorem 



(NIT) by Curri, Sudarshan and Mukunda |3E]: In a QN + 1-dimensional phase-space a Lagrangian for a 
Poincare-invariant dynamics can only be formulated in the case of non-interacting free particles. In order 
to reduce the 8A^-dimensional phase-space to the commonly used 6N + 1 dimensions, Lorentz-covariant 
constraints have to be introduced. These constraints have to yield reasonable equations of motion in the 
non-relativistic limit of a dilute gas. These constraints impose a time-ordering on the dynamics of the 
system. They consist of N mass-shell constraints: 

^ = pi - m? - Vt = i = 1, . . . , N (2.64) 

and N — 1 constraints which connect the relative times of the particles 

Xi = T,9ijP ij Qi j = Q i = l,...,N (2.65) 

with 

exp(q?/L) 

Qij = Qi~ Qj , Pij =Pi+ Pj , 9a = f (2-66) 

The iVth constraint then serves to attach the individual eigentimes to a common time-evolution parameter 
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In the case of finite range interactions, acausalities may occur. The constraints ( 2.65|) , however, suppress 
correlations of particles with strongly differing eigentimes. 

Furthermore a time-ordering for binary collisions must be imposed (via the time-evolution parameter r) . 
This time-ordering, unfortunately, depends on the reference frame. One has to assume that the number 
of collisions is large enough for the dynamics of a single binary collision not to influence macroscopic 
properties of the system. 

The Hamiltonian is defined as linear combination of all Poincare-invariant constraints: 

N N-l 

H = ^XiHi+^SfMXi (2-67) 

i=i i=i 

The canonical equations of motion are then given as 



with the coefficients A,-: 



ft = _§* = (2.69) 
dr d qj fri dq 3 

Xi^-^Sm i = l,...,N (2.70) 

The A^th constraint xn connects the eigentimes of the individual particles to the time-evolution parameter 
t. Although this parameter is used to describe the time-evolution of the system, it should, however, not 
be interpreted as the system time. The matrix is then given by: 

(S-%- = {^Xi) i,j = l,...,N (2.71) 

Here, the main drawback of the RQMD ansatz appears: In order to solve the equations of motion one 
needs to calculate the coefficients Aj. For their calculation the matrix S* -1 must be inverted. Since the 
number of elements of S' 1 is quadratic in the number of particles N, the inversion of S^ 1 (which can 
only be done numerically) is very time-consuming, especially for heavy collision systems at CERN/SPS 
energies. 
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Chapter 3 

The UrQMD-Model 



3.1 Initialization 

This section describes the initialization of projectile and target nuclei in the UrQMD model. Projectile 
and target are modeled according to the Fermi-gas ansatz. The nucleons are represented by Gaussian 
shaped density distributions: 

<Pj( Xj ,t) = (^) exp j-afo - rj(t)) 2 + -pjtyxA . (3.1) 

The wave-function of the nucleus is defined as the product wave-function of the single nucleon Gaussians: 

$ = \{<f 5 {x h p h t) . (3.2) 

j 

Each initialized nucleus must meet the following constraints: 

• J2iQi — 0, i.e. it is centered in configuration space around 0, 

• J2i v i — 0, i.e. the nucleus is at rest 

• its binding energy should correspond to the value given by the Bethe-Weizsacker formula, 

• the radius should yield the following mass dependence 

R(A) ~ r • (3.3) 

and have a reasonable surface-thickness, 

• in its center, the nucleus should have nuclear matter ground-state density. 

In configuration space the centroids of the Gaussians are randomly distributed within a sphere. The finite 
width of Gaussians results in a surface region beyond the radius of that sphere. Therefore its radius is 
reduced by half a layer of nucleons from the original nuclear radius of equation 



R(A) = r ( - 



a+ (a* - iy 



(3.4) 
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The parameter r is a function of the nuclear matter ground state density p used in the UrQMD model: 



r o = h— • (3-5) 

The relatively small number of nucleons to be distributed over the volume of the nucleus may result in 
large fluctuations in the mean density of the nucleus. Therefore the phase-space density at the location 
of each nucleon is evaluated after its placement. If the phase-space density is too high (i.e. the respective 
area of the nucleus is already occupied by other nucleons), then the location of that nucleon is rejected 
and a new location is randomly chosen. 

The initial momenta of the nucleons are randomly chosen between and the local Thomas-Fermi- 
momentum: 

= hc^pf , (3.6) 

with p being the corresponding local proton- or neutron-density. 

A principal disadvantage of this type of initialization is that the initialized nuclei are not really in 
their ground-state with respect to the Hamiltonian used for their propagation. The parameters of the 



Hamiltonian (see section |3.2|) are tuned to the equation of state of infinite nuclear matter and to properties 
of finite nuclei (such as their binding energy and their root mean square radius). If, however, the energy 
of the nucleons within the nucleus is minimized according to the Hamiltonian in a self-consistent fashion, 
then the nucleus would collapse to a single point in momentum space because the Pauli-principle has not 
been taken into account in the Hamiltonian. 

A viable solution to this problem is the inclusion of fermionic properties of the nucleons via the antisym- 
metrization of the wave-function of the nucleus. This ansatz has first been implemented in the framework 
of the Fermionic Molecular Dynamics (FMD) ||114j| . The FMD equations of motion are computationally 
very expensive - it has therefore not been possible so far, to calculate systems heavier than Ca+Ca in the 



framework of FMD [116 



An alternative recipe is the use of a so-called Pauli-Potential ||118|| in the Hamiltonian. This potential, 
which is repulsive in configuration and momentum space, allows to maintain the product ansatz for the 
wave-function of the nucleus. A self-consistent minimization (e.g. via a Metropolis-algorithm [122]) of 



the energy of the nucleus results in a reasonable ground-state due to the Pauli-Potential mimicking the 
fermionic properties of the nucleons. A drawback of the Pauli-Potential is, however, that the kinetic 
momenta of the nucleons are not anymore equivalent to their canonic momenta, i.e. the nucleons carry 
the correct Fermi-momentum, but their velocity is zero. Furthermore, the Pauli-Potential leads to a 
wrong specific heat and changes the dynamics of fragment formation. A big advantage of the Pauli- 
Potential is that the initialized nuclei remain absolutely stable whereas in the conventional initialization and 
propagation without the Pauli-Potential the nuclei start evaporating single nucleons after approximately 
20 - 30 fm/c. 



3.2 Equations of motion 

This section describes the real part of the nucleon-nucleon interaction as it is implemented into the UrQMD 
model [123]. The interaction is based on a non-relativistic density-dependent Skyrme-type equation of 



state with additional Yukawa- and Coulomb potentials. Momentum dependent potentials are not used 
a Pauli-potential, however, may be included optionally. 



26 



The nucleon- or baryon-density can be obtained from the Gaussian ( |3.1| ) : 

2 ~ 

Qj( x 3>t) = {-^) exp{-2a(^ -r^t)) 2 ). (3.7) 



V 7T 

where Xj denotes the quantum mechanical position variable, while rj(t) is the classical parameter of the 
Gaussian. The Skyrme-Potential (momentum-dependence and spin-exchange has been neglected) has the 
form: 

ySk = E ' s ( x i - x k) + ^2 E ' s ( x i - x k)K x i - x i) , (3-8) 

j,k ' j,k,l 

where in order to avoid self-interactions, all terms where at least two indices are identical are discarded in 
the primed sum. This potential consists of a sum of two- and a three-body interaction terms. The two-body 
term, which has a linear density- dependence models the long range attractive component of the nucleon- 
nucleon interaction, whereas the three-body term with its quadratic density-dependence is responsible for 
the short range repulsive part of the interaction. Using the Gaussian (O) as the wave-function of the 
nucleon we obtain for the two-body Skyrme potential of particle j: 

N 

Vj k2 = J2' dxj dx k ip*{xj) ip* k {x k ) t x 5(xj - x k ) <pj(xj) Pk(xk) 
k J 

= *iE' (-) 2ex p{- a ( r i - r k?} 

= hdtlrs) ■ ( 3 - 9 ) 

In the last line the interaction density was introduced. This density has the same form as the nucleon 
density (|3.7f) obtained from the Wigner-transform of the Gaussian (|3.1|) , but omits the nucleon at the 
location j and its Gaussian has twice the width of that used in equation (|3.7|). The three-body potential 
for particle j can be derived in an analogous fashion: 



l N r 

v ? k3 = ^ E 7 dx i dx k dx i rf( x i) v^^fc) tifa) 

z - kl J 



xt 2 5(xj - x k )5(xj - xi) (fij(xj) (p k (x k ) ipi(xt) (3.10) 
*2^j E ' (^l) ex P {~\ a (( r J - r kf + ^k - n) 2 + (r, - rj) 



In the case of infinite nuclear matter the individual relative distances should be close to their average 
value. Therefore the relative distance between particle k and I may be substituted by the average of the 
other two relative distances: 



t? 3 - ^E'^j'expl-a^-r^ + ^-n) 2 )} . (3.11) 

Using the definition of the interaction-density given in equation ( |3.9|) , the quadratic density dependence of 
the three particle term ( |3.11| ) may be generalized to arbitrary exponents for the density. This is of great 
importance for the implementation of a so-called soft equation of state. Then, however, the interpretation 
of V^ Sk3 as three particle interaction is no longer valid: 

Uf 3 « t 2 3-§ {off - t 7 ( 7 + ir'Hgfy . (3.12) 
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In the UrQMD model expression ( |3.12[ ) is always used, even for the case 7 = 2. 

The Yukawa-, Coulomb- and (optional) Pauli-potentials may be written in the form of two-particle 
interactions: 



= V Uk 1 _ ' " ' X) (3-13) 

\'i ' j I 

^coui = r^-n (3-14) 




(Tj and Tj denote the spin and isospin of particle j and Zj represents its charge. 

In infinite nuclear matter the contribution of the Yukawa-potential to the total energy has a linear 
density-dependence, just like the two-body Skyrme-contribution. Therefore all parameter sets which satisfy 
the following relation for the parameter t\ yield the same equation of state in infinite nuclear matter: 

^1 + 27rK Yuk 7y = const. (3.16) 

In finite nuclei the usage of a Yukawa-potential has the advantage that the parameters can be tuned to 
the proper surface properties of the nuclei without changing the equation of state. 



parameter 


without Pauli-potential 


with Pauli-potential 


a (fur 2 ) 


0.25 


0.1152 


ti (MeV fm 3 ) 


-7264.04 


-84.5 


t 7 (MeV fm 6 ) 


87.65 


188.2 


7 


1.676 


1.46 


y Yuk (MeV fm) 


-0.498 


-85.1 


7y (fm) 


1.4 


1.0 


yPauli ( Me y) 




98.95 


<?o ( fm ) 




2.16 


po (MeV/c) 




120 



Table 3.1: Parameters of the hard equation of state implemented in the UrQMD model, with and without 
Pauli-potential. 

The classical UrQMD Hamiltonian which governs the motion of the parameters r 3 - and of the 
wave-functions is thus given by: 

N 1 if IV 

U rpkin i \ ~* \ ( 7-iSk2 , 771 Yukawa , 7-iCoulomb , rpPaulA 

UrQMD = 1^ &j +0 2^ 2^ \ h jk + h jk + h jk + &jk I 

J =1 1 j=l fe=l 

1 N N N 

+ EEE^f (3-17) 

u 7=1 fc=i 1=1 



2s 



The individual contributions are defined as: 



kin 



with 



E 



7->Sk2 



7-.Sk3 
&jkl 



771 Yukawa 



77'Coulomb 
h 3k 



rriPauli 



p] + m 2 , 

3 

/ , ( - ) exp{-ar 2 fc } 




exp{-a(r 2 fc + rj)} 
f 1 

exp 




erf 



\fotr jk 



— exp 



2 7 y v / a V '"' 
—j= + y/ar jk ] 



^'^' e -erf (Var jfc 



Vr 



Pau 



ar 



x exp 



jk 



Pjk 



2«g 2 + 1 2^ 



and 



Pjk = \Pj-Pk\ ■ 

So far, only a hard equation of state has been implemented into the UrQMD model. 



(3.18) 
(3.19) 

(3.20) 



(3.21) 
(3.22) 



(3.23) 

(3.24) 
The respective 



parameters are listed in table |3J. 



Unfortunately, the generalization of such two body forces to the relativistic region is not a simple task 
and is not incorporated in the present UrQMD model. In principle, the interaction must be mediated by 
fields, which are propagated according to wave equations or one must make use of the so called constraint 
Hamiltonian dynamics 



, Many other models avoid the propagation of fields using local density 
approximation. Above 2 AGeV, we therefore resort to a cluster decomposition in phase space, i.e. potential 
interactions are only enforced for particles with relative momenta smaller than 2 GeV/c. 



3.3 The collision term 

The UrQMD collision term contains 55 different baryon species (including nucleon, delta and hyperon 
resonances with masses up to 2.25 GeV/c 2 ) and 32 different meson species (including strange meson res- 
onances), which are supplemented by their corresponding anti-particle and all isosp in-projected states. 
The baryons and baryon-resonances which can be populated in UrQMD are listed in table [3.2| , the respec- 
tive mesons in table [37|. The states listed can either be produced in string decays, s-channel collisions 
or resonance decays. For excitations with higher masses than 2 GeV/c 2 a string picture is used. Full 
baryon/antibaryon symmetry is included: The number of implemented baryons therefore defines the num- 
ber of antibaryons in the model and the antibaryon-antibaryon interaction is defined via the baryon-baryon 
interaction cross sections. 

Elementary cross sections are fitted to available proton-proton or pion-proton data. Isospin symme- 
try is used when possible in order to reduce the number of individual cross sections which have to be 
parameterized or tabulated. 
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nucleoli delta lambda sigma xi omega 



-^938 


Ai232 


Ani6 


S1192 


"1315 


-^1440 


^1600 


^-1405 


^1385 


"1530 


-^1520 


^1620 


^-1520 


^1660 


"1690 


^1535 


A 1700 


A 1600 


^1670 


"1820 


^1650 


^1900 


A 1670 


S1750 


"1950 


-^1675 


^1905 


A-1690 


^1775 


"2030 


^1680 


A1910 


A18OO 


S1915 




-^1700 


A1920 


A18IO 


S1940 




Nino 


Al930 


Ai820 


^2030 




N\720 


^1950 


A 1830 






-^1900 




Al890 






^1990 




A210O 






^2080 




A2110 






-^2190 










^2200 










^2250 











Table 3.2: Baryons and baryon-resonances included into the UrQMD model. Through baryon-antibaryon 
symmetry the respective antibaryon states are included as well. 



o-+ 


1— 


0++ 


1++ 


1+- 


2++ 




(1 — )** 


7T 


P 


a 




h 


a 2 


P1450 


P1700 


K 


K* 


K* 


K* 




K* 


K* 

-"■1410 


K* 

-"■1680 


V 




fo 


fi 


h 


/ 2 


^1420 


^1662 


7/ 





fo 


fi 


K 


fi 


01680 


01900 



Table 3.3: Mesons and meson-resonances, sorted with respect to spin and parity, included into the 
UrQMD model. 
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3.3.1 The collision criterion 



In the UrQMD model hadron-hadron collisions are performed stochastically, in a similar way as in the 



original cascade models ||65|| . The cross section is interpreted geometrically as an area. Two particles 



collide if their distance d*™™ fulfills the relation: 



dtmns < d = J — , cr tot = a (y/s , type) . (3.25) 

The total cross section cr tot depends on the cm. energy y/s and on the species and quantum numbers of the 
incoming particles. <i tr ans is the relative distance between the two colliding particles (in three dimensional 
configuration space). At the point of closest approach this distance is purely transverse with regard to the 
relative velocity vector of the particles. 

During the calculation each particle is checked at the beginning of each time step whether it will collide 
according to criterion ( |3.25D within that time step. After each binary collision or decay the outgoing 



particles are checked for further collisions within the respective time step. This procedure assumes that 
all particles have the same clock, i.e. that a 6iV + 1 dimensional phase-space is used. 

In relativistic heavy ion collisions, however, the relative distance between the particles depends on the 
reference-frame of the calculation. Therefore the time-order of the binary collisions and their cross sections 
also depend on the reference frame. One possible ansatz to overcome this unphysical frame-dependence 
is the use of a 8N dimensional phase-space []37||, although even here the time-order of the collisions is not 



unique. In 8N dimensional phase space each particle has its own eigentime. The connection between a 
time-evolution parameter r (the system clock) and the individual eigentimes is constructed via 2N — 1 
Lorentz-covariant constraints, which reduce the 8N dimensional phase-space to a 6iV + 1 dimensional 
phase-space. The particles move in Minkowsky-space along 4-dimensional curved trajectories. At the 
beginning of each time step (in units of r) criterion ( |3.25| ) is applied to scan for collisions. Now, however, 
dtranq is defined as the covariant relative distance: 



trans 



\ K ~ r\ V ' & + PmY - (in ~ QmY , (3.26) 

\ V {Pn+PmY J 



using 4- vectors for the locations q v and momenta p v of the particles. 

In the UrQMD model a 6N + 1 dimensional phase-space has been used. We therefore use a different 
ansatz to minimize the frame-dependence of the collisions: The impact parameter of two colliding particles 
- here also called Gorans _ is well defined in the local rest-frame of the two particles. It corresponds to the 
relative distance between the two particles at the time of closest approach. In order to compute the time 
of closest approach of the two particles one always transforms into the local rest-frame of the two particles 
(using a QN + 1 dimensional phase space - this implies that no time-coordinates are transformed). With 
q i being the locations and p' \ the momenta in the local rest-frame one obtains for the squared impact 
parameter d t 2 rans : 

4ans = <? ~ 4 = fa - <? 2 ) 2 - i(qi ' iP J " P ' 2)f , (3.27) 

VP 1 P 2) 

The constraint of always using the local rest-frame of the colliding particles ensures that the cross section 
of the two particles is always calculated in the same fashion and does not depend on the reference frame 
of the nucleus-nucleus collision. 

The time of closest approach r co ii (i.e. the collision time), however, still depends on the reference- frame 
of the nucleus-nucleus reaction. This dependence cannot be avoided since the system clock is linked to 
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that reference-frame. Therefore the time-order of the individual binary collisions strongly varies with the 



respective reference-frame [124]- Using the locations r*j and the momenta p { in the reference frame of the 



nucleus-nucleus collision one obtains for the time of closest approach for the two colliding particles: 

(ri - r 2 ) ■ (p 1 /E 1 - p 2 /E 2 ) 

Tcoll — / TT, 1 „ , 2 • 

(p 1 /E 1 -p 2 /E 2 ) 2 

We have studied the computational frame dependence in the system S+S at 200 GeV/nucleon: using 
the above described algorithm, particle multiplicities and collision numbers vary by less than 3% between 
the laboratory frame and the CM frame. 

The frame-independent definition of the cross section (via the impact parameter in the two-particle 
rest frame) is an important factor in ensuring the approximate reference-frame independence. 

3.3.2 Cross sections 
Nucleon-nucleon interactions 

In UrQMD cross sections are a function of the incoming and outgoing particle types, their isospins and 
their cm. energy. They may either be tabulated, parameterized according to an algebraic function or 
extracted from other cross sections via general principles, such as detailed balance or the additive quark 
model. 



The total and elastic proton-proton and proton-neutron cross sections are well known [|125|| . Since their 
functional dependence on y/s shows at low energies a complicated shape, UrQMD uses a table-lookup for 
those cross sections. Figure EO shows the comparison between measurements for the proton-proton elastic 



and total cross sections and the respective UrQMD table lookup. The same comparison for the proton- 



neutron case is shown in figure |3.2| . At low energies large differences are visible between the proton-proton 
and the proton-neutron cross sections. Therefore the proper treatment of isospin - especially at low 
energies - is of major importance; simple averaged nucleon-nucleon cross sections should not be used. The 
neutron- neutron cross section is treated as equal to the proton-proton cross section (isospin-symmetry). 
In the high energy limit (-^i > 5 GeV) the CERN/HERA parameterization for the proton-proton cross 
section is used ||125| . 



Particle production in UrQMD either takes place via the decay of a meson- or baryon-resonance or 
via a string excitation and fragmentation. Up to incident beam energies of 8-10 GeV/nucleon particle 
production is dominated by resonance decays. Production cross sections for the excitation of individual 



resonances can be calculated in the framework of OPE or OBE models [|126|1 . Regarding the number of 
implemented resonances in UrQMD and considering the limited applicable energy-range for cross sections 
calculated within OPE and OBE models the calculation of all implemented resonance excitation cross 
sections in the framework of these models is not practical. We therefore employ in UrQMD an effective 
parameterization based on simple phase space considerations; free parameters are tuned to experimental 
measurements. The cross section has the general form: 

*t,2^AV~s) ~ (2S 3 + l)(2S, + l) { p^ T ^- 2 \M(m 3 ,m,)\ 2 . (3.29) 

(pi,2> (v^r 

The matrix element \M.{m 3l m.4) | 2 is assumed to have no spin-dependence but may depend on the masses 
of the outgoing particles. The cross section depends also on the momenta of the in- and outgoing particles 
in the two-particle rest-frame (pi j). If the outgoing particles are stable particles with a well-defined mass, 
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Figure 3.1: UrQMD parameterization of the total 
and elastic proton-proton cross section. The data 
has been taken from [ 125| |. A table-lookup has 



been used at low energies to properly describe the 
data. 
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Figure 3.2: UrQMD parameterization of the to- 
tal and elastic proton-neutron cross section. The 
data has been taken from [ 125|| . A table- lookup 
has been used at low energies to properly describe 
the data. Note the large differences between the 
proton-proton and the proton-neutron data which 
emphasize the importance of correct isospin treat- 
ment 



(P3,4) 1S defined as: 



(P3,4(\/s)) = PCMs(Vs) 



— =J (s - (m 3 + w 4 ) 2 ) (s - (m 3 - m 4 ) s 

l\J s 



(3.30) 



((pi.,2) is defined as above exchanging m 3 , m 4 by mi, m 2 , respectively). If, however, 3 or 4 are resonances, 
the width of their mass distribution must be taken into account. Then relation ( |3.30 ) must be modified by 
an integral over the mass distribution A r (m) of the respective resonance (in the case both are resonances): 



(P3,4,(Vs)) = J J Pcms(Vs, m 3 , 777,4) A 3 (m 3 ) A 4 (m 4 ) dm 3 dm 4 . 



(3.31; 



The mass distribution A r (m) is normalized to unity. The upper and lower boundaries for the integration in 
equation (|3.31| ) are determined by the production thresholds for the respective resonances (e.g. in UrQMD 
this threshold is itln + vn^ = 1.07 GeV for non-strange baryon resonances) and the maximum available 
energy (e.g. in the case of pp — > pA the maximum available energy is \/s — m^). For the solution of 
this integral we assume the mass distribution to have the shape of a free Breit-Wigner distribution with 
a mass-dependent width: 



1 



Tim) 



N(m r -m) 2 + r(m) 2 /4 
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lim A r 

r^o 



777 



8{m r 



777 



(3.32) 
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with the normalization constant 



N = / '— , ol dm. (3.33) 

J (m r -m) 2 + r(m) 2 /4 v ; 

— oo 

Alternatively one can also choose a Breit-Wigner distribution with a fixed width, the normalization con- 
stant then has the value N = 2tt. 

In UrQMD, the excitation of non-strange resonances is subdivided into 6 classes: NN — > iVAi232, 
NN -> NN*, NN -> NA*, NN -> A 12 3 2 A 1232 , JVJV -> A i232^V* and A^A" — > Ai 2 3 2 A*. Here the Ai 2 3 2 is 
explicitly listed, whereas higher excitations of the A resonance have been denoted as A*. For each of these 
classes specific assumptions are made with regard to the form of the matrix-element |_M(m 3 , m 4 )| 2 . Free 
parameters are tuned to experimental measurements, when available. We make the following assumptions 
for the matrix elements: 



1. NN -> A^Ai 232 excitation: 



|_M(v^m 3 ,m 4 )| 2 = A — ^ , (3.34) 

((v/i^-mi) +m 2 A T 2 A 

with m A = 1.232 GeV, T A = 0.115 GeV and A = 40000. Note that this form of the matrix element 



has been adjusted to fit the data shown in figure qA 



2. NN -> NN*, NN -> JVA*, JViV -> A 1232 A^ and JViV -> A 1232 A* excitation: 

\M(m 3 ,m 4 )\ 2 = A- ^ , (3.35) 

(m 4 — m 3 ) z (m 4 + m 3 ) z 

with A = 6.3 for JViV -> A" A"*, A = 12 for iVJV -> A^A* and A = 3.5 for JViV -> Ai 232 A^*. In the case 
of A^A" — > Ai 232 A* there are insufficient data available, therefore we use the same matrix element 
and parameters as in the case of A^A" — > Ai 2 3 2 A r *. Since m 3 ^ m 4 is valid for all above cases, the 
matrix element cannot diverge. 

3. NN -> AA excitation: 

\M(m 3 ,m 4 )\ 2 = A , (3.36) 

with A = 2.8. 

E. g. the total cross section for the excitation of a N* resonance can be obtained by summing equa- 
tion (|3.29| ) over all A^* resonances implemented in the UrQMD model. Figure |3.3| shows the total cross 



sections for the excitation of the different resonance classes A^A" — > A^Ai 232 , A^A" — > NN*, NN — > NA*, 
NN — > Ai 232 Ai 232 , NN — > A 1232 N* and AW — > A i232 A*. Also plotted is the cross section for the excita- 
tion of a string which is defined as the difference between total inelastic cross section and the sum of all 
the resonance excitation cross sections: cr s t r ing = cr tot — a e \ — J2 Cj nel . For y/s > 6 GeV the string excitation 
cross section yields the strongest contribution to the total inelastic cross section. 

Figure |3.4| shows the fit of the UrQMD pp — > A^A 1232 cross section to experimental measurements 



127 ]. The measurements refer to the A + + n exit channel and have been rescaled to match the full 
isospin-summed cross section. In the case of the exclusive A 4232 cross section the quality of the data and 
thus also the quality of the resulting fit is very good. For higher resonance excitations this is unfortunately 
no longer the case and additional measurements are needed to clarify the situation. One has to keep in 
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Figure 3.3: Inelastic proton-proton cross section 
in UrQMD, subdivided into the classes NN — > 
iVAi 232 , NN -> NN*, NN -> NA*, NN -> 
Ai232A 123 2, NN -> A 1232 N*, NN -> A 1232 A* 
and A^A" — > N string. For high energies, string- 
excitation is dominant. 



Figure 3.4: UrQMD fit for the exclusive A 1232 
production in proton-proton reactions compared 
to data |27|. 



mind, however, that the experimental extraction of exclusive resonance production cross sections is only 
possible via two- or three-particle correlations (e.g. a pion-nucleon correlation in the case of the A) which 
introduces large systematic errors, especially for broad resonances. 



In figure |3]5| the UrQMD cross sections for the processes pp — > pp\ im -, pp — > PPi^q VV ~^ PPiem an d 



pp — > PP1700 are compared to data ||127|| . One single parameter has been used to fit all four cross sections. 
The data situation is not as good as in the case of the Ai 232 resonance, some ambiguities are visible which 
results in the quality of the fit being not as good as in the previous case. The parameters for the other 
classes are fitted in the same fashion. 

After the species of the excited resonance has been assigned, its mass must be determined: This 
happens stochastically between the minimum allowed mass for the resonance and the maximum available 
free energy according to the mass distribution ( |3.32|) . 



In the UrQMD model, the full decay width T to t(M) of a resonance is defined as the sum of all partial 
decay widths and depends on the mass of the excited resonance: 

N br 

T tot (M) = J2 r iA M ) ■ (3-37) 

br={i,j} 

The partial decay widths Tij(M) for the decay into the exit channel with particles % and j is given by: 

M \ M M R ))) 1 + . 2( «)- 

here Mr denotes the pole mass of the resonance, T 1 ^ its partial decay width into the channel i and j at 
the pole and I the decay angular momentum of the exit channel. All pole masses and partial decay widths 
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Figure 3.5: UrQMD parameterization for exclusive pp* 
describe all available cross section data 



cross sections. Only one parameter was used to 
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at the pole are taken from the Review of Particle Properties [|125|| . Tij(M) is constructed in such a way 



that Tij(Mf>) = T % j( is fulfilled at the pole. In many cases only crude estimates for T l j[ are given in [|125| 



- the partial decay widths must then be fixed by studying exclusive particle production in elementary 
proton-proton and pion-proton reactions. 

All masses, full widths and decay probabilities used in UrQMD are listed in table |3.4| for non-strange 
baryon- resonances, in table |3~5| for single-strange baryon-resonances, in table [3.6| for double-strange baryon- 
resonances, and in tables |3.7| and |3.8| for meson-resonances. Due to experimental uncertainties there 
frequently exists a large interval in which the respective parameters may be tuned. Therefore the given set 
of parameters is not unique. However, the known uncertainties have been exploited to gain an optimized 
fit to the available inclusive and exclusive meson production data. 

An important component in the computation of mass-dependent decay widths is the momentum of the 
decay products in the rest frame of the resonance, {pij(M)), which has been defined in equations (|3.30| ) 
and ( |3.31| ). In order to avoid recursion (the full decay width as defined in (|3.37|) is needed for the evaluation 
of ( |3.31| )), the full and partial decay widths are evaluated at the beginning of the calculation using a mass 
distribution A r (M) with a fixed mass-independent width and then tabulated for further use. 
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Table 3.4: Masses, widths and branching ratios for non-strange baryon-resonances in UrQMD. Masses 
are given in GeV and the widths in MeV. All parameters are within the range given by the Review of 



Particle Properties [125] and have been tuned to exclusive particle production channels. 



Figure |3.6| shows the total and partial decay-widths of the iV* 535 resonance as a function of its mass. 
This resonance is particularly interesting since it dominates the production of the r\ meson at SIS-energies. 
The grey-shaded area represents the experimental uncertainty of the full width at the resonance pole 
125|| . The opening of the Nt] decay-channel at its threshold energy is clearly visible. Figure |3.7] shows 



the respective probabilities for the different decay channels. Here, the grey-shaded area depicts the Breit- 
Wigner mass-distribution of the N* 535 resonance. Obviously the resonance can also be easily populated 
below the r\ production threshold - due to limited phase-space in a heavy-ion reaction the integrated 
decay-probability of a N± 535 into a nucleon and an 7] meson may lie well below the free branching ratio 
given in the Review of Particle Properties ]|125|] , 

Unfortunately, equation ( |3.38| ) cannot be easily extended to include three- or four-body decay channels. 
In order to treat all decay channels on an equal footing in UrQMD, the outgoing particles of a three- or 
four-body decay are combined into two "effective" particles which are used to compute the respective 
partial decay-widths. iV-body phase-space, however, is treated explicitly. 

All resonances decay isotropically in their rest frame. For a two-particle exit channel the momenta are 
given by equation ( |3.30|) . If a resonance is among the outgoing particles, its mass must first be determined 
according to a Breit-Wigner mass-distribution. If the exit channel contains three or four particles, then 
the respective N-body phase-space must be taken into account for their momenta ] 128| . 



37 



resonance 


mass 


width 


NK 


NK* 92 


E-7T 


E*7T 


A7 


Ar? 


Auj 


A-7T 


£»7 


A*7T 


AK 


A* 

"1405 
A* 

"l520 


1.407 
1.520 


50 
16 


0.45 




1.00 
0.43 


0.11 


0.01 














A* 

"1600 


1.600 


150 


0.35 




0.65 


















^1670 


1.670 


35 


0.20 




0.50 






0.30 












A* 

"1690 


1.690 


60 


0.25 




0.45 


0.30 
















A* 

"1800 


1.800 


300 


0.40 


0.20 


0.20 


0.20 
















^1810 


1.810 


150 


0.35 


0.45 


0.15 


0.05 
















^1820 


1.820 


80 


0.73 




0.16 


0.11 
















^1830 


1.830 


95 


0.10 




0.70 


0.20 
















^1890 


1.890 


100 


0.37 


0.21 


0.11 


0.31 
















^-2100 


2.100 


200 


0.35 


0.20 


0.05 


0.30 




0.02 


0.08 










A* 

-"-2110 


2.110 


200 


0.25 


0.45 


0.30 


















^1385 


1.384 


36 






0.12 










0.88 








^1660 


1.660 


100 


0.30 




0.35 










0.35 








^1670 


1.670 


60 


0.15 




0.70 










0.15 








^1750 


1.750 


90 


0.40 




0.05 












0.55 






^1775 


1.775 


120 


0.40 




0.04 


0.10 








0.23 




0.23 




^1915 


1.915 


120 


0.15 




0.40 


0.05 








0.40 








^1940 


1.940 


220 


0.10 


0.15 


0.15 


0.15 








0.15 




0.15 


0.15 


^2030 


2.030 


180 


0.20 


0.04 


0.10 


0.10 








0.20 




0.18 


0.18 



Table 3.5: Masses, widths and branching ratios for single-strange baryon- resonances in UrQMD. Masses 
are given in GeV and the widths in MeV. All parameters are within the range given by the Review of 
Particle Properties |125| and have been tuned to exclusive particle production channels and to the total 
kaon-nucleon cross section. 
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Table 3.6: Masses, widths and branching ratios for double-strange baryon-resonances in UrQMD. Masses 
are given in GeV and the widths in MeV. All parameters are within the range given by the Review of 
Particle Properties ||125|| and have been tuned to exclusive particle production channels. 
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Table 3.7: Masses, widths and branching ratios for meson-resonances in UrQMD, part I. Masses are 
given in GeV and the widths in MeV. All parameters are within the range given by the Review of Particle 
Properties | |125| . Additional branching ratios can be found in table . 
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Table 3.8: Masses, widths and branching ratios for meson-resonances in UrQMD, part II. Masses are 
given in GeV and the widths in MeV. All parameters are within the range given by the Review of Particle 
Properties [125|. 
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Figure 3.6: Total and partial mass-dependent Figure 3.7: Mass-dependent branching ratios 
decay widths of the N* 535 resonance. The grey- for the different exit-channels of a N* 535 decay, 
shaded area marks the size of the error-bar for the The grey-shaded area symbolizes the Breit-Wigner 

mass-distribution of the resonance. The resonance 
may also be populated in a mass-range in which 
the ^-production is strongly suppressed or even en- 
ergetically impossible. 



total decay-widths according to [ |125|| . The opening 
of the Nrj decay-branch above its energetic thresh- 
old is clearly visible. 



The implementation of the resonances and their inelastic production cross sections can be evaluated 
by studying exclusive meson production in proton-proton reactions: Figure |3.8| shows a comparison for 
exclusive ir° production in proton-proton reactions between the UrQMD model and data [ 127| |. Especially 
at low energies close to the production threshold, such exclusive cross sections are very sensitive to the 
partial decay width and pole mass of the resonances from which the respective particle emerges. Above 
~ 5 GeV all inclusive particle production cross sections are dominated by string decays, which, however, 



have no influence on the exclusive p + p — > p + p + ir channel. 



Baryon-antibaryon annihilation 

The physics of baryon-antibaryon annihilation is still a subject of active research and not yet fully un- 



derstood. A number of different models exist - for an overview we refer to [ |129| , |130|| . In UrQMD the 



elementary annihilation cross section is fitted to the available data: Figure [3.9| shows the UrQMD param- 
eterization in comparison to data for the antiproton-proton total, elastic and annihilation cross section 



125, 127\. 



For the annihilation cross section we employ a parameterization first suggested by Koch and Dover 
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a 



.pp 



a 



N S 



A 2 s 



s 



A 2 s 



+ B 



(3.39) 



with (Jq 



120 mb, s = Am 2 N , A = 50 MeV and B= 0.6. 
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Figure 3.8: Exclusive ir° production cross sec- 
tion in proton-proton reactions within the UrQMD 
model (solid line) compared to data (solid squares) 
|127|| . The comparison allows to evaluate the 
implementation of baryon-resonances into the 
UrQMD model. 



antiproton-proton cross section 




Piab (GeV/c) 

Figure 3.9: Parameterization of antiproton- 
proton cross sections in UrQMD compared to data 
||125|| . The "diffractive" cross section is the differ- 
ence between the total cross section and the sum 
of elastic and annihilation cross sections. 



The antiproton-neutron annihilation cross section is very similar to the antiproton-proton annihilation 
cross section ||132|| and is therefore treated identically. 

The total and elastic proton-antiproton cross sections are treated according to the CERN/HERA 
parameterization: 

a(p) = A + B p n + C In 2 (p) + D \n(p) , (3.40) 



with the laboratory-momentum p in GeV/c. The respective parameters are listed in table 3.S. 





A 


B 


n 


C 


D 


Otot 
0"el 


38.4 
10.2 


77.6 
52.7 


-0.64 
-1.16 


0.26 
0.125 


-1.2 
-1.28 



Table 3.9: Parameters for the CERN/HERA parameterization for the total and elastic proton-antiproton 
cross sections. This parameterization is used in UrQMD for momenta pi^ > 5 GeV/c. 



For momenta p\ a b < 5 GeV/c, UrQMD uses another parameterization to obtain a good fit to the data: 

' 75.0 + 43.1p- 1 + 2.6jr 2 -3.9p : 0.3 < p < 5 
271.6 exp(-l.lp 2 ) : p < 0.3 



Cel 



31.6 

78.6 



18.3P" 1 - l.lp- 2 - 3.8p 



0.3 < p < 5 
p < 0.3 



(3.41) 
(3.42) 
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For low lab- moment a the annihilation cross section is dominant. The sum of annihilation and elastic 
cross section, however, is smaller than the total cross section: 



Act = CT tot - CT c i - CT ann (3.43) 

This difference is called "diffractive" cross section in UrQMD, cr diff = Act, and is used to excite one (or 
both) of the collision partners to a resonance or to a string. In the string case the same excitation scheme 
as for proton-proton reactions is used. For high energies the "diffractive" cross section is the dominant 
contribution to the total antiproton-proton cross section. 

The extrapolation from nucleon-antinucleon annihilation to generic baryon-antibaryon annihilation 
poses a large problem since the respective cross sections have not been measured and the overall un- 
certainties in the understanding of the basic annihilation mechanism leave ample room for speculation. 
Due to different masses and quark-contents, modifications to the nucleon-antinucleon cross sections are 
inevitable. 

Assuming a weak t-dependence of the matrix-element, a phase-space correction analogous to that for 
exchange reactions was suggested in ||133|| , which leads to the following rescaling of the nucleon-antinucleon 



cross section: 

0bb Onn • (3-44) 

Pbb 

A disadvantage of this ansatz is that the cross section diverges for low energies. 

The same -y/s- dependence is used in UrQMD for all baryon-antibaryon cross sections (i.e. the pa- 
rameterization of the proton-antiproton cross section). The different quark-content of the colliding (anti- 
)baryons is taken into account by rescaling the cross section with the ratio of the total baryon-baryon and 
nucleon-nucleon cross sections as obtained from the Additive Quark Model (AQM, see also the following 
section): 



CT 



AQM 
BB 



°bb{Vs) = —aqm a ^(Vs) • (3-45) 



°"nn 



By scaling the cross section with the ratio of the respective AQM cross sections, the strangeness content 
of the colliding baryons is taken into account. 

The final state of a baryon-antibaryon annihilation is generated via the formation of two meson-strings. 
The available cm. energy of the reaction is distributed in equal parts to the two strings which decay in the 
rest frame of the reaction. On the quark level this procedure implies the annihilation of a quark- ant iquark 
pair and the reordering of the remaining constituent quarks into newly produced hadrons (additionally 
taking sea-quarks into account). This model for the baryon-antibaryon annihilation thus follows the 
topology of a rearrangement-graph. Another possibility would be the annihilation of a diquark-antidiquark 
with the subsequent fragmentation of a single meson-string. This ansatz, which is being used in RQMD 
1 57] follows the topology of an annihilation-graph. 

Details in the treatment of the baryon-antibaryon annihilation cross section may have a large influence 
on the final yield of antiprotons and antihyperons: If the baryon-antibaryon annihilation cross section is 
parameterized as the proton-antiproton annihilation cross section but then rescaled to equivalent relative 
momenta in the incoming channel, changes in the order of 50% - 300% may occur. At CERN/SPS energies 
the H yield in central Pb+Pb reactions at 160 GeV/nucleon would decrease by a factor of 3. The p and Y 
yields would then be diminished by 50% and 25%, respectively. 
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Additive Quark Model 



For all interactions for which no experimental data exist (e.g. hyperon-baryon resonance scattering) ad hoc 
assumptions must be made. Ignoring all unknown cross sections and setting them to zero would contradict 
the experimentally observed hadron universality, which states that particle production in highly energetic 
jets does not strongly depend on the incoming collision partners (taking all conservation laws into account). 
A relatively simple prescription for obtaining unknown hadron-hadron cross sections is the Additive Quark 
Model: In this model the cross section only depends on the quark-content of the colliding hadrons [|134j| : 

o-tot = 40 • -j • (1 - 0.4xJ) • (1 - 0Ax s 2 ) (3.46) 

a el = 0.039 -att- (3.47) 

Here «m is the number of colliding mesons and x\ is the ratio of strange quarks to non-strange quarks 
in the i-th hadron. Both relations are phenomenological and do not contain any energy- or momentum- 
dependence. For high energies they agree well with experimentally known hadron-hadron cross sections. 
The missing momentum- or energy- dependence, however, leads to a breakdown of this prescription for 
cross sections close to threshold. 



Detailed Balance 

The principle of detailed balance is based on the time-reversal invariance of the matrix element of the 
reaction. It is most commonly found in text books in the form: 

CT /-n = 4-^-/ > ( 3 - 48 ) 

Pf 9f 

with g denoting the spin-isospin degeneracy factors. UrQMD applies the general principle of detailed 
balance to the following two process classes: 

1. Resonant meson- meson and meson-baryon interactions: Each resonance created via a meson-baryon 
or a meson-meson annihilation may again decay into the two hadron species which originally formed 
it. This symmetry is only violated in the case of three- or four-body decays and string fragmentations, 
since N-body collisions with (N> 2) are not implemented in UrQMD. 

2. Resonance-nucleon or resonance-resonance interactions: the excitation of baryon-resonances in 
UrQMD is handled via parameterized cross sections which have been tuned to data. The reverse 
reactions usually have not been measured - here the principle of detailed balance is applied. Inelastic 
baryon-resonance deexcitation is the only method in UrQMD to absorb mesons (which are bound in 
the resonance). Therefore the application of the detailed balance principle is of crucial importance 
for heavy nucleus-nucleus collisions. 

Equation ( |3.48|) , however, is only valid in the case of stable particles with well-defined masses. Since in 
UrQMD detailed balance is applied to reactions involving resonances with finite lifetimes and broad mass 
distributions, equation ( 3.48Q has to be modified accordingly. For the case of one incoming resonance the 



respective modified detailed balance relation has been derived in ||92|| . Here, we generalize this expression 



for up to two resonances in both, the incoming and the outgoing channels. 
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The differential cross section for the reaction (1 , 2) — > (3 , 4) is given by: 

^ = i^i^n • < 3 *» 

here the pi in the 5-function denote four-momenta. The 5-function ensures that the particles are on mass- 
shell, i.e. their masses are well-defined. If the particle, however, has a broad mass distribution, then 
the ^-function must be substituted by the respective mass distribution (including an integration over the 
mass) : 

da- = -dn {[ m - 2 ^ . (3.50) 

12 64vr 2 sp 12 (m - M,) 2 + T 2 /4 2vr V ; 



Incorporating these modifications into equation fl3.48| ) and neglecting a possible mass- dependence of the 
matrix element we obtain: 

dfflJ (pj 2 ) (2Si + l)(25i + 1) . . . ,,,,,, d<r« 

1ST = W) (2S 3 + l)(2g 4+ l) E_0i-*™.ll J«) — . (3.51) 

Here, Si indicates the spin of particle % and the summation of the Clebsch-Gordan-coeflicients is for the 
isospin of the outgoing channel only. For the incoming channel, isospin is treated explicitly. The summation 
limits are given by: 

J_ = maxflji - j 2 \, \ja ~ (3-52) 
J+ = min(ji + j 2 ,h+h) • (3-53) 

The integration over the mass distributions of the resonances in equation ( |3.51| ) has been denoted by the 
brackets (), e.g. 



pIa ^ (PjU) = J J PcMs(Vs,m 3 ,m 4 ) A 3 (m 3 ) A 4 (m 4 ) dm 3 dm. 



4 , 



which is identical to equation ( 3.3 1| ) . Correspondingly, the mass distribution A r (m) is given by equa- 
tion ( ^32|) . 



The most frequent applications of equation (|3.51 ) in UrQMD are the processes A1232 N N N and 
Ai 2 32 A 1232 -> NN. 

Meson-baryon and meson-meson interactions 

Due to the large pion-nucleon cross section at low cm. energies, resonant meson-baryon and meson-meson 
cross sections are among the most important in UrQMD. Up to cm. energies of 2.2 GeV meson-baryon 
and meson-meson interactions in UrQMD are dominated by resonance scattering, i.e. the formation of 
an intermediate resonance. For example, the total meson-baryon cross section for non-strange particles is 
given by 

— 2S + 1 

°m B (Vs) = (3B,m B ,j M ,m M \\J R ,M R ) ,^ (0Q ~7\ 

R=A,N* i Zb B + J-Ji^M + -U 

PCMS {Mr - ^i) 2 + % 
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Figure 3.10: Total ir~p cross section in UrQMD, 
compared to data ||125|| . Up to pi a b = 2 GeV the 
total cross section consists solely of resonance scat- 
tering. For higher momenta elastic scattering and 
string excitation are employed to fill the total cross 
section. 



Figure 3.11: Totalis p cross section in UrQMD. 
The data is from [125]. Resonant K~p scattering 



is only possible via A and £ resonances, whereas 
in the case n~p only non-strange resonances con- 
tribute. 



with the total and partial A/s-dependent decay widths Y tot and T R ^ MB (see equations ( |3.37| ) and ( |3. 3S| )). 
Therefore, the total pion-nucleon cross section depends on all pole masses, widths and branching ratios of 
all N* and A* resonances listed in table WA. 

Figure |3.10| compares the total 7r~p cross section of UrQMD to data ||125|| . Up to pi a b = 2 GeV the total 
cross section is very well described by resonance scattering. For higher momenta the resonance description 
breaks down and the total cross section is described as a superposition of elastic scattering and inelastic 
string excitation. A comparison between the total ir~p and the total K~p cross sections (see figures |3. 10 



and |3.11| ) is of particular interest: In the first case only A^*- and A-resonances contribute to the cross 
section, whereas in the latter case only hyperon resonances may be excited. 

The same basic principles apply in the meson-meson case: Figure |3.12j shows the total 7r + 7r~ cross 
section, which is dominated by the p-resonance. However, in the case of meson-meson scattering a constant 
elastic cross section of 5 mb is added in order to fully reproduce the measured cross section. 

The cross section for an exclusive exit channel M B — > R — > M' B' may be obtained by replacing 
the total width Y tot in equation (|3.54|) by the respective partial decay- width Tr^m'B'- Figure |3.13| shows 
the cross section for the process 7i~p — > r]n compared to data. Such exclusive processes allow for the 
fine-tuning of the resonance parameters listed in table |3.4j , which is very important since the parameters 

^-production is a rather selective probe: Apart from the 
contribute to this cross section. 



given m 
Af* 535 resonance 



125] often are only crude estimates. 



also the N* 650 , the N^ no , the A"* 710 and the N% 080 
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Figure 3.12: Total 7r + 7r cross section in UrQMD Figure 3.13: Exclusive pion-induced 77 production 



compared to data | 127 |. The cross section is dom- 
inated by the p-resonance. 



in UrQMD compared to data [127]: Apart from 



the N£ 535 resonance, also the Ni* 65Q , the A^qq, the 
A-* 710 and the A2 080 contribute to the cross section. 



3.3.3 Angular distributions 

The UrQMD approach uses an analytical expression for the differential cross-section of in-medium NN 
elastic scattering derived from the collision term of the RBUU equation [|135| , [76|] to determine the scattering 
angles between the outgoing partners of all elementary hadron-hadron collisions. Here it is assumed that 
the angular distributions for all relevant two-body processes are similar and can be described approximately 
by the differential cross-section of in-medium NN elastic scattering. This cross section is 

cr NN ^ NN (s, t) = [D(s, t) + E(s, t) + (s, t < — ► it)], (3.55) 

[2ir) z s 

D(s,t) = J S ° Nn) * (t - 4m* 2 ) 2 + -. ig%N \\ 2 (2s 2 + 2st + t 2 - 8m* 2 s + 8m* 4 ) 
2(t — mr c Y (t — m^) z 

- - s (t - S - mi) w + •) + *»> - «)i + 2(t - ST_ ml) (» - ^ 

x(s - Qm z ) - — — 5r(4m - s - t)m H 

(t — m z )(u — m z ) 

u 2 t 2 - 4m* 2 s - 10m* 2 t + 24m* 4 (t + s) 2 - 2m* 2 s + 2m*% 



+ ( 



4(t — m 2 )(u — mlj) 4(i — m 2 )(« — m 2 ) 

\2r 



^NN&NNj 



r 3m* 2 (4m* 2 - s - t)(4m* 2 - t) 3t(t + s)m* 2 

2(t — m. 2 )(w — m 2 ) 2(i — m 2 r )(n — m 2 ^ 
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ff , 2r 3m* 2 (t + g -4m* 2 )(t + g -2m* 2 ) 3m* 2 (t 2 - 2m* 2 t) 
+ ^^H (t-m 2 )(.-m 2 ) + (t-m 2 )(u-m 2 ) J ' ^ 

where the function D represents the contribution of the direct term and E is the exchange term. The 
coupling strengths are g^ N = 6.9, g^ N = 7.54 and g^ = 1.434 and Mandelstam variables are given by: 

s = (P + V2? = [E*(p) + E*(p 2 )} 2 - (p + P2 ) 2 , (3.58) 
t= (p-p 3 ) 2 = ^(s- Am* 2 ) (cos 9- 1), (3.59) 
u = (p-p 4 ) 2 = Am* 2 - s -t, (3.60) 

with 9 denoting the scattering angle in the cm. system. The in-medium single-particle energy is given by 



E*(p) = ^p 2 + m* 2 . (3.61) 

The formula for the differential cross section of in-medium NN elastic scattering can be extended to 
all elementary hadron-hadron collisions if it is scaled by the replacement 



.a 



s -»• s - (ml + m* 2 y + Am* , (3.62) 

where ml and m\ denote the effective masses of the incoming hadrons. Furthermore, we take into account 
the effects stemming from the finite size of hadrons and a part of the short range correlation by introducing a 
phenomenological form factor at each vertex. For the nucleon-nucleon-meson vertex we take the commonly 
used form 



A 2 



NN A 



k\-t 



(3.63) 



Here is the cut-off mass of the meson A. These cut-off masses are A^ =1200 MeV, A u =808 MeV and 
A^ =500 MeV. 

The total energy and the masses of the incoming hadrons serve as an input for calculating the angular 
distribution. 

Since at the current stage UrQMD only uses free cross sections and free on-shell particles the effective 
in-medium quantities E* and m* are replaced by the respective free quantities E and m. 

It is worth to stress again that equation (|3.55|) is only used for the angular distributions of all elementary 
two-body processes but not for the corresponding total cross sections. 



3.3.4 Resonance lifetimes 

Recently, the treatment of resonance lifetimes in microscopic transport calculations has received much 
attention [ |136|| . The standard approach to resonance lifetimes tr in transport calculations is the application 
of Tfi = 1/Tji in conjunction with a Monte-Carlo sampling of the exponential decay law. The total decay 
width of the resonance T R is either taken to be constant or mass dependent (see e.g. figure |3.6p. The 



resulting lifetimes in the case of the A 12 32 resonance can be seen in the left frame of figure |3. 14 . 

However, this treatment is highly questionable. Time delays and resonance lifetimes in the context of 
semi-classical transport calculations can be investigated via the scattering of an incident wave packet on 
a spherical volume. The details of this investigation are presented in the work of Danielewicz and Pratt 
(DP) ||136[ |. Here, only a brief summary with emphasis on its application in a microscopic transport model 
is given. 
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Figure 3.14: Mass-dependent lifetime of the A(1232) resonance as used in transport model calculations. 
The current prescription of the lifetime computed as inverse (optionally mass-dependent) width is shown 
on the left. The r.h.s. displays lifetimes according to Danielewicz and Pratt [|136| . 



As a result of the scattering, both, the forward going and the scattered part of the wave packet suffer 
time delays. In terms of phase shifts they are given by: 

, dtf, . (2f + l)2cos(2ft) dS l 

ATs= dE ' A ^ = D(2/ + l)(l-4sin^))d^ (3 - 64) 

If one is only interested in the thermodynamic quantities of the system, one can define an ergodic 
constraint which relates the time delays to its densities of state. The investigation of this system, which 
in the context of transport model calculations may be thought of as a "thermal" or infinite matter limit 
yields the following results: 

- At> is consistent with the motion of particle through a mean field: 
U = p ■ ReJF(O) (impulse approximation) 

- all interaction effects (to lowest non-trivial order in density) can either be absorbed into a mean field 
or into a modified scattering delay Ar' s 

At' = tttttt; (3-65) 

5 2sin 2 (5 z )d£ v ; 

Transport model calculations dominantly treat resonances as Breit- Winger resonances. In that case the 
form of the phase shift Si is well known and the application of the formalism yields: 

Ats = (E-/J 2 2 +r V 4 (3 - 66) 

A ° R {E ~ Er)2 ~ F2/4 (*fi7) 
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At' s = 1/Tr ("thermal" or infinite matter limit) (3.68) 

During these delay-times, particle identity should change to ensure the correct density of states. However, 
Atf is negative near Er. When treating the forward time delay as a mean field interaction this negative 
time delay corresponds to an acceleration of the respective particle. If, on the other hand, the time delay 
is treated as such with the scattered particle in a resonance state during that time-span, then Arp and 
(in order to remain consistent) At$ must be rewritten in a way that no negative delays occur: 

A UrQMD _ (3 QQ) 



.UrQMD _ a R (E — E 



ATF ~ ^T((E-E R )* + I*/4) (3 - 70) 

These time delays are plotted in the right frame of figure |3.14| . The 1/T(E) divergence of Arp should 
not pose a serious problem since it appears at an energy far off the pole at which the cross section for 
resonance production usually vanishes. It can therefore be omitted through the proper energy dependence 
of the ratio a R /ap. However, the determination of the cross section o> is not straightforward^, since it is 
experimentally not directly observable via 4£. For the application in a transport model one can attempt 
to determine the strength of <jp via the infinite matter limit: an infinite matter calculation employing At$ 
and Arp with their respective cross sections must yield the same result as a calculation with At' s and 
the scattering cross section alone - if the forward cross section has been chosen properly. This approach 
does of course not yield a unique choice for op but provides a range of ap parameterizations which are 
consistent with the well known "thermal" or infinite matter limit. 

Different ansatzes for op are currently under investigation in the framework of UrQMD. Apart from 
this novel approach it is also possible to perform calculations with tr = 1/Tr, both, for fixed and mass- 
dependent widths. 
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3.3.5 String-excitation and -fragmentation 

The strong interactions of hadrons and nuclei are described by Quantum Chromodynamics. Due to 
asymptotic freedom of QCD the coupling constant a s (Q 2 ) becomes small at small distances or large Q 2 . 
It gives a possibility to apply the perturbation theory to the processes with large momentum transfer. On 
the other hand at large distances (~ 1/Aqcd) the coupling constant is not small and non perturbative 
effects, which are responsible for the confinement should be important. This large distance dynamics 
is very essential for understanding the processes with small momentum transfer, which give a dominant 
contribution to the high energy hadronic interaction. 

If one considers the hadrons as bubbles in QCD vacuum liquid with the chromodynamical fields of 
quarks which do not penetrate in the vacuum medium, then the interactions of such bubbles at high 
energies leads to the production of new objects - color tubes, or strings. The interaction between such 
objects can take place with color transfer when quarks from different hadrons are interacting and without 
color transfer when hadrons are interacting only by exchange of the momenta. In both cases it leads 
to increase of energy and at some moment it will be energetically favorable to break the color tube 
by producing gg-pair from the vacuum. This process repeats until many white bubbles - hadrons will 
be produced. According to the uncertainty principle the time needed for production of a hadron with 

^^Note that when following the mean field approach for the forward time delay, the form of the mean field can directly be 
determined from the experimentally accessible forward scattering amplitude .F(O) 
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momentum p is r « p/m 2 , so in the cm. system the last will be produced the fastest hadrons containing 
the spectator quarks (inside-outside cascade). Each produced q and q have small relative momenta in 
their rest frame or small rapidity distances. As a result due to Lorentz invariance of the picture, finally 
produced hadrons at high energies will be uniformly distributed in rapidity and will have limited transverse 
momenta. 

Hereby for high energy reactions we use a similar picture corresponding to the case when two hadrons 
are interacting only by momentum transfer. It imitates the processes with the double and single diffraction. 
In baryon-baryon (meson-meson) interaction strings between quark q v and diquark qq v (antiquark q v ) from 
the same hadron are produced. The hadron-hadron interactions at high energies are simulated in 3 stages. 
According to the cross sections the type of interaction is defined: elastic, inelastic, antibaryon-baryon 
annihilation etc. In the case of inelastic collision with string excitation the kinematical characteristics of 
strings are modeled in the following way: The hadron momentum transfer px is simulated according to a 
Gaussian distribution 

f(p T ) cc -L=eM-P 2 T /v 2 ) (3-71) 
V vrcr^ 

where a = 1.6 GeV/c. The other interacting hadron gets the same momentum transfer but in the opposite 
direction. The excited strings have the continuous mass distribution f(M) oc 1/M 2 with the masses M, 
limited by the total collision energy yfs: Mi + M 2 < y/s. The rest of the a/s is equally distributed between 
the longitudinal momenta of two produced strings: pi\\ = ~P2\\- The energy of the ith string Ei is defined 
by the longitudinal momentum of the string pi\\, momentum transfer and the mass of the string Mf. 

El = p%+p 2 T + Ml 2 = 1,2 (3.72) 
Ei + E 2 = v^- (3.73) 

The longitudinal momenta of the constituent quarks are chosen according to the structure functions of 
hadrons: 

f(x q ) = fo(x q ) a - 1 (l-x q f- 1 (3.74) 

with a = 0.5 and (3 = 2.5 for valence quark in nucleons. 

The transverse distribution of the constituent quarks was generated according to the same Gaussian 
distribution as for the momentum transfer Eq. ( |3. 71| ) . The diquark transverse momentum is equal in 
magnitude, but of opposite direction. 

The second stage of h-h interactions is connected with string fragmentation. UrQMD uses the 
Field- Feynman fragmentation procedure when the strings decay independently from both ends. It includes 
energy, momentum and quantum number conservation laws and the possibility of converting diquarks into 
mesons via diquark breaking. 

Figure 3.15| schematically sketches the fragmentation excitation of a baryon-string: Two quark- 



antiquark pairs are created (a uu and an ss pair). The leading diquark combines with a newly produced 
s quark to form a hyperon, the newly produced s quark combines with a newly produced u quark to 
form a kaon and the newly produced u quark forms together with the leading quark a pion. In order 
to produce baryon-antibaryon pairs in a string-excitation, also the creation of diquark-antidiquark pairs 
from the Dirac-sea must be possible. The suppression factor for the flavors of the quark-antiquark pair 



produced in string decays can be defined by a Schwinger-like formula ||137f] : 

,2\ 



A4| 2 ~exp(-— ) , (3.75) 



51 



s s u u 

(oo o) — ► (oo o) (o •) (• o) 

Figure 3.15: Particle production in a string-excitation. Two quark- ant iquark pairs are spontaneously 
created in the color flux-tube between the constituent diquark and the constituent quark. 



where A4 gives the matrix-element for the production of a quark- ant iquark pair. The string-tension k is 
(k ~ 1 GeV/fm) and m is the mass of the quark-antiquark pair |93|, |138| , |139|| . 

The suppression factors can also be tuned to the production probabilities for certain meson- and baryon- 
species in elementary proton-proton collisions: The s-quark suppression factor is very sensitive with regard 
to kaon-production, whereas the diquark suppression factor plays an important role for the antinucleon 
production. The strange diquark suppression factor can be tuned by studying antihyperon production. 
The standard values for the different suppression factors in UrQMD are: 

u : d : s : qq = 1 : 1 : 0.35 : 0.1 . (3.76) 

In the UrQMD string-excitation scheme it is not only possible to excite the ground-states of baryons 
and mesons, but also their excited states, which have the same quark-content. If for baryons, the quark- 
content does not determine whether the state belongs to the lowest octet or to the lowest decuplet, then 
octet or decuplet are chosen with equal probability. The probabilities for the higher excited multiplets are 
tuned to multiplicities measured in proton-proton reactions. In the case of mesons the multiplet must also 
be determined before a species can be assigned. The probability of choosing a certain multiplet depends 
on the spin of the multiplet and its average mass: 

P JP c~^±±- . (3.77) 
(m)jpc 



The respective values for the meson multiplets implemented in UrQMD are listed in table 3.10 



multiplet J PC 


-+ i— o ++ 1++ 2++ l +_ 1""* l" ** 


probability 


0.102 0.190 0.056 0.124 0.197 0.127 0.110 0.095 



Table 3.10: Excitation probabilities for the different meson-multiplets in UrQMD. The probabilities 
depend on spin and average mass of the multiplet. 



The excitation and fragmentation of a string follows an iterative scheme: 

string hadron + new string , 

i.e. a quark-antiquark (or a diquark-antidiquark) pair is created and placed between leading constituent 
quark-antiquark (or diquark-quark) pair. Then a hadron is formed randomly on one of the end-points of the 
string. The quark content of the hadron determines its species and charge. In case of resonances the mass 
is determined according to a Breit-Wigner distribution. Finally, the energy-fraction of the string which 
is assigned to the newly created hadron must be determined: After the hadron has been stochastically 
assigned a transverse momentum, the fraction of longitudinal momentum transferred from the string to 
the hadron is determined by the fragmentation function. 
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For the following discussion, it is convenient to introduce light-cone variables in configuration space (z^) 
and in momentum space (p^- They depend on space-time (t,z) and energy-momentum (E,p) coordinates 
and are defined as: 

z ± = t ± z and p ± = E±p . (3.78) 
where For the stochastic fragmentation of a string one defines the Lorentz-invariant quantities: 

x ± = Pha^ron ( < X ± < 1) , (3.79) 
Ptotal 

which represent the longitudinal momentum-fraction of the string which is transferred to the new hadron. 
As examples we show a meson string fragmenting from the r.h.s.: 

p + (q qq q) = x + p + (qq) + (p + — x + p + ) qq , (3.80) 

string meson string 

and a baryon string fragmenting from the l.h.s: 

p~ (qqqqq) = x~p~ (qqq) + (p~—x~p~) qq . (3.81) 

string baryon string 

This iterative fragmentation process is repeated until the remaining energy of the string gets too small for 
a further fragmentation. 

The fragmentation function f(x, mt) represents the probability distribution for hadrons with the trans- 
verse mass m t to acquire the longitudinal momentum fraction x from the fragmenting string. One of the 



most common fragmentation functions is the one used in the LUND model [139]: 



f(x,m t ) ~ -(1 - x) a exp (-b-rrif/x) . (3.82) 

x v ' 

In UrQMD, different fragmentation functions are used for leading nucleons and newly produced particles, 
respectively (see figure |3.16| ): 

/ {x-Bf\ 

f(%)mic = exp I ~2Jp — j ^ or l eacmi S nuc l eons (3.83) 

/(^)prod = (1 — x) 2 for produced particles (3.84) 

with A = 0.275 and B = 0.42. The fragmentation function /(x) pro d, used for newly produced particles, 
is the well-known Field- Feynman fragmentation function [|140| , |141| . At string break-up, gg-pairs have 



zero transverse momentum in the string frame, but the transverse momenta of a single quark, p t , and the 
corresponding antiquark, — p t , is taken according to eq. (|3.71| ). 

The main criteria for the selection of the fragmentation functions and the fine-tuning of their parameters 
is that the particle multiplicities and momentum distributions in elementary proton-proton reactions must 
match the data. 



Figure |3.17| shows the Feynman-x (xf = 2piab/y/s) distributions for emitted protons and lambdas in 



proton-proton reactions at Ei a b w 200 GeV. UrQMD fits quite well the data ||142|| , which is important for 
the predictive power with regard to stopping and particle production in nucleus-nucleus collisions. The 
UrQMD estimate of particle multiplicities in comparison to data is shown for proton-proton reactions at 



200 GeV in table |3.11| , and in figure [3.18j at y/s = 27 GeV. The overall agreement is very good, however, 
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Figure 3.16: UrQMD fragmentation functions: 
nucleons containing leading constituent quarks 
have a Gaussian fragmentation function whereas 
a Field-Feynman fragmentation function is used 
for newly produced particles. 



Figure 3.17: Comparison of the Feynman-rc distri- 
bution for lambdas and protons between UrQMD 
and data [ 142|| in proton-proton reactions at ~ 
200 GeV. 



currently the 0-meson production is underestimated by almost a factor of 2 and (anti-)lambda and E° 



production is overestimated by 50%. Also shown in Fig. |3.19| is the excitation function of anti-proton 
multiplicities in inelastic proton-proton collisions. 

The last central issue to be discussed in this section is the definition of the formation time: The 
formation time of a hadron is the time the constituent quarks of the respective hadron need in order 
to bind together and form the hadron and to tunnel out of the vacuum. For composed particles, there 
are two possibilities to define the formation point: i) "constituent" formation point has the coordinates 
of the string break-up points; it) "yo-yo" formation point has the coordinates of the quark trajectories 
intersection. In UrQMD, the "yo-yo" formation time definition is used. The formation time depends on 
the momentum and energy of the created hadron. 

Figure |3.20| schematically shows the fragmentation of a meson string into several mesons. At the 
vertices V n and V n _\ quark- ant iquark pairs are created. The antiquark of vertex V n then forms a new 
meson together with the quark of vertex V n -\. Configuration space (as depicted in figure |3.20|) and 



momentum space are linked by the following relations | 93[| 



Vn ^ [tn— 1 tn 



and for the light-cone variables: 



Pn 
Vn 



K(Z, 



n-1 



-n-1, 



(3.85) 

(3.86) 



(3.87) 
(3.88) 
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particle 


data 


UrQMD 




2.62 ±0.06 


2.47 


7T + 


3.22 ±0.12 


3.07 


7T° 


3.34 ±0.24 


3.01 


K+ 


0.28 ±0.06 


0.23 


K~ 


0.18 ±0.05 


0.13 


K° 
n S 


0.17 ±0.01 


0.18 


A + S° 


0.096 ±0.01 


0.17 


A + S° 


0.013 ±0.004 


0.03 


P 


1.34 ±0.15 


1.32 


V 


0.05 ±0.02 


0.05 



Table 3.11: 
GeV. 



UrQMD particle multiplicities compared to data |143| for proton-proton reactions at 200 



The formation point and time of the hadron is defined as the point in space-time where the quark trajec- 
tories of the quark- ant iquark pair forming the hadron meet for the first time. In figure |3.20| this point is 
denoted as Pf. The formation time is then given by: 

tf,n = 2^-1 + t n + Zn-1 ~ Z n) = \( Z t,-\ + Z n ) ( 3 - 89 ) 

If the momenta of the produced particles are numbered consecutively from right to left, one obtains in 
momentum space: 

tf,n = Y K ^ M + En ~ Pn ~ 2 E ^ j ' ( 3 - 90 ) 

with L = M/k and T = 2L = 2M/k. 

Depending on the hadron-species, formation times are in the order of 1 - 2 fm/c. During the formation 
time, the cross sections of the leading hadrons containing leading constituent quarks are reduced since 
prior to the formation of the hadron only the constituent quark part of the cross section is present: 

a qh w -0Bh for baryons , (3.91) 
3 

Oqqh ~ jjo-Bh for baryons , 

cr qh ~ -&Mh for mesons . 

Newly created hadrons without any leading constituent quarks of the initial hadron have zero interaction 
cross section during their formation time. 
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Figure 3.18: Data JL44], B EH vs - UrQMD 
particle multiplicities in proton-proton reactions at 
a cm. energy of 27 GeV. 



Figure 3.19: Proton-proton excitation function of 
anti-proton multiplicities as obtained in UrQMD. 
The experimental data are taken from [ |147| 



3.3.6 Color fluctuations, color opacity and transparency 

Quantum Chromodynamics (QCD) has important - and so far not considered - implications on high energy 
processes [148 1 : 



• Hadrons consist of configurations of very different spatial size; 

• At high energies, hadronic quark-gluon configurations can be considered frozen. Due to the long 
coherence length at these energies geometrical color optics can be applied; 

• Small object have reduced interactions; 

• Small objects will expand at high energies; 

• Large-size objects can be captured in the central zone of a heavy ion collision. 

These effects have diverse consequences. One of them is color transparency: Suppose a small sized 
object is created; this object does not interact very much due to color screening. At sufficiently high 
energies, it can pass through an entire nucleus without fluctuating into a large configuration. Thus it does 
not interact - the nucleus is transparent. Nucleus-Nucleus collisions provide a tool to investigate the effects 
of color transparency e.g. in the production of J ||149|| . However, the complementary color opacity effect 



(large size configurations) can cause stronger stopping and significant fluctuations in the transverse energy 
for central reactions. 

A first step to investigate these QCD effects within a microscopic transport model can be made by 
using fluctuations in the elementary hadron hadron reactions. Thus one needs the probability -P(cx) for 
different sized configurations. It is convenient to consider moments of the distribution: 

<a°>= J daa°P(a) = 1 , (3.92) 
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Figure 3.20: Schematic picture of a fragmenting meson string in configuration space. The first newly 
produced meson on the r.h.s. forms at the point Pf at which the quark lines first meet. 



< a 1 >= J daa 1 P{a) = a , (3.93) 

: (3.94) 

where a denotes the average cross section, and -P(er) the probability of a certain cross section. The 
second moment < a 2 >can be determined from diffractive dissociation experiments. In addition further 
information can be obtained from QCD, which implies: 

P(a) oca N "~ 2 , (3.95) 

for a — > 0, with N q being the number of valence quarks. Thus, we get for the nucleon and pion distributions: 

-P/v(°") oc o nucleon , (3.96) 
-P/v(cr) oc const, pion . (3.97) 



from these arguments, we can construct P(<j)- As shown in Fig. |3.21| , one gets a broad distribution for 
proton projectiles and an even broader one for pions. 

By exploring these color coherent phenomena, we have an effective new method to investigate QCD 
and the nuclear structure. This helps to disentangle and select distinct properties of interacting hadronic 
systems. This has been so far neglected in microscopic calculations of heavy ion collisions. 



57 



0.03 



0.025 




10 20 30 40 50 60 70 80 90 100 
er(mb) 



Figure 3.21: Cross-section probability for pions P n (a) and nucleons P/v(cr) as extracted from experimental 
data. P^io = 0) is compared with the perturbative QCD prediction. (Figure taken from ||148|| ). 
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Chapter 4 



Probing hot and dense nuclear matter with 
relativistic heavy ion collisions 



4.1 Baryon stopping 

Head-on collisions of two nuclei can be used to create highly excited nuclear matter || |J. In the course 
of the reaction, a part of the respective longitudinal momenta is converted into transverse momentum and 
secondary particles. This causes the creation of a zone of high energy density. Nuclear shock waves have 
been suggested as a primary mechanism of creating high energy densities in such collisions [0, |I|, |150[ . 



The term nuclear stopping power ||151|| characterizes the degree of stopping which an incident nucleon 
suffers when it collides with the nuclear matter of another nucleus. It can be studied by measuring the 
rapidity distribution of the net-baryons present in the reaction. Obviously, the average collision number 
per baryon increases with the mass number of the colliding nuclei. Thus, the heaviest systems available, 
such as Pb+Pb or Au+Au, are best suited for the creation of strongly stopped matter and high energy 
densities. 

The shape of the baryon rapidity distribution can give clear indications on the onset of critical phe- 
nomena: Due to the strong dependence of the baryon rapidity distribution on the baryon-baryon cross 



section [152, |153|, 154 1, a rapid change in the shape of the scaled dN/d(y/y p ) distribution with varying 



incident beam energy is a clear signal for new degrees of freedom entering the reaction (i.e. decon- 
finement) or for phenomena such as critical scattering jL55]. The width of the dN/d(y /y p ) distribution 



for baryons is inversely proportional to their cross section. In [[OxJ it was pointed out that at ener- 
gies E > 200 GeV/nucleon large fluctuations of size in the nucleon wave function, leading to the color- 
transparency or color-opacity effect, respectively, should be observed. 

The degree of stopping can furthermore be used to estimate the achieved energy density in the course 
of the collision. One often assumes that particles produced at y = ycM originate from the central reaction 
zone at z = and the initial proper time r . The rapidity distribution of these produced particles can 
then be used to estimate the initial energy density in the central reaction zone [ffl: 



rriT dN 
t A dy 



(4.1] 

y=ycM 



Here A is the transverse overlapping region area in the collision and the transverse mass of the 
produced particles. The proper production time is estimated to be of the order of 1 fm/c. Estimates for 
the CERN/SPS energy region were in the order of 1 to 10 GeV/fm 3 pfl. 
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Pb+Pb -^(+)-(-)+X,158 GeV/Nukleon 




1 2 3 4 5 6 

Figure 4.1: Rapidity distributions of net-protons 
for Pb+Pb collisions at the SPS (160 AGeV) in 
comparison to preliminary NA49 data ||157|| . 



Pb+Pb -^h +X,158GeV/Nukleon 




1 2 3 4 5 6 



Figure 4.2: Rapidity distributions of negatively 
charged hadrons (i.e. ir~ , K~ and p) for Pb+Pb 
collisions at the SPS (160 AGeV) in comparison to 
preliminary NA49 data ||157|| . 



Hadronic transport models are well able to predict or reproduce the measured rapidity distributions 
over a vast energy range, if baryon and meson rescattering and particle production via string decay 
| |107| , |159| , p.60| , [L61[ | are incorporated. The VENUS model including possible formation and decay of multi- 

A simple "first collision model" without rescattering 
However, it has been shown recently that a 



quark clusters ||162|1 seems to give similar results 
19" 



163| , |164| , |165|| does not suffice to reproduce the data. 



simple geometrical model of nucleus-nucleus collisions (LEXUS) [|166|| describes rather well the rapidity 
and transverse momentum distributions of protons at SPS energies. It is based on linear extrapolations of 
nucleon-nucleon collisions, following the concepts of Glauber ||167|| and the rows-on-rows model of Hufner 
and Knoll [|168|| . LEXUS ignores rescattering of secondary particles or any other collective effect completely. 



This is visible in the resulting numbers of produced (anti-)strange particles which differ strongly from 
experimental values. Similar in spirit are investigations, where the final state of nucleus-nucleus collisions 
is described as a superposition of isotropically decaying fireballs [ 1159 ] . Here, the longitudinal motion is 
extrapolated from p+p data while the transverse motion is due to random walk collisions, fitted to p+A 
reactions. Like in [|166|| the model is designed to serve as a baseline for further analyses. Unfortunately some 
genuine collective effects which arise in AA collisions may counteract themselves, leaving the discussed 
observables unchanged. With respect to the net-baryon distribution, these studies show that gross features 
alone are not necessarily decisive signatures of new physics. 

Let us first explore the stopping behavior of UrQMD in central Pb+Pb collisions at the SPS (160 
AGeV) as depicted in Fig. |4.1| in comparison to recently measured data by the NA49 collaboration [ p.57] . 
Overall good agreement is observed. A large amount of longitudinal momentum has been deposited, 
leading to excessive particle production at central rapidities in line with preliminary data from NA49 [ |157| 
(cf. Fig. [Op. 

Figure [4.3| shows a UrQMD calculation of the proton rapidity distribution for three presently used 
heavy ion accelerators. In all cases, a system as heavy as Au+Au or Pb+Pb exhibits a central pile- 
up at mid-rapidity. However, the physical processes associated are different: The average longitudinal 
momentum loss in the SIS energy regime is mainly due to the creation of transverse momentum whereas 
at the AGS/SPS energies abundant particle production eats up a considerable amount of the incident 
beam energy. The form of the distributions change from a Gaussian at SIS energies to a plateau at AGS 
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Figure 4.3: Rapidity distributions of net-baryons 
for Au+Au collisions at SIS (1 GeV/nucleon), 
AGS (10.6 GeV/nucleon) and CERN/SPS energies 
(160 GeV/nucleon) calculated with the UrQMD 
model (higher curves). All distributions have been 
normalized to the projectile rapidity in the cen- 
ter of mass frame. We have also plotted rapidity 
distributions of strange baryons decaying into A 
(lower curves). 
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Figure 4.4: Formation time dependence of 
stopping power for the reaction S+S at 200 
GeV/nucleon. Solid dots represent experimental 
data from | |158|| . The cross sections of particles 
within their formation time introduce large uncer- 
tainties for their rapidity distribution. 



and finally to a slight two-bump structure at CERN energies. 

At CERN/SPS energies baryon stopping is dominated not only by rescattering effects but also by the 
formation time after hard collisions in which strings are excited. Within their formation time, baryons 
originating from a leading (di)quark interact with (2/3) 1/3 of their free cross section and mesons with 1/2 
of their free cross sections. The influence of the formation time is shown in figure [4.4| for the system S+S 



at 200 GeV/nucleon. The default calculation (including formation time) reproduces the data ||158|| fairly 
well whereas the calculation with zero formation time (dotted line) exhibits total stopping. A calculation 
with zero cross section within the formation time, however, exhibits transparency. 

In order to study this effect more closely, the i/s distributions of elementary baryon-baryon collisions 
are analyzed for Au+Au collisions at AGS and S+S collisions at SPS energies. Figure [15| shows the 
respective distribution for Au+Au. The collision spectrum is dominated by BB collisions with full cross 
sections and exhibits a maximum at low energies. Approximately 20% of the collisions involve a diquark, 
i.e. a baryon originating from a string decay whose cross section is reduced to 2/3 of its full cross section. 

In figure O, the same analysis is performed for S+S at 200 GeV/nucleon. In contrast to the heavy 
system at AGS the collision spectrum exhibits two pronounced peaks dominated by full BB collisions, one 
in the beam energy range and one in the low (thermal) energy range. Approximately 50% of the collisions, 
most of them at intermediate yfs values, involve baryons stemming from string excitations whose cross 
sections are reduced by factors of 2/3 (referred to as diquarks) or 1/3 (referred to as quarks). The peak 
at high a/s values stems from the initial hard collisions whereas the peak at low energies is related to the 
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BB collisions in Au(10.6GeV/c)Au 

800 I 
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Figure 4.5: distribution for baryon baryon 

collisions in a central Au+Au reaction at the AGS. 



BB collisions in S(200GeV/c)S 




s 1/2 [GeV] 

Figure 4.6: -E^m distribution for baryon baryon 
collisions in a central S+S reaction at SPS energies. 



late, thermal stages of the reaction. 

Figure [4.7| shows the importance of fluctuations of the net-baryon distribution in a single event for 
highest energies. Although the average net-baryon dN/ dy is less than 5 at mid-rapidity, regions of high 
positive or negative net-baryon density occur stochastically in a limited region of phase space. 



4.2 Particle production 



4.2.1 Subthreshold particle production from resonance matter 

Subthreshold particle production is kinematically forbidden in first nucleon-nucleon collisions. Therefore 
it provides a good tool to extract information about collective effects and the high density phase. 

As an example, we consider in the following the subthreshold production of antiprotons: The process 
has first been observed experimentally for Si+Si at BEVALAC ||170|| . Experimental studies have also been 
done recently with a 2.0 GeV/u and 1.8 GeV/u Ne beam from the heavy ion synchrotron (SIS) at GSI 

EH]. EI- 

Since the antiproton yields cannot be described in terms of a first collision model 



173) 



various mecha- 



nisms of antibaryon enrichment have been proposed, e.g. hadronic multi-step processes [|131| , |1 74| , |175 
strong mean field effects ||177|| and meson-meson interaction like pp 



pp 



178 



Subthreshold antiproton production has been investigated with respect to a reduction of the effec- 
tive nucleon mass in the framework of RBUU ||1 75 , pL 80|| , QMD ||181|| and RVUU ||182|| . The important 



contribution of high lying baryon resonances without consideration of in-medium masses has also been 
studied f"183| ]. This we will discuss in the following. Note that the threshold behavior of p-production in 
iViV-reactions for — 4m > 1 GeV is not experimentally determined, so any model calculation depends 
on extrapolations. 

Schematically there are three different energy regions in hadron-hadron interactions: the (quasi-)elastic 
scattering, the region of resonance production and formation and the high energy region characterized by 
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Figure 4.7: Event-averaged rapidity density of net-baryons, hyperons and anti-hyperons in collisions of 
Pb + Pb at 6.5 TeV and b = 0, calculated with Fritiof 7.02. For the mid-rapidity range, the net-baryon 
distribution of one single event is also depicted (figure taken from f34|). 

abundant production of particles. In the resonance region the kinetic energy of the reaction partners is 
sufficient to form excited states of the ingoing hadrons which subsequently decay again into the stable 
hadrons. The most prominent examples are 



ttN — ► A 



1232 



2(770) (4.2) 



NN — > NA 



1232 



The formation of resonances influences the particle production rates and their momentum distributions, 
e.g. the pt spectra. The particle chemistry may change, too, distorting signatures from earlier reaction 
stages, e.g. due to reactions like 

irY < — ► Y* < — ► KN 

tttt < — ► / (980) < — ► KK ^ > 

Any of these secondary interactions can be of utmost importance for particle production as discussed in 
|160| , [174 , |175| , |184|1 , the dominant process being the annihilation of produced mesons on baryons which 



leads to the formation of s channel resonances or strings. Those resonances formed are not only responsible 



for strangeness enrichment at the AGS [|184j , [L85|| , but they may be further excited in a subsequent collision 
to a mass larger than 3mjv which allows for NN creation. 

If the projectile nucleon gets excited in the first collision, but does not lose momentum for target 
excitation (which may happen in a subset of cases), a further excitation in a subsequent collision then may 
help to bring the projectile beyond the B creation threshold. 

On the one hand, the larger the mass of the incoming hadrons, the better the chance to overcome the 
B creation threshold in a subsequent collision, because the second step of excitation can be smaller. On 
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Figure 4.8: RQMD calculation of the collision spectrum for the reaction Ni + Ni at 2 GeV/u (min. bias). 
The contributions are classified according to the heaviest baryon involved: N (solid line), A1232 (dashed 
line), iV* 440 (dotted line) and all higher resonances (dashed-dotted line). Shown is the total spectrum (left) 
and the tails above the p-production threshold (right). 



the other hand, the life time goes down with increasing mass reducing the chance to hit another target 
nucleon. Therefore the multi-step process of p production has to be viewed as a complicated interplay 
between excitation and decay processes. 



Calculations in the framework of RQMD ||186|| exhibit a scenario of cumulative scattering of participants 
and secondaries, a multi-step process of successive binary collisions like ttN, pN, rjN*, etc, as well as AA, 
etc. Baryons and mesons get more and more excited and decay or rescatter subsequently. Possibly, the 
invariant mass of a single resonance-resonance collision exceeds the threshold for antiproton production. 
The collision spectra of AA-reactions at 2 GeV/u show that meson-baryon-interactions play a decisive role 
for the excitation of these high masses. We do not find evidence for any antiproton production via the 
meson-meson channel (especially pp) at SIS-energies. 

According to RQMD, BUU and IQMD calculations, the enrichment of nuclear matter with A resonances 
is only about 15% for medium systems at SIS energies |187| , |188| , [L89[| , higher resonances are even less 
populated. Nevertheless, they turn out to be most important for the production of heavy particles, since 
the collisions at high center of mass energy preferentially involve these highly excited baryons. Fig. 4.8 
shows the baryon-baryon collision spectrum of the reaction Ni+Ni at 2 GeV/u. Above a/s = 4 ■ ton only 
1/3 of the 55-reactions are NN-, A^A 1232 - or A 1 232A 12 32-collisions. 50% involve resonances higher than 
A'iW)- They act as energy depots and dominate the p-production. 

In the other microscopic models mentioned above the contribution of high mass resonances are ne- 
glected, but the antiproton production is strongly enhanced due to a reduced effective mass. In particular 
the p annihilation probability differs drastically between the models [|183|| . This question might be resolved 



by the measurements of anti-flow, as was suggested in 
incompatible scenarios ( see also section [4.2.6| ) . 



Thus, different microscopic models give totally 
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4.2.2 Temperatures and single particle spectra 

The hot and dense reaction zone consists of incident nucleons and produced particles. The fireball model 
considers the hadrons as a gas in thermodynamic equilibrium. For temperatures above 50 MeV, Fermi- 
Dirac and Bose- Einstein distribution functions for baryons and mesons may be approximated by a Maxwell- 



Boltzmann distribution ||191| , |192|| , with the temperature T as parameter. Equilibration is thought to be 
visible predominantly in the transverse degrees of freedom; therefore, transverse momentum or transverse 
mass distributions are used to extract temperatures. 
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Figure 4.9: Mass and collision system dependence of the inverse slope parameter T measured by the 
NA44 collaboration p3| . 

The inverse slope T is obtained from a fit to the transverse mass ttit spectrum using 

1 dN 



cx e 



m T /T 



, (4-4) 

Heavy ion collisions at the AGS and, even more at the SPS, exhibit an intricate scenario. Here, typical 
MB collision energies are in the range of 2-10 GeV in a baryon rich regime. Thus two problems have 
to be addressed: The influence of the modeling of meson baryon scattering as well as the creation of 



collective transverse flow (first measured at the SPS from the NA44 collaboration, see Fig. |4.9|) which 
lead to increasing apparent temperatures "T". To study these aspects, the UrQMD model was used in 
two different modes: 

1. A complete downward extrapolation of the high energy longitudinal color flux picture to \fs = 2 GeV 
(referred to as forward-backward peaked (f-b) variant). 

2. An upward extrapolation of the resonance behavior, i.e. an isotropic angular distribution of all 
outgoing particles in MB reactions (indicated as (iso)). 
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Figure 4.10: Mean p t at midrapidity (\y\ < 
0.5) function of particle mass in central 

Pb(160GeV)+Pb reactions. Figure taken from 



Figure 4.11: Mass and collision system dependence 
of the inverse slope parameter T in comparison 
to measured NA49 data [|195|| . Figure taken from 

pa. 



To gain insight into the creation of transverse flow, Fig. |4.10| discusses mean transverse momenta of 
particles for different mass bins at \y\ < 0.5 in central Pb+Pb collisions at the SPS. Circles show the 
p(160 AGeV)+p events, open squares indicate the (f-b) Pb+Pb events and the full black squares show the 
isotropic model (iso). In addition, we fit the resulting mass spectra with a simple fireball plus flow model. 
The full black line is the mean transverse momentum < pT(m,To, (3?) > for different masses m calculated 
from a thermal distribution with temperature To = 160 MeV without any additional flow (/3t = 0). The 
grey dashed line shows the < pT(m,To, (3t) > from an expanding thermal source (To = 170 MeV) with an 
additional transverse flow velocity of (3t = 0.34c. This yields in the saddle point approximation a rough 
value for the apparent temperature T ||196|| : 



T « T + m/?| 



(4.5) 



The proton proton reactions and the f-b scenario show no significant difference in mean p? to a non 
expanding thermal source. Only the isotropic model produces a transversally expanding hadronic source 
with significant additional transverse flow velocity of 0.34c. Note that the statistical error bars are only 
stated for the (f-b) Pb+Pb reaction, in the (iso) and p+p case they are of approximately the same order. 

This behavior of the (iso) prescription is also seen if we study the apparent temperatures (inverse 
slopes) at midrapidity (\y\ < 0.5) as a function of particle mass m as depicted in Fig. [4.11| . Proton proton 
collisions (crosses) show about the same freeze-out temperature for all particles from 130 MeV mass (pions) 
to 2 GeV mass. This is mostly the same in the forward backward scenario (high energy limit extrapolated 
downwards) in Pb+Pb. The inverse slopes increase only slightly with particle mass. Only the isotropic 
treatment of intermediate energy meson baryon collisions in line with the resonance picture yields an 
increasing apparent temperature with the particle mass as indicated by the preliminary NA49 data ||195 |. 
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The mass dependence of the apparent temperature is sensitive to the detailed modeling of the meson 
baryon rescattering process above the resonance region and below the high energy domain. Strong trans- 
verse expansion only builds up in the 'resonance model' in contrast to the high energy model which creates 
much less transverse flow. 

Thus, collective (radial) flow | |197| , 



as well as feeding from resonances, strongly influences the shape 
of particle spectra [|199| , [200|| . For light composite particles, the influence of collective flow is visible in a 
strong shoulder-arm shape of transverse momentum spectra. This can be seen in figure 4.12 . In order 
to account for flow effects, the spectra have to be fitted with a thermal distribution including collective 
flow. The temperature T and the transverse flow velocity (3 t are the fit-parameters. The shapes of the 
flow velocity profile and density profile at freeze out enter as additional degrees of freedom in the analysis. 
Usually a box shaped density profile and a linearly increasing transverse velocity profile are assumed. 




Figure 4.12: RQMD prediction of transverse momentum spectra at mid-rapidity for protons and light 
composite particles in central Si and Au collisions at the AGS. The bold solid histogram is a calculation 
with potentials, the thin histogram the respective calculation without potentials and the smooth solid line 
depicts a Boltzmann parameterization adjusted to the high momentum part of the spectra. The strong 
shoulder arm structure is due to collective flow. The figure has been taken from ||201 |. 
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In microscopic calculations, temperatures and flow velocities can also be extracted by subdividing the 
system into cells and analyzing the local transverse and longitudinal velocity distributions. Temperatures 
extracted via a global two parameter fit are more than a factor of two higher than temperatures obtained 
from such a microscopic analysis, at least at beam energies in the 100 MeV/nucleon to 10 GeV/nucleon 
regime | |128 . 

The reason for this discrepancy lies in the assumed shape of the freeze out density and velocity profiles. 
Analyses based on microscopic transport calculations show an extreme sensitivity of the extracted values of 
Pt and T on the shapes of the profiles ||123| , |201|| . Whereas a linearly increasing transverse freeze-out velocity 
profile seems a tolerable assumption the shape of the freeze-out density profile has a rather Gaussian shape 
(centered at r t — 0), far from the usually assumed box-like distribution. When using realistic density and 
velocity profiles one finds the high m t components of the particle spectra dominated by large collective 
effects (i.e. the high expansion velocity). This results in a lower value for the temperature T. A microscopic 
analysis of spectra of protons and light composite particles at AGS energies showed fit and T depending 
on the mass of the particle ||201 . 



Spectra of light secondary particles, such as pions, at SIS and low AGS energies are dominated by 
resonance decays: Here we will use pion spectra in Au+Au reactions at 1 GeV/nucleon as example for the 
importance of resonances. They show how the shape of the spectra and the origin of the spectra can be 
explained in a fully dynamical and non-thermal scenario: 

Figure |4.13| shows a comparison of inclusive 7r° spectra for Au+Au and Ca+Ca (minimum bias and 
y c .m. =t 0.16) between the IQMD model [ 111 , [202j| and data published by the TAPS collaboration |[203|| . 
The data have been rescaled by a factor of 0.6 in oder to obtain a minimum bias distribution from the 
measurement. The IQMD model overpredicts the yield for the heavy system Au+Au by approximately 
20% which is consistent with the findings of the KaoS collaboration for charged pions (a comparison has 
been published in ref. [ I 1 1[). However, the 7r° yield of the light system Ca+Ca is overpredicted by almost 
a factor of 2. For both systems IQMD calculations fail to reproduce the high transverse momentum tails. 

Newer transport models perform much better when compared to the data [|112| , [204|| . Figure |4.14| shows 
a comparison to the same data with the UrQMD model. Note that the slopes agree extremely well with 
the experiment and the yield for the system Ca+Ca deviates less than 20% from the data. However, the 
yield for the heavy Au+Au system is underpredicted for low transverse momenta. A reason for this may 
lie in the strength of the A1232 — N absorption cross section, which is calculated via the principle detailed 
balance from the NN — > iVAi232 cross section (see figure |3.4j) . The UrQMD fit to the A1232 is very good. 
The absorption channel, however, diverges for low energies (due to the detailed balance principle) and 
is therefore extremely sensitive to even the minutest changes in the cross section parameterization. A 
reparameterization of the exclusive A1232 production cross section (within the constraints given by the 
data) may therefore improve the pion yield in heavy collision systems. 

The importance of higher baryon resonances can be seen when comparing the IQMDQ pion spectrum at 
1 GeV/nucleon, figure |4.13| , with the respective spectrum generated by the UrQMD model in figure [4.14 . 
The enhancement at high p t in the UrQMD calculation - which is due to the incorporation of higher 
baryon resonances | |205| | - is clearly visible and UrQMD agrees with the data even at the highest momenta. 

Especially the measurement of those high energy pions is of great interest. They correlate directly to 
early freeze out times and heavy resonances [[206 1: 

Transverse momenta pt of pions in semi-central Au+Au collisions at 1 GeV/nucleon are shown in 
figure [4. 15| as a contour plot versus their freeze out time (the time of their final interaction in the heavy 
ion collision): High p t pions are produced almost exclusively in the early reaction phase with freeze out 

1 Note that the IQMD model includes only the A1232 and no other resonances are considered. 
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Figure 4.13: Comparison of inclusive n spectra 
for Au+Au and Ca+Ca (minimum bias) colli- 
sions between the IQMD model and data measured 
by the TAPS collaboration. A hard EoS including 
momentum dependence is used. Figure taken from 
TT1. 



Figure 4.14: Same analysis as in figure [4.13| , but 
with the UrQMD model: The influence of higher 
resonances on the high energy part of the spectrum 
is clearly visible. 



times smaller than 20 fm/c, well in the high density phase of the reaction. The scaling of the contour lines 
in figure |4.15| is linear, far higher transverse momenta are obtained than depicted by the contour lines. We 
can establish a correlation between high p t pions, early freeze out times and the hot and dense reaction 
phase. A correlation between high energy pions and early freeze out times has also been observed in BUU 
calculations of La+La collisions at 1.35 GeV/nucleon ||207 , 208]. However, the global, fixed event averaged 
time cut of 20 fm/c for La+La at 1.35 GeV/nucl. employed in ref. |[2Q7| , [208j| is a rough estimate of the 
nucleonic freeze out time (while in our approach the actual time (event by event) for pion freeze out is 
employed): The densities in the collision center are then already far lower than ground state density. 

The average freeze out density for all pions lies well below nuclear ground state density; this is well 
known and has severely diminished the usefulness of pions as probes for the hot and dense reaction phase 
58, |209| , |210|| . High p t pions, however, freeze out at far higher densities: Their average freeze out density 
lies between 1.2 and 1.5 p/po with some of them even freezing out at densities near the maximum density 
obtained in the collision. 

Figure |4.16| shows the transverse mass spectrum of tt + in Au+Au collisions at 1 GeV/nucleon within 
the IQMD model [206]. The spectrum has been decomposed into contributions of light (m^ < 1172 MeV), 
medium (1172 MeV < rrij\ < 1292 MeV) and heavy (ijia > 1292 MeV) A resonances. The correlation 
between high m t pions and heavy resonances is obvious. Note that since the IQMD model is limited to 
the A1232 resonance and the inelastic nucleon nucleon cross section is projected onto the excitation of this 
particular resonance this statement has to be interpreted in a more general manner: at high transverse 
masses heavy resonances such as the -Vj* 440 , iV-j^Q anc ^ -^1*535 should be substituted for the heavy A(1232) 
resonances of the IQMD model and contribute strongly to the transverse momentum or mass spectrum. 
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Figure 4.15: Transverse momentum versus freeze 
out time for pions in semi-central Au+Au colli- 
sions at 1 GeV/nucl. in the IQMD model ||206 . 
High p t pions freeze out predominantly in the hot 
and dense reaction phase. The scaling of the con- 
tour lines is linear; far higher transverse momenta 
occur than is depicted by the contour lines. 



Figure 4.16: Transverse mass m t distribution of 
freeze out 7r + with cuts on the mass of the emit- 
ting A resonance. The calculation are done for 
Au+Au collisions at 1 GeV/nucleon within the 
IQMD model. The connection between high mt 
and large resonance mass is obvious: The produc- 
tion of freeze out pions with high m t is dominated 
by the decay of heavy A resonances. 



4.2.3 Hadrochemistry: multiplicities and ratios 



Hadron abundances and ratios have been suggested as possible signatures for exotic states and phase 
transitions in dense nuclear matter. In addition they have been applied to study the degree of chemical 
equilibration in a relativistic heavy-ion reaction. Bulk properties like temperatures, entropies and chemical 
potentials of highly excited hadronic matter have been extracted assuming thermal and chemical equilib- 
rium 0, 0, g |2Tl], |212], |2T|, §, gl|, [m| [2l|, [217], [213, [21|, [2Sg, g. Recent SPS data on hadron yields 
and ratios have been fit either in the framework of a hadronizing QGP droplet [ |211| , 222, 223, [224|| or of 
a hadron gas in thermal and chemical equilibrium [ |221| | (including elementary proton-proton interactions 
144j| ). It has been shown that the thermodynamic parameters T and fis imply that these systems have 
been either very close to or even above the critical T, line for QGP formation ||225| , 211 1 . However, a 
UrQMD calculation ||226|| agrees with many of the data (n/p, d/p, p/p, A/A, H/A etc.) without assuming 
thermal and chemical equilibrium. Large discrepancies to the data (> 50 %) are only found for the K s j A 
and ratios (see figure f4.17| ). The microscopic transport model calculation shows that secondary in- 
teractions significantly change the primordial hadronic cocktail of the system. The ratios exhibit a strong 
dependence on rapidity [ |226| . 

The assumption of global thermal and chemical equilibrium - as being employed by the static thermal 
models - is not justified: Both, the discovery of directed collective flow of baryons and anti-flow of mesons 
in Au+Au reactions at 10.6 GeV/nucleon ||239|| and Pb+Pb reactions at 160 GeV/nucleon energies ||240| 



(see also section |4.3|) as well as transport model analyses presented in section [4.2.7| , which show distinctly 
different freeze-out times and radii for different hadron species, indicate that the yields and ratios result 
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Figure 4.17: comparison between the UrQMD 
model and data for the system S+Au(W,Pb) at 
200 GeV/nucleon. Also shown is a fit by a sim- 
ple hadronization model ||211| . Both models agree 
quite well with the data compilation ||221| . The 
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Figure 4.18: IQMD prediction for the 7r~ to 7r + ra- 
tio versus transverse momentum p t in Au+Au col- 
lisions (minimum bias) at 1 GeV/nucl. with a hard 
EoS and momentum dependence ||1 1 1|| . The solid 
line shows the full calculation including Coulomb 
forces. The dashed-dotted line shows a calculation 
without Coulomb forces acting upon the pions. 



from a complex non-equilibrium time evolution of the hadronic system. A thermal model fit may therefore 
not seem meaningful. 

On the other hand, particle ratios also contain valuable information about the time-evolution and the 
dynamics of the collision system: The inclusion of Coulomb forces and energy dependent cross sections is 
also important in the domain of particle production. Figure |4.18| demonstrates this by showing the IQMD 
prediction 1 11 1 1 ] for the 7r~ to 7r + ratio versus transverse momentum p t for Au+Au at 1 GeV/nucleon. 
The solid line shows a full calculation including Coulomb forces. For high p t the tt~/tt + ratio decreases 
towards 1, whereas for low pt it increases to 2.5 - considerably higher than the value of 1.9 predicted by 
the A-isobar model. The dashed line shows a calculation without Coulomb forces. This ratio remains 
constant at ~ 1.9. The (small) remaining variations might be due to the different energy dependence of 
the 7r + — p and tt~ — p inelastic cross sections. The IQMD prediction is compatible to recent measurements 



by the KaoS collaboration ||241|| . A detailed comparison of the data to the UrQMD model has yet to be 
performed. 

Note that the pt-integrated values agree very well with the A-isobar model and thus hide any further 
information possibly contained in the p t dependent ratio. The A-isobar model | |242| ! was originally devel- 
oped to describe p-wave pion-nucleon scattering and is based on considering the A 12 32 as an independent 
baryon species. Here it is applied in a statistical limit, i.e. it is the basis for the assumption that all pions 
are produced via the A1232 and the resulting multiplicity ratios only depend on the number of incident 
protons and neutrons and can be calculated via combinatorics applying the respective Clebsch-Gordan 
coefficients for the A 12 32 excitation and the subsequent decay into pions. 
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Figure 4.19: By using projectile and target with different isospin content, the mechanism of stopping, i.e. 
the amount of penetration between projectile and target nucleons, can be measured by plotting the neutron 
to proton ratio (or related quantities) versus the rapidity. Three possible scenarios are schematically 
depicted in the figure together with their respective signature. 



Particle ratios like the proton to neutron or the experimentally accessible 3 He to triton ratio may 
also be used to gain experimental information on the mechanism of stopping (i. e. the amount of inter- 
penetration between projectile and target nucleons). So far the analysis of the stopping mechanism has 
only been accessible to transport model calculations. We propose an experiment with which regions in 
rapidity dominated by projectile or target nucleons can be identified (even at central collisions!) and with 
which therefore the question of the stopping mechanism can be investigated. If projectile and target have 
a large difference in their total isospin (i.e. they differ strongly in their neutron to proton ratio) then the 
comparison of the neutron to proton ratio vs. rapidity before the reaction with the respective ratio after 
the reaction yields direct information on the degree of penetration between projectile and target nucleons. 
Figure [4.19| depicts three extreme cases: (a) full stopping, (b) complete mixing of projectile and target 
nucleons and (c) transparency. Note that the configuration in momentum space is opposite to that in 
coordinate space. 

We have used the QMD/IQMD model in order to test the feasibility of the suggested observable. As 
a test system we selected 50 Cr+ 48 Ca which has a large difference in the neutron to proton ratio at almost 
identical mass. The calculations were performed with impact parameters b< 2 fm. 

Figure [4.20| a shows the neutron to proton asymmetry vs. rapidity for 50 Cr+ 48 Ca and the symmetric 
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Figure 4.20: Neutron to proton asymmetry as obtained in IQMD model for 50 Cr+ 48 Ca and 50 Cr+ 50 Cr at 
1 GeV/nucl. (a), for 50 Cr+ 48 Ca at 150 MeV/nucl. vs. 1.5 GeV/nucl. (b), for the same system without 
vs. with cluster subtraction (c) and with an enhanced (factor 5) cross section (d). 



system Cr+ Cr, both at 1 GeV/nucleon. The horizontal lines mark the values expected for case (2) of 
fig. |4.19| for the isospin asymmetric (dotted) and symmetric system (dash-dotted). The calculations clearly 
indicate that projectile and target move through each other although the nucleon rapidity distribution 
(which usually characterizes stopping) shows a peak at cm. rapidity. Figure |4.20| b shows the energy 
dependence of the asymmetry. At 1.5 GeV/nucl. the transparency is far more pronounced than at 150 
MeV/nucl. (which is within error bars almost compatible with case (2) of Fig. f4.19p . However, Figs, (a) 
and (b) only show nucleons without taking clusterization into account. Fig. [4.20| c compares an analysis 
based on (all) nucleons with the respective analysis in which clusters were subtracted (250 MeV/nucl. 
incident beam energy). The qualitative trends remain unchanged. In order to investigate the sensitivity 
of the asymmetry to the stopping power we have performed a calculation with a 5 times higher nucleon- 
nucleon cross section. Fig. |4.20| d shows the respective calculation at 250 MeV/nucl. and compares it with 
the default calculation. A sign reversal in the slope of the asymmetry is visible. 

Experimentally it may be easier to study the 3 He to triton ratio instead of the neutron to proton ratio. 
At beam energies above 1.5 GeV/nucl. the isospin content of the system is transferred to pionic degrees 
of freedom and the neutron to proton ratio looses its effectiveness. 
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Figure 4.21: Equation of State e = e(T) obtained 
in UrQMD box calculations (squares) and in ideal 
hadron gas (solid line) at normal nuclear matter 
density and zero strangeness. 
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Figure 4.22: The delta to nucleon ratio vs. tem- 
perature as obtained in UrQMD infinite matter 
calculation. The gray dashed area shows the delta 
to nucleon ratio calculated from the law of mass 
action in a Boltzmann approximation, taking fluc- 
tuations in the delta mass into account. 



4.2.4 Equilibrium properties of hadronic matter 

The possibility of thermal and chemical equilibration of hadronic matter is studied using UrQMD box 
calculations with periodic boundary conditions for normal baryonic density and zero strangeness at different 
energy densities from e = 0.2 to 2.0 GeV/fm 3 . Starting from a state of pure nucleons with random 
momentum distribution, one observes that particle multiplicities equilibrate after some time. At low energy 
density e, equilibration times for strange particles are much larger than for non-strange particles, but this 
difference decreases with growing e. In equilibrium, the energy spectra of the different hadronic species are 
nicely fitted by a Boltzmann distribution exp (-E/T) with the same temperature T. Figure [4.21 shows 



the variation of the temperature T versus energy density e at ps = Po and ps = 0, obtained in UrQMD. 
The solid line shows the results of an ideal gas model including the same hadrons as UrQMD (tables 
and |3.3|). The UrQMD Equation of State shows a Hagedorn-like shape with a limiting temperature 



T rs 130 MeV, while the hadron gas EoS shows a continuous rise of the temperature with e. Both models 
are however in good agreement at low energy densities e < 0.3-0.35 GeV/fm 3 . Because the UrQMD model 
uses a stochastic collision term, and no hard core repulsion is considered for the different particles (cascade 
mode), excluded volume correction is not included in the ideal gas formulation. This correction, however, 



could be significant ||243|| and could be treated in the Chappman-Enskog approximation. 

In particular, the delta to nucleon ratio is consistent with the theoretical expectation of a hadron gas. 
This can be seen in figure [4.22| which displays the delta to nucleon ratio obtained at various energy densities 
ranging from 0.16 to 0.3 GeV/fm 3 . The gray shaded area shows the delta to nucleon ratio calculated from 
the law of mass action in a Boltzmann approximation, taking fluctuations in the delta mass into account. 
At higher energy densities, despite the smaller temperature in UrQMD compared to ideal gas, the number 
of pions and other mesons is larger than in the Ideal Gas. The difference between the two models resides 
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in the string degrees of freedom included in UrQMD ||244| . At high energy densities, they keep a large 
fraction of the system energy. Moreover, the inclusion of strings in UrQMD leads to a violation of the 
detailed balance at high energies, assumed in statistical model. Box calculations in a simplified version 
where no strings and no many-body decays were included, show a very good agreement with the ideal gas. 
They nicely reproduce both temperature slopes and different hadronic multiplicities ||244 . 




Figure 4.23: Number of particles as a function of cms time, obtained in the central zone (5x5x5 fm 3 ) of 
central (b =0) Au+Au collisions at 10.7 AGeV. Solid lines correspond to all hadrons in the zone, dashed 
lines - to only formed hadrons. Points are the predictions of ideal gas calculated at the same energy, 
baryonic and strangeness densities as in UrQMD full sample in the zone. 

The possibility of local equilibration in Au+Au collisions at AGS energies (10.7 AGeV) is studied by 
analyzing the time evolution of hadronic matter in the central reaction zone [ 245|| . It is found that between 



t cm = 10 and 20 fm/c, there is an isotropic velocity distribution of hadrons in the central region. In this 
time interval, there are practically no strings left in this zone. Figure 14.23] shows the time evolution of 
different particle multiplicities in the central region. The solid lines correspond to the total number of 
particles for the different species, whereas the dashed lines show the multiplicities of only already formed 
particles (cf. section |3.3.5|) . The dots depict the results of the statistical model, obtained at the same e, 



pB and ps as in the central zone of UrQMD calculation. 

Baryonic density pb drops from « 2p at t cm =10 fm/c to 0.5 p at tcm =15 fm/c. Energy spectra of the 
different species are nicely reproduced by Boltzmann distributions with the temperature predicted by the 
statistical model (fig. |4.24j ). The statistical model reproduces also the multiplicities of baryons and mesons 
obtained in UrQMD (fig. [4 . 2 5| ) . One can see from fig. |4.21| that for the energy densities reached in the 



central zone of Au+Au collisions at AGS at these times (e.g. e = 293 MeV/fm , Pb = Po and ps ~ for 
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Figure 4.24: Energy spectra of N (•), A (□), tt (A) and K + (y) in central zone of Au+Au collisions 
at 10.7 AGeV at t cm =13 fm/c fitted by Boltzmann distributions with parameters T=128 MeV, /ie=534 
MeV, ijls — 112 MeV, predicted by the statistical model (lines) and obtained in UrQMD box calculations 
(open symbols). 



t cm =13 fm/c) the equation of state of the UrQMD model agrees quite well with that of the ideal resonance 
gas. The comparison with the equilibrated state obtained in box calculations at the same e, ps and ps 
(fig. 4.24 44.251) also indicates that in the central zone of Au+Au collisions, one probably reaches local 
chemical and thermal equilibrium. However, the existence of collective flow [233, 246 1 as well the analysis 
of the freeze-out properties of nucleons, hyperons and various meson species in section |4.2.7| indicate that 
the final state does not originate from an equilibrium configuration. One therefore has to be very careful 
when trying to draw conclusions from the existence of a local zone of hadrochemical equilibrium to the 
final state of the reaction. 
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Figure 4.25: Yields of hadrons obtained in central zone of Au+Au collisions at t cm =13 fm/c (*), in 
statistical model (Q) and in UrQMD box calculations (□). 



77 



4.2.5 Resonance matter 



The possibility of producing A-matter (or in more general terms: resonance matter) and density isomers 
has been already discussed early on ]14|, g. Recently this topic has received renewed attention [|247| , |248 



249| , |250| , |251| , [252| , [253|| : At beam energies above a few hundred MeV/nucleon, the nucleons can be excited 
into A-resonances. If the density of these resonances is as high as the nuclear matter ground state density, 
then a new state of matter, A-matter, has been created. One of the potential signals for the presence of 
A-matter is the creation of pions as decay products of the A- resonance. 

How can A-matter be produced? In order to address this question we first limit ourselves to Au+Au 
collisions at 1 GeV/nucleon. In that case the hadrochemistry is well described by a simple closed system 
of nucleons, pions and the A 12 32 resonance: 

Figure 4.26| shows the pion - nucleon cycle in the IQMD model, which is limited to nucleons, pions 
and the A1232 resonance. The scheme describes (for impact parameters b < 5 fm and averaged over 60 
fm/c possible processes linked to the creation of A-matter. The probabilities in the boxes always refer to 
the vertices they are directly connected with. A-resonances are initially produced via inelastic nucleon 
nucleon scattering. The produced resonances can either be reabsorbed via inelastic scattering or decay 
by emitting a pion. The pion can then either freeze out or interact with a nucleon to form a A again. 
In case the A has been absorbed the corresponding high energetic nucleon might have a second chance 
of becoming a A by inelastic scattering. It could also transfer energy via elastic scattering onto another 
nucleon which then could scatter inelastically and form a new A. A nucleon interacts in the average about 
three times before it freezes out. This value fluctuates considerably, depending on whether the nucleon 
was in the participant zone (geometrical overlap of the colliding heavy ions) or in the spectator zone of 
the collision. 

Unfortunately, the probability for a nucleon to undergo inelastic scattering and to form a A during 
the heavy ion collision at 1 GeV/nucleon is as low as 10%. The main process for sustaining A-matter is 
the A — > Nir — > A loop, which, however, first has to be fueled by the NN — > AiV process. The average 
pion passes approximately three times through this loop (it has been created by the decay of a hard A) . 
However, 30% pass more than 6 times through the loop. For nucleons the probability of forming a soft 
A i.e. via ttN — > A is almost twice as high (A-matter pump) than the probability of forming a hard A 
via NN — > NA. This picture is however much more complicated when going towards higher energies and 
higher resonances are considered. We show for example in fig. [4.271 , the NJ 535 -pump (time integrated) as 
obtained in UrQMD calculations for Ni + Ni collisions at 2 GeV/nucleon. One sees from the figure that 
only half of //-mesons produced at this energy are coming from N^ 535 decays. The other half comes mainly 
from meson-baryon collisions like Nn — > Ar] or from other resonances like the NJ 710 . 

Figure [4.28| a shows the IQMD time evolution of the total baryon, nucleon and A densities in units of 
p/po for Au+Au collisions at 1 GeV/nucleon. The densities are calculated in a sphere of 2 fm radius around 
the collision center. Between 5 and 20 fm/c more than 20 A-resonances can be found in the whole system: 
This time interval coincides with the hot and dense reaction phase. At 10 fm/c up to 55 resonances are 
present in the total reaction volume (keep in mind this is not in the 2 fm test sphere). A A multiplicity of 
> 40 can be sustained for an interval of 10 fm/c, 6 times longer than the lifetime of a free A-resonance. 
However, this is not pure A-matter. in the small test volume shown in figure |4.28| a the resonance density 
is 0.5 po and the nucleon density is 2.2 po : the A-contribution is 20% in the test volume which contains, 
as a matter of fact, only 2.5 resonances. The total multiplicity of A resonances is just about 10% of the 
total nucleon multiplicity. 

However, it is obvious that the other As can be distributed all over the reaction volume. Figure 
4.28| b shows the A density distribution as experienced by the A's in the system at 5, 10 and 20 fm/c. 
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The densities were calculated by summing over all contributing Gaussians of all As in the system at the 
locations of the respective As. We would like to point out that the mean A-density experienced by the As 
is about 0.25 p . Less than 1% of the As experience A densities around 0.5 p . As show collective flow in 
the reaction plane. Its measurable signature (the pion p x (y) distribution in central collisions) is discussed 
in section [4.3. 1| . 

At AGS energies, however, RQMD calculations predict an excited state of baryonic matter, dominated 
by the A1232 resonance. Analysis show a long apparent lifetime (> 10 fm/c) and a rather large volume 
(several 100 fm 3 ) for this A matter state in central Au+Au collisions at the AGS [25jJ (see figure |4.29|) . 
At higher energies (SPS), UrQMD calculations shown in fig. |4.30| , indicate that in the early stages of 
the reaction, about 75% of baryons are excited resonances (among which 55% are nucleon and delta 
resonances and 20% are hyperons). This fraction is even higher at mid-rapidity (about 85%). The excited 
hyperons are created in the first 3 fm/c and remain almost constant for the whole propagation. This 
fast strangeness equilibration is different from the time scales for strangeness equilibration predicted by 
UrQMD box calculations at an energy density equivalent to that obtained in the central region of a Au+Au 
collision at 160 GeV/nucleon |[255 |. Because of the large fraction of hyperons during the whole propagation, 



there is a finite probability for the formation of hyperon clusters (hyperon-matter). 
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Figure 4.26: Pion - nucleoli cycle in the IQMD model for Au+Au collisions at 1 GeV/nucleon. The 
scheme describes (for b < 5 fm and time-averaged) all possible processes (in the IQMD model) linked 
to the creation of A-matter. The probabilities in the boxes always refer to the vertices they are directly 
connected with. The main process for sustaining A-matter is the A — > Nir — > A loop, which, however, first 
has to be fueled by the NN — > AN process. 
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Figure 4.27: N^ 535 -pump as obtained in UrQMD calculations for Ni + Ni collisions at 2 GeV/nucleon. 
The scheme describes (for b < 5 fm and time-averaged) all possible processes linked to the creation of 
NJ 535 . The probabilities in the boxes always refer to the vertices they are directly connected with. 
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Figure 4.28: IQMD time evolution of the total baryon, nucleon and A-resonance density in units of p/po 
(a) and A density distribution the respective As experience for 5, 10 and 20 fm/c (b) for Au+Au collisions 
at 1 GeV/nucleon. The densities in the upper frame (a) are calculated in a sphere of 2 fm radius around 
the collision center. The densities in the lower frame (b) were calculated by summing over all contributing 
Gaussians of all A's in the system at the locations of the respective A's. 
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Figure 4.29: Time evolution of particle multiplicities (scaled with the number of incident nucleons) for 
central Au+Au collisions at 1 GeV/nucleon (SIS) and at 10.6 GeV/nucleon (AGS). At SIS energies, only 
about 10% of the nucleons are excited to resonances whereas at AGS energies the degree of excitation 
exceeds 50%. For a time-span of up to 10 fm/c the baryons are in a state of A— matter. The figure has 
been taken from [1254 
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Figure 4.30: Fraction of excited resonances vs. time for central Au+Au collisions at SIS (solid line), AGS 
(dotted line) and SPS (dashed-dotted line) energies. 
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4.2.6 Antimatter and strange matter 

As mentioned above, the dominant reaction mechanism in the early stage of a reaction is the excitation 



of collision partners to resonances or strings |[37| . Then secondary interactions, i. e. the annihilation of 
produced mesons on baryons, lead also to the formation of s channel resonances or strings, which may 
explain the strangeness enrichment ||184j , |185| , |160|| and (for masses larger than 3rriN) allow for NN creation 
|190|| . It is clear, however, that any subsequent interaction of the newly produced particle can change the 
final yields and spectra, which has to be taken into account. E.g. the escape probability for p's from 
the exploding nuclear matter enters into microscopic models usually via the free NN annihilation cross 
section. 

The two counter-acting effects (production vs. absorption) may be measured by the directed "antiflow" 
of antimatter. Varying the formation time r, i. e. applying r = 6 fm/c instead of the default string 
fragmentation formation time (r( anti _) baryons « 1.5 fm/c) reduces the absorption from « 90 % to ~ 50 % 
and the flow (to one half) for heavy systems. 

The observable asymmetry for antimatter can be quantified by dividing the yields of the upper and lower 

hemispheres for target and projectile rapidity separately. The ratio R 

then gives the probability for a single particle to follow the collective (directed) flow. The following table 
shows the sensitivity of the azimuthal asymmetry R and the antiflow-parameter (the mean p x , integrated 
over the whole target or projectile rapidity) on the impact parameter for Au+Au reactions at 10.7 GeV/u 
within the RQMD model. The anti-flow of antiprotons appears to be strongest for semi-central collisions 
of b ~ 5 — 7 fm, while for protons the maximum p x is at considerably smaller 6-values. The latter is 
due to the pressure (i. e. the EoS), the former one due to absorption and geometry. A similar absorption 
mechanism is at work for all hadrons with valence anti-quarks, e.g. 7r's and anti-kaons, resulting in mesonic 
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Now, let us turn to even more complex and sensitive observables, e.g. clusters of antinucleons. The 
antideuteron formation rate (d/p 2 -ratio) in Si+Au collisions at the AGS is sensitive to the shape and 
size of the antinucleon source [[257| , |258|| . It is different from the nucleon source. Annihilation strongly 
distorts a homogeneous coordinate space distribution of antibaryons. Asymmetric phase space distributions 
of anticlusters result for finite impact parameter b, visible as directed transverse anti-flow (fig. |4.32|) . 
Antimatter cluster formation is strongly suppressed. 

In analogy to the deuteron case (see figure |4.12|) , the (anti-)deuteron formation is calculated by project- 
ing the (anti-)nucleon pair phase space distribution on the (anti-)deuteron wave function via the Wigner- 
function method. The (anti-)deuteron Wigner-density (/>j ) is again given by the Wigner-transformed 



Hulthen wave function. Thus the recently measured Si(14.6AGeV)Au antideuteron data at the AGS ||259 
can be reproduced! 

Fig. |4.33| shows the calculated coordinate space distribution of mid-rapidity antideuterons in a cut 
perpendicular to the reaction plane for the Au(10.7 AGeV)+Au, b = 5 fm reaction. Anti-matter is only 
emitted from the surface of the fireball - i.e. the region of hot and dense matter - in clear contrast to the 
source of baryons, which spreads over the whole reaction volume. This can be understood as a consequence 
of antibaryon absorption, which produces deep cuts in the momentum- and even more so in the coordinate 
space distribution of the antibaryons at freeze-out ||190|, |258||. p's in-plane are shadowed more effectively 



due to the presence of spectator matter. This leads to an antimatter cluster squeeze-out. 
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Figure 4.31: Mean directed transverse momentum 
(in-plane) of various hadrons in Au+Au collisions 
at 10 GeV/u function of rapidity (RQMD cal- 
culation). Fig. taken from ||190|| . 
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Figure 4.32: Mean directed transverse momentum 
(in-plane) of anti-clusters in Au+Au collisions at 
10 GeV/u function of rapidity (RQMD calcu- 
lation). Fig. taken from |[256|| . 
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Figure 4.33: Coordinate space distribution of mid- 
rapidity anti-deuterons as obtained in RQMD cal- 
culations for the reaction Au(10.7 AGeV)+Au at 
freeze-out in a cut perpendicular to the reaction 
plane. (Fig. taken from |260|| ) 



Figure 4.34: Impact parameter dependence of the 
anti-deuteron to the anti-proton ratio squared (full 
circles) and of the deuteron to the proton ra- 
tio squared (open squares) both at mid-rapidity. 
Open circles: the result of a simple momentum 
coalescence model. (Fig. taken from ||26U[|) 



Let us explore the reaction volume dependence of the (anti-)deuteron formation by varying b as shown 
in fig. [4.34| (d/p 2 -ratio, full circles). The deuteron formation rate (d/p 2 -ratio, open squares) is proportional 
to V~ l . Those rates measure the average phase-space distance of N and iV's, respectively. 

In central collisions up to 95% of the initially produced antinucleons are reabsorbed. 

The d formation rate is predicted to be roughly two orders of magnitude lower than the formation 
rate of deuterons. This difference vanishes when going to high impact parameters or small systems (like 
Si+Al). The assumption of the independent production of both antinucleons becomes less justified in very 
small systems. Calculating the antideuteron formation rate within a momentum coalescence model with a 
Ap-parameter of 120 MeV (<i/p 2 -ratio, open circles) one finds only small sensitivity to the chosen impact 
parameter. This proves the pure coordinate space nature of the effect. 

Fig. |4.35| shows the azimuthal distribution of mid-rapidity p, d, i in peripheral Au+Au collisions. The 
momentum distribution of p in the p x — p^-plane is nearly isotropic. It reflects the geometry of the almond 
shaped reaction zone. In contrast the d distribution is shifted towards = 90° in line with the described 
coordinate space distribution of the d. This looks like a squeeze-out effect. It is even more pronounced for 
t. 

The predicted reaction volume dependence and the squeeze-out of anti-fragments reflect the spatial 
distribution of antibaryons. This suppression of antideuterons may mask a possible antimatter cluster 



enhancement in a quark-gluon plasma event [261 



An enhancement of the production of anti-nucleons due to the presence of relativistic mean fields 
was predicted in [|177] , |262|| . It has been calculated in the framework of RBUU [|179| , |180|| , QMD ||181|| and 
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Figure 4.35: Correlations of antiprotons lead to a squeeze-out for d and t in peripheral Au(10.7AGeV)+Au 
reactions due to a higher anticluster formation probability at the ends of the almond shaped fireball (RQMD 
calculation). Fig. taken from ||260 . 



RVUU ||182|| that enhanced antiproton production due to reduced effective masses may explain subthreshold 
antiproton measurements at SIS energies. On the other hand, these transport model calculations show 
that most of the additional anti-nucleons get reabsorbed. The absorption rate is much higher than in 
a cascade without mean fields [183], because this additional production mechanism is most effective in 



regions of high baryon density and because of the additional attractive interaction between nucleons and 
anti-nucleons. 

Enhanced production of anti-nucleons due to mean fields or a quark-gluon plasma phase is conceivable 
at higher energies, i.e. at AGS and SPS, too. However, UrQMD calculations of a Si-Al collision at 
Eiaf, = 13.7 GeV/nucleon ||263|| without mean fields indicate that even at these energies, an enhanced 
production of anti-nucleons will be difficult to observe. Figure |4.36| shows the time evolution of 16 anti- 
nucleons (8 • p, 8 • n) that are added by hand after 4 fm/c in the central region of the collision. One can 
see that nearly no anti-nucleon survives the further evolution of the system (less than 4%). In Au+Au 
collisions, where higher baryon densities are predicted, all additional anti-nucleons are reabsorbed after 
1-2 fm/c. 

We have learned that iV's suffer strong final-state interactions. These interactions have in principle two 
components which can be related to the N self-energy in matter: collisions and annihilation on baryons [[26 -| 
(imaginary part, semi-classically given by 2 Im V = avp) and a piece in the real part (KeV = £nnP in the 
impulse approximation). In the semi-classical limit the real part of the self-energy can be approximated 
by potential-type interaction PSR |266|| or a mean field. 

Here we will focus on the effect of the real part. The motivation is that the long-range force of baryons 
acting on a p is expected to be stronger than for protons since the Lorentz-scalar and the Lorentz-vector 
parts of a meson exchange potential now have the same sign. The influence of baryonic mean-fields on 
baryons and mesons is well established. Therefore there should also be some influence on p's. 

The strong imaginary part of the optical AW-potential has been taken into account by using the 



geometric annihilation cross section. The real part may also be substantial |[267|| . However, it would be fatal 
to "improve" the description by just adding a mean field interaction, since a part of the real potential has 
already been taken into account by the (parameterized) free elastic and inelastic cross sections. Fig. |4.37 
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Figure 4.36: Time evolution of 16 anti-nucleons (8-p, 8-n) in a Si-Al collision at E iab = 13.7 GeV/nucleon 
calculated in the UrQMD-model (see text). 
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shows this effect for our model potential. Due to the strong attraction for the BB case, a reduced 
geometric annihilation cross section suffices to account for the measured free annihilation probability in 
binary pp reactions. To avoid the inevitable double counting in fully dynamical calculations this must 
not be neglected. In the case of a homogeneous nuclear medium the attractive potentials can partly 
cancel, leaving only the geometrical cross section effective for N annihilation. The importance of the real 
potential for the annihilation cross section results in an additional uncertainty of the Y and A annihilation 
probabilities, since the real part of the self-energy is even less known than for p's. This should be taken 
into account for the analyses of (preliminary) reports of surprisingly high A/p ratios at AGS and SPS 
experiments H, |6|, [27(1. 

As a first step, we restrict our study on the residual potential interaction with its finite range to the 
dynamics after the hadronic freeze-out. The success of Dirac equation optical model calculations for pA 
scattering |[271|| leads us to using these relativistic potentials with Yukawa functions as interaction form 
factors — applying G-parity transformation — for the p case: The mass parameters are [iy = 3.952 fm~ 
and us = 2.787 fm~ , the coupling constants are gy = 2674.5 MeV fm and gs = 2158. 2MeV fm. In line 



with [271 1 Gaussians are used as baryon profiles with a mean square radius of 0.8 fm. The central part of 



an effective Schrodinger equivalent potential (SEP) is constructed from the above potentials: 

U SE = ±(2EU v + 2mU s -U$ + U$) , (4.6) 

where E is the total energy of the incident particle. 

For small NN distances the real potential should not show an effect, since the huge imaginary part 
absorbs the particular p anyhow. We have chosen 1.5 fm as cutoff distance corresponding to an averaged 
pp-annihilation cross section of ~ 70 mb. 

The Schrodinger equivalent potential with the above parameters results in a mean p-potential in nuclear 
matter of about -250 MeV (p^ = GeV/c), increasing with energy towards -170 MeV at p^ = 1 GeV/c. 
The actual (averaged) SE-potential of the p's at freeze-out is about -70 MeV. 

Due to the strong momentum dependence of the pp-annihilation cross section, p's with low transverse 



momentum are suppressed in RQMD 1.07. Figure 4.38 shows the invariant multiplicity of p's with p t < 



200 MeV for central collisions of Au+Au at 10.7 AGeV. Calculations with and without potentials are 



compared to preliminary data of the E878 collaboration [272]. Besides the proposed model interaction we 
calculated the effect for the same potentials, arbitrarily reduced by 50 %. Still the dip at mid-rapidity 
vanishes although the change is less pronounced. 

Note that the p r integrated spectrum is not considerably affected by the potential interaction, but 
the final phase space distribution of the p's at \ow-p t deviates substantially from the standard RQMD 
calculation: A clearly non-thermal spectrum with a dip at mid-rapidity for p t = gets distorted by 
the potential interaction during the last stage of the collision. Up to now, the possible influence on 
the formation of anti-clusters, which have just proved to be a delicate probe of the final phase-space 
distribution, is unclear. 

Mean-fields are not only expected to influence the phase space distributions of hadrons, but also might 
have a strong impact on the production probabilities (see above). Relativistic meson- field models, which, 
at high temperature Tq « 180 MeV qualitatively simulate chiral behavior of the nuclear matter, exhibit 
a transition into a phase of massless baryons |[275|| . Every (anti-)baryon species (hyperons included ||25|| ) 
shows approximately the number density of normal nuclear matter (po ~ 0.16 fm -3 ) near Tq. Thus, 
the fraction of (anti-)strange hyperons increases by 1-2 orders of magnitude at Tq. Several hundred 



(anti-)baryons, the majority being (anti-)hyperons, may then fill the hot mid-rapidity region. Fig. [4.39 



shows this transition for (anti-)nucleons and (anti-) hyperons at /x g =100 MeV. Above Tq all effective masses 
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Figure 4.37: pp annihilation cross section as a func- 
tion of the lab-momentum. Parameterization of 
the free measured cross section (full line) and the 
corrected value, if the described potential inter- 
action is added (dashed line). (Fig. taken from 

m 



Figure 4.38: Invariant rapidity-distribution of the 
p's with p t < 200 MeV for Au+Au {b < 4 fm) 
at 10.7 GeV/u. Shown is the RQMD calculation 
(dashed line), with the additional optical potential 
(full line) and the 50 % reduced potential (points). 
Preliminary data (full circles) from [272]. (Fig. 
taken from g7J] ) 
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are small and the relative yields in the medium are dictated by the isospin degeneracy, thus favoring (anti- 
)hyperonic matter. This scenario could enhance the production probability of (anti-)hyperons. Even the 
formation of meta-stable exotic multi-strange objects (MEMO's |]25|| ) is conceivable via this mechanism. 
However, the high temperature will suppress the formation of clusters of mass A by e~ A ^ m ~ m ^ T . 

It was speculated that strange matter might exist either as meta-stable exotic multi-strange objects 
(MEMO's p5| ) or in form of strangelets. MEMO's are baryonic clusters which contain hyperons, thus 
constituting lumps of hypermatter. Strangelets, on the other hand, can be thought of as (meta-)stable 
multi-strange quark clusters p8|, [276|| . 

The possible creation — in heavy ion collisions — of long-lived remnants of the quark-gluon-plasma, 
cooled and charged up with strangeness by the emission of pions and kaons, was proposed in p3| , |277 , 
278| , |279|| . Strangelets can serve as signatures for the creation of a quark gluon plasma. The detection of 
strangelets would verify exciting theoretical ideas with consequences for our knowledge of the evolution 
of the early universe |30], |31fl ) the dynamics of supernova explosions and the underlying theory of strong 
interactions. Currently, both at the BNL-AGS and at the CERN-SPS, experiments are carried out to 



search for MEMO's and strangelets, e. g. by the E864, E878 and the NA52 collaborations [280, 269 



The production of strangelets and anti-strangelets might be possible via the strangelet distillation 
mechanism |[278| , [279|| if low values of the bag-constant B 1 ^ < 180MeV are assumed. However hadronic 
ratios at SPS energies can be reproduced with the same hadronization model only under the assumption 
of a bag constant of B 1/A > 230MeV |M] (see also |223|, |224], M ). This would exclude the formation of 
meta-stable strangelets within this model. 

Be also reminded that the question, whether strangelets or MEMO's can exist as bound states at all, is 
very speculative and thus still a controversial point, on which we did not focus here. Special (meta-)stable 
candidates for experimental searches are the quark- alpha |32| with Ab = 6 and the H-Dibaryon with 
A B 
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4.2.7 Particle freeze-out 

The previous sections have dealt with the description of the final state of heavy ion reactions. Different 
analyses of data at AGS and CERN/SPS energies have shown that the experimental results may be 
described as a transversally and longitudinally expanding chemically and thermally equilibrated hadron- 
gas ||217| , |220| , |221| , |284|| . The assumption of a chemically and thermally equilibrated source, however, 
implies a uniform freeze-out for all particles; i.e. all particle species have the same freeze-out times, radii 
and densities. Such a scenario is easily tested in the framework of a microscopic calculation. This section 
therefore deals with the investigation of freeze out properties of different particle species in the framework 
of microscopic transport models and investigates model-independent methods to gain access to freeze-out 
information contained in the final state of the heavy ion reaction. 

Figure [4.40| shows the time-evolution of the baryon density at the collision center and of the total 
strangeness per particle (i.e. the sum of s and s quarks per particle) for Au+Au reactions at 2, 10.6 
and 160 GeV/nucleon. The time-axis has been scaled in order to synchronize the time-evolution to the 
time of maximum baryon-density. At 2 GeV/nucleon only about 50% of the total strangeness has been 
produced. It saturates first in the late stages of the reaction at « 1.4p/p . At AGS energies already 60% 
of the total strangeness has been produced at the point of maximum density - saturation again takes 
places at « 1.5p/po- Between SIS- and AGS-energies no significant differences are visible for the time- 
evolution of strangeness saturation. At SIS- and low AGS-energies the dominant production mechanisms 
for strangeness are multi-step excitation of heavy resonances and their subsequent decay into a hyperon 
and a kaon. At AGS-energies string excitation may also contribute to strangeness production, but is not 
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Figure 4.39: Densities of (anti-)nucleons and (anti-)hyperons as functions of temperature for fixed 
/i q = 100 MeV and strangeness fraction / s = 0, calculated with a relativistic meson-baryon field theory 



(RMF) which implements hyperon-hyperon interactions (figure taken from [34]). 



yet a dominant factor. The situation changes for SPS-energies: Here already 75% percent of the total 
strangeness is produced at the time of maximum baryon density and saturation takes place at ~ 3p/po- 
The dominant production mechanism for strangeness at SPS energies is string-excitation. 



The time-evolution depicted in figure |4.40| characterizes the chemical freeze-out of strangeness in heavy 
collision systems. We shall now turn to the kinetic freeze-out, i.e. the time of the last interaction of the 
particle (i.e. scattering or its production via the decay of a resonance): 

Freeze-out radii are closely correlated to freeze-out times. Although the latter are definitely not mea- 
surable (for the radii, two-particle correlation techniques may yield similar information) they are the 
cleaner observable for a transport model analysis. Freeze-out radii may be contaminated by low energy 
particles produced close to the geometric center of the reaction during the late dilute expansion stage - 
such particles would have small freeze-out radii but large freeze-out times. 

The mass dependence of freeze-out times is studied in figure |4.41| . It shows the freeze-out time distri- 
butions for 7r, f] and u mesons in central Au+Au reactions at 2 GeV/nucleon. The heavier the meson, 
the earlier its freeze-out. The reason for this correlation is that heavy mesons are solely produced by the 
decay of massive baryon-resonances which can only be excited in the early hot and dense reaction zone. Of 
course freeze-out times do not only depend on the mass of the produced particle, but also on its interaction 



cross section: Figure ]4.42| shows freeze-out time distributions for K + , K~ and hyperons. Despite their 
large mass difference, the freeze-out times for hyperons and K~ are very similar, whereas the K + freeze 
out much earlier. This is due to the different quark content of K + and K~\ Both, K~ and hyperons 
contain a s-quark and have relatively large hadronic interaction cross sections. The K + , however, contains 
a s-quark and has a far lower cross section. 

So far, we have only studied the freeze-out behavior of mesons and hyperons at relatively low bom- 
barding energies around 2 GeV/nucleon. In that domain, baryon-baryon and meson-baryon interactions 
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Figure 4.40: Time evolution of the mean baryon density and total produced strangeness per participant 
baryon (UrQMD calculation). The mean baryon density is defined as the average over the baryon densities 
computed at the locations of all individual baryons. The total strangeness is the sum of all produced quarks 
and antiquarks. The time axis has been rescaled in order to overlay the respective times of highest baryon 
density (dashed grey line). 
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Figure 4.41: Normalized freeze-out time distri- 
butions for ir, t] and u mesons in UrQMD cen- 
tral Au+Au reactions at 2 GeV/nucleon. Heavy 
mesons freeze out earlier than light ones. 



Figure 4.42: Normalized freeze-out time distri- 
butions for K + , K~ and hyperons in UrQMD cen- 
tral Au+Au reactions at 2 GeV/nucleon. Particles 
with s-quarks freeze out earlier than those with s- 
quarks. 



are dominant and available phase-space may strongly influence the freeze-out behavior of the particles 
under investigation. Let us now turn to CERN/SPS energies - here meson production is so abundant that 
meson-meson and meson-baryon interactions dominate over baryon-baryon interaction and threshold or 
phase-space restrictions do not anymore apply. 

Unfortunately, neither the freeze-out density, nor the freeze-out time is directly observable. However, 
figure 4.43 shows that we can establish a correlation between high transverse momenta and early freeze-out 
times, at least in heavy colliding systems. In figure [4.43| the freeze-out time of pions is plotted versus their 
transverse momenta for p+p, S+S and Pb+Pb reactions at SPS energies. Naturally, the proton-proton 
system does not show any correlation, whereas in the heavy Pb+Pb system a strong p r dependence of 
the freeze-out time is visible. Selecting particles with high transverse momenta thus yields a sample of 
particles with predominantly early freeze-out times and high freeze-out densities. 

Do all hadron species exhibit a uniform freeze-out behavior - or does each species have its own compli- 
cated space-time dependent freeze-out profile, as observed at SIS-energies? Figures |L44| and |4.45| show the 
freeze-out radius and freeze-out time distributions for pions, kaons, antikaons and hyperons at mid-rapidity 
in central Pb+Pb reactions at 160 GeV/nucleon. The distributions have been normalized in order to com- 
pare the shapes and not the absolute values. In contrast to the situation at 2 GeV/nucleon, all meson 
species show surprisingly the same freeze-out behavior - the transverse freeze-out radius and freeze-out 
time distributions all closely resemble each other. Only the hyperons show an entirely different freeze-out 
behavior. Whereas the common freeze-out characteristics of the mesons seem to hint at a thermalization, 
the hyperons show that even at SPS energies there exists no common global freeze-out for all hadron 
species. 

Since the freeze-out distributions have a large width, the average freeze-out radius clearly does not 
define a freeze-out volume and therefore estimates of the reaction volume or energy density based on 
average freeze-out radii have to be regarded with great scepticism. 

In this section we have investigated freeze-out properties of hadrons produced in relativistic heavy-ion 
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Figure 4.43: Freeze-out time of pions as a func- 
tion of transverse momentum for p+p, S+S and 
Pb+Pb reactions at CERN/SPS energies. For 
heavy systems early freeze-out is correlated to high 
Pt- 



Figure 4.44: Normalized transverse freeze-out ra- 
dius distribution for pions, kaons, antikaons and 
hyperons. Whereas all mesons show similar distri- 
butions, the hyperons freeze out at larger radii. 



reactions. At SIS energies each particle species has its own complex freeze-out characteristics, governed 
by interaction cross sections, available phase-space and resonance lifetimes. At CERN/SPS energies, 
similar freeze-out distributions for various meson species hint at a system closer to thermal equilibrium. A 
comparison with the hyperon freeze-out distributions, however, indicates that no global thermalization and 
freeze-out has been achieved. Therefore, thermal model fits to hadron ratios and spectra - even though 
they work pretty well - may result in a misleading interpretation of the final state of a relativistic heavy 
ion reaction. 



4.2.8 Dilepton production 

One basic problem of heavy ion physics is actually how to investigate the hot and dense phase that is 
presumably formed in a collision of two relativistic nuclei. Since hadrons may interact several times, 
the intermediate stages are so to speak "shadowed" by the freeze-out distributions. Dileptons are of 
great interest, because they do not interact with the hadronic matter. Furthermore, they are emitted by 
various mechanisms at all stages of a heavy ion collision. Thus, the dilepton signal yields time integrated 
information on the reaction dynamics. 

Dileptons were initially proposed as highly penetrating probes of the QGP state | |286| , |287| , |288| | . They 
can be produced via a virtual photon in qq — > 7*g annihilations and in the QCD Compton process 
gq — > 7*g. If one assumes the QGP to be an equilibrated gas, then dileptons (as well as real photons) can 



probe its thermodynamic conditions [289]. However, due to large background contributions of hadronic 



radiation at low masses and the Drell-Yan process ||290|| at high masses, dileptons stemming from the QGP 
will possibly be overshined. 

Recent interest has focussed on low mass dileptons. They are emitted in hadronic decays and collisions. 
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Figure 4.45: Normalized freeze-out time distribu- 
tion for pions, kaons, antikaons and hyperons. As 
with the freeze-out radii, the times for the meson 
species are very similar. The hyperons again show 
a different behavior. 



Figure 4.46: quark condensate —(qq) func- 
tion of temperature T and baryon density p/po- 
The figure has been adapted from ||285| . 



Especially the Dalitz decays of light mesons and the direct decays of vector mesons are supposed to 
contribute significantly for invariant dilepton masses Mu < lGeV. In this mass region also various 
bremsstrahlung mechanisms might be important [51], [29 1| , |292|| as well as a number of direct production 
channels, e.g. irp — ► lift ||293| . Because the vector mesons have rather small life-times, they can supposedly 
resolve the rapid changes during the hot and dense phases of the collision. Their detection can proceed 
via low-mass dileptons. 



296|| , the QCD condensates (e.g. 



The dependence of (qq) on the 



In conjunction with the chiral symmetry restoration 
(qq)) should lower their values at high temperatures and/or densities, 
temperature T has been studied in the framework of lattice QCD [ |297| ] and chiral perturbation theory 
298|| . Up to 0.7 - 0.8Tc, (qq) remains nearly constant and then its absolute value decreases rapidly (see 
figure ^4.46|) . The behavior of the quark condensate at finite baryon densities is described in a model 
independent fashion by the Hellman-Feynman-theorem ||29S |. A model calculation of the dependence of 
(qq) on both, the baryon density p/p and temperature T, can be seen in figure |4.46| - the drop of (qq) 
with p and T is quite analogous to the temperature and density dependence of the nucleon effective mass 
in the a — to model as noted in ||275| . 

The connection between the QCD condensates and phenomenological hadronic quantities (masses, 
width etc.) can be drawn via the QCD sum rule technique |[300| , |301| , |302|| . Under simplifying assumptions 
on the spectral function, such calculations find lowering vector meson masses as indicators of chiral sym- 
metry restoration [|303| , |304|| . More refined calculations, which try to evaluate the medium contributions to 
the self energy of the vector mesons, find broadening spectral functions in matter |S5], |87|, ^ [Jt], [9T| . 

Various experiments have focussed on low mass lepton pairs: the DLS spectrometer at the BEVALAC 
301 HI, the CERES 



3071 , |308| and HELIOS [gOg |32j experiments at the SPS at CERN. The dilepton 



spectrometers HADES at SIS (GSI) |3ig and PHENIX at RHIC (BNL) |JTTJ are under construction. 
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Experiments at CERN have reported an enhanced production of dilepton pairs in nucleus-nucleus 
collisions over a broad invariant mass region around M ~ 0.5 GeV [312]. The enhancement is relative to 
the known sources as measured in pp or pA collisions after scaling to the nucleus-nucleus case (see l.h.s. 



of fig. f4.47|) . The data have been however reproduced by assuming density dependent masses of vector 
mesons as a consequence of partial chiral symmetry restoration | |313| , |314| , |315| , [5Tj] , |317| | (r.h.s. of fig. 
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Figure 4.47: Inclusive e + e~ mass spectra in 200 GeV/nucleon S+Au collisions as measured by the 
CERES collaboration [[307]]. The figures have been taken from ||312|| . The shaded area depicts hadronic 
contributions from resonance decays. The l.h.s. shows a comparison with calculations based on a purely 
hadronic scenario ||317| , |318| , |313j , |315| , |3 1 9|| whereas the r.h.s shows calculations including either a QGP 



phase transition or medium dependent vector meson masses [|313| , |315|| 



However, even bare hadronic transport model calculations, without any mass shift included, miss only 
the data in the 400 to 600 MeV bins (by 2 to 3 standard deviations) ||316| , |317| . Since hadronic transport 
models did - so far - neglect some contributions, e.g. from bremsstrahlung, it has yet to be determined 
whether partial restoration of chiral symmetry is the only possible explanation of these interesting new 
data. Calculations evaluating in-medium spectral functions, due to the coupling of the p with nucleon 
resonances and particle-hole excitations, also achieve a satisfactory reproduction of the CERES data 189] . 
without requiring a dropping p-mass. One has to conclude that, up to now, comparisons of dilepton spectra 
with hadronic models do not give evidence for (partial) restoration of chiral symmetry [ |320| | . Another, 
up to now neglected contribution to the low mass dilepton spectrum in nucleus-nucleus collisions, are 
secondary Drell-Yan processes, which are described in greater detail in the next sub-section. 

The production of low mass dileptons is included in the UrQMD model via Dalitz decays of 7r°, rj, u, rj' 
mesons and of the A (1232) resonance, direct decays of neutral vector mesons and incoherent pn 
bremsstrahlung ]|321|| . This approach has been used to analyze data of the DLS collaboration at the 
BEVALAC accelerator ][305| , [306|| . The recent AA data all show a large enhancement as compared to 
model calculations (see fig. |4.48| ). The pd spectra imply a high pn — > rjX cross section, but TAPS 
data delimit this as an explanation for the enhancement. So far, these data resist any explanation by 
medium-dependent spectral functions ||321|, |322|| . Note that the recent data strongly exceed earlier pub- 



lished measurements of the same collaboration [ 323 which have been reproduced by various transport 
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models |324| , |325| , |322|| and are also favored by the UrQMD model. A satisfying explanation of the 
recent DLS data is still missing. 
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Figure 4.48: UrQMD dilepton cocktail plot for Ca+Ca in comparison to recent DLS data [|305|1 . The upper 
solid curve is the sum of all contributions. The detector acceptance has been modeled with the DLS Filter 
4.1 and a mass resolution of 10% is adopted. 



Drell-Yan pairs from secondary collisions 



Reports on intermediate mass muon pairs ||328| , |309| , |326|| produced in heavy ion collisions have recently 
attracted much attention. The measured spectrum is enhanced relative to the expected 'cocktail' dilepton 
sources (Fig. 4.49| , left), one of them being the Drell-Yan pair production ||29( 



Rate estimates of this process (e.g. ||329|| ) are essentially based on extrapolations of p+A reactions 



where a linear scaling with A is observed for the Drell-Yan production cross section p3 



a 



pA 



(T A. 



This linear scaling can be understood in the Glauber picture of the hadron-nucleus cross section, con- 



structed at high energies using the AGK cutting rules [331 



But also the hot matter produced in ultrarelativistic heavy ion collisions contributes to the total 
dilepton radiation ||332| , |333| , |334| , |335| , |336| , |337|| . It is shown that emission from thermal sources such 
as a quark-gluon plasma or hadron gas, whether in or out of equilibrium, cannot be neglected relative 
to nucleon-nucleon contributions to the Drell-Yan process for masses below the J/ip peak — at least at 
higher energies than those presently available. 

Another source of dileptons which may already play a role at SPS energies and could account for 
the observed intermediate mass dilepton enhancement has been proposed in |[327|| : Drell-Yan production 
by interactions involving produced, or secondary, hadrons. In order to investigate secondary dilepton 
production at SPS energies UrQMD has been employed to obtain a realistic collision spectrum of secondary 
hadrons. The differential Drell-Yan cross section is computed at leading order (LO) using the standard 



equation ||338 



d 2 a 
dm 2 dy 



(AB -> UX) 



A-rra 2 
9m 2 s 



e l [? A (za, m 2 )q B (x B , m 2 ) + q A {x A , m 2 )q B (x B , m 2 ) 



(4.7) 
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Figure 4.49: Left: Comparison of measured dimuon yield with expected sources in Pb+Pb collisions. 
Figure taken from ||326| . Right: Calculation of secondary Drell-Yan processes based on a UrQMD simu- 



lation. Shown is the ratio between the total dilepton mass spectrum and the nucleon-nucleon scattering 
alone, at y cms = 0.5, for Pb+Pb at 160 A GeV assuming different hadronic formation times rp (Figure 
taken from 11321 



where q(x,m 2 ) and q(x,m 2 ) denote the quark and antiquark densities (according to [|339| , |340|1 ), y/s is the 



center-of-mass energy of the colliding hadrons, m is the invariant mass of the lepton pair, xa = y/r e y 
and xb = \fr e~ y with r = m 2 /s, and y is the dilepton rapidity in the cms frame. The lepton pair 
production cross section is calculated for each hadron-hadron collision and weighted by the inverse of the 
total hadron-hadron cross section. The distributions from these elementary hh-collisions are then summed. 
It is clear that in pion-nucleon collisions, valence quark-antiquark annihilation can play a significant role 
in the Drell-Yan process. Pion-nucleon dilepton production cross sections are consequently higher than 
nucleon-nucleon cross sections, especially when mj \fs ~ 0.1. It is shown that these secondary collisions 
can serve as an important source of m ~ 2 GeV dileptons due to the availability of valence antiquarks in 
mesons and antibaryons (Fig. |4.49| right). 



The standard Drell-Yan process corresponds to the interaction of fully formed hadrons. However, it 
was shown ||341| , |342| , |343|| that, during the early stages of the system evolution, partons can scatter and 



annihilate before they have come on mass-shell. To estimate the importance of these "primordial" or 
"pre-resonance" qq annihilations, the contribution of such processes have been calculated assuming that 
the asymptotic parton distribution functions are also valid for the primordial states. This is done very 
simply in the UrQMD calculation by decreasing the formation time of the produced hadrons within string 
excitations, t f , from the "default" value of around 1 fm/c. 

The importance of this primordial (pre-resonance) contribution to the dilepton mass spectra is shown 
in Fig. |1.49j (right). The secondary dilepton yield for rp = Oincreases by a factor of ~ 5 at all masses 
compared to the calculations with the default Tp ~ 1 fm/c. With a value of Tp = 0.5 fm/c the enhancement 
in the range 1.5 < m < 2.5 GeV shows quite good agreement with the mass data which are about a factor 



of two higher than naively expected [|326], |328|, |309 



4.2.9 Charmonium production and suppression 

J/^'s, x' s an d ^"s, consisting of a cc pair, are formed in pp-collisions by the fusion of two gluons, or to a 



larger extend by the annihilation of a quark-antiquark pair | 345 |. In a QGP of sufficiently high density and 
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Figure 4.50: J/tjj production cross sections versus Ax B. Figure taken from ||344 | 



temperature cc pairs cannot bind together - the presence of gluons and quarks in the plasma causes the 
color charge of quarks to be screened. Alternatively, hard gluons may be able to dissociate the cc pair (QCD 
photo-effect) in the QGP. Suppression of quarkonia-mesons has been suggested as a possible signature for 
deconfmement ||346|| . Lattice SU(3) gauge theory ||347| , p48|| indicates that the necessary conditions for the 
screening are already fulfilled shortly above Tc (the screening length Ad is inversely proportional to the 
temperature). Inclusion of dynamical fermions decreases the screening length ||349| , |350| . 

However, the suppression of quarkonia meson production is no unambiguous signal for deconfmement. 
They are no weakly interacting probes which can traverse hadronic matter unhindered. On the contrary, 
heavy quarkonia mesons suffer from low production cross sections and rather large hadronic dissociation 
cross sections. I.e. charmonium is already suppressed in pA-collisions, where the conditions for a QGP 
are not fulfilled. Figure |4.50| depicts the suppression of J/\I/'s in pA- and AB-collisions. The line shows 
the suppression due to the nuclear constituents of the projectile and target nuclei. One can see that this 
effect reproduces the data from pp- up to SU-collisions, but to explain the new lead on lead data one needs 
additional suppression. Purely hadronic scenarios have been suggested [ 254" , 351 , 3"5*2 , |353| , 354 1, which 
can explain J/\I/ and suppression up to PbPb-collisions without the need for a deconfined phase. The 
most successful hadronic scenarios for J and suppression explicitly take interactions with produced 
particles into account and are referred to as comover models - the comovers being the produced mesons 
which are energetic enough to destroy a J/\l/ in the case of a collision. For the all mesons can be 
comovers, but there is a threshold for the reaction J/\l/ + X — > DD, because 2m£> — nij/^ m 640 MeV. 
So, "slow" pions (in the J/^f rest frame) cannot destroy a J/\l/. 

Hadronic transport model calculations which incorporate the full collision dynamics and go far beyond 
the commonly used simplified version of the Glauber theory yield conflicting results ||355| , |263|| . These 
transport model calculations are very sensitive to certain input parameters such as the formation time 
of the J/\l/ and the comovers. The HSD transport model [[355|| can fully reproduce the NA50 lead data 
while assuming a fixed formation time of 0.7-0.8 fm/c for both, J/^> and comovers. A second calculation 



101 



with the HSD model using a formation time of fm/c for the J/\l/ reproduces also the data [ 356|| . The 
UrQMD model, however, uses for the comovers a variable formation time emerging from the Lund string 
fragmentation formalism (here the formation time depends on the hadron mass) and zero formation time 
for the J/ty. The assumption of zero formation time is valid if the J/\l/ is considered as a pre-resonance 
ccg — g state with a hadronic dissociation cross section of 7 mb. However, the UrQMD model does not 



reproduce the additional suppression of the Pb+Pb experiment |263|. The question of formation time 



might be a central issue since in the color octet model the dissociation cross section is actually higher 
during the lifetime of the pre-resonance cc$ — g state ||357| , |358|| than after hadronization. Furthermore the 
amount of comover-charmonium interaction will crucially depend on the formation time of the comovers. 

One has to bear in mind, however, that hadronic transport models do not contain partonic degrees of 
freedom explicitly and are therefore incapable of correctly describing hard processes like charm production 
or the Drell-Yan process. From the measurements of Drell-Yan lepton pairs we know that the cross section 
in pA-collisions is proportional to A, as predicted by the Glauber model. Hadronic models on the other 
hand take into account the energy loss of incident nucleons due to the production of secondaries (which 
then could play an important role as comovers). If, however, the elementary hard scattering cross section 
is strongly energy dependent - as the charm production at SPS energies and below - the effect of nuclear 
stopping leads to a considerable underestimation of the AB cross section. So, at the moment it is not 
clear, how it is possible to include hard processes consistently in a hadronic transport model. 

Another point under discussion is whether charmonium states might be created as a pre-resonance 
cc§ — g state [|343|| , which is not an eigenstate of the QCD Hamiltonian. Thus, one would expect a time 
evolution of this state and its cross section until it is a physical state, i.e a J . The pre-resonance should 
evolve during the formation time which would lead to different effective absorption cross sections in small 
and heavy nuclei. 

To summarize this section, we can say that the suppression of charmonium states in heavy ion collisions 
and even in proton-nucleus is not yet understood. Although, the assumption of a pre-resonance state with 
a constant cross section seems to reproduce the data, it contradicts principle ideas of quantum mechanics. 



4.3 Collective flow 



The excitation function of transverse collective flow is the earliest predicted signature for probing the 
formation of compressed nuclear matter |150 



It has been shown that the excitation function of flow is 



sensitive to the EoS and can be used to search for abnormal matter states and phase transitions 

In the fluid dynamical approach the transverse collective flow is directly linked to the pressure of the 
matter in the reaction zone. With P(p, S) being the pressure (depending on the density p and the entropy 
S) the generated collective transverse momentum can be written as an integral of the pressure over surface 
and time 



P± = JJ A P( P ,S)dAdt. 



where dA represents the surface element between the participant and spectator matters and the total 
pressure is the sum of the potential pressure and the kinetic pressure: The transverse collective flow 
depends directly on the equation of state, P(p, S). 

Collective flow has been predicted by nuclear fluid dynamics (NFD) [ffl, |150| , |360| , |361| , |362| , |TBJ. It is well 
established experimentally at the BEVALAC ||363 , 



365] for charged particles by the Plastic-Ball and 



streamer chamber collaborations and at GSI by the FOPI collaboration ||366|| . Microscopic models such as 
VUU (Vlasov-Uehling-Uhlenbeck) and QMD (Quantum Molecular Dynamics) have predicted smaller flow 
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than NFD. These models, however, show good agreement with the experimental findings [^, |152| , |154 



202 , 367 1 . One has to distinguish between different signatures of collective flow: The undirected radial 



flow which is best observed in ultracentral collisions and the bounce-off [|360|| of compressed matter in 



the reaction plane as well as the squeeze-out [361 of the participant matter out of the reaction plane in 
semi-central to semi-peripheral collisions. 

The most strongly stopped, compressed matter around mid-rapidity is seen directly in the squeeze-out 
7% . A strong dependence of these collective effects on the nuclear equation of state is predicted R202|| . For 



higher beam energies, however, projectile and target spectator decouple quickly from the reaction zone, 
giving way to a preferential emission of matter in the reaction plane, even at mid-rapidity ||368|| . 



Apart from the above discussed directed flow, the so-called "radial" , i.e. undirected, flow component 
can be used for simplicity (spherical symmetry) [ |197| , |198 |. It changes drastically the interpretation of 



particle spectra used for temperature extraction which may drop by as much as a factor of 2 (see the 
l.h.s. of figure |4.51| ). The mean (undirected) transverse velocity ((/?*)) or momentum ({pt)) can be used 



as a measure for this undirected flow in central collisions. The r.h.s. of figure |4.51| shows a UrQMD 
excitation function for (p t ). Up to an incident beam energy of 2 GeV/nucleon, (p t ) rises steadily and 
then almost saturates. The inclusion of potentials does not play a major role for this observable. The 
saturation at 2 GeV/nucleon (which is also observed experimentally, see the l.h.s. of the figure) can be 
explained in the following scenario: Up to 2 GeV/nucleon the longitudinal momentum of the incoming 
nucleons is predominantly transferred into transverse degrees of freedom; particle production only accounts 
for a minor fraction of the total energy of the system at the end of the reaction. Above 2 GeV/nucleon, 
however, particle production gets more important and finally most of the energy of the final state is 
stored in produced particles, leading to a saturation in the amount of energy carried by the "primordial" 
constituents of the collision system, the nucleons. The increase of energy stored in produced particles can 



be seen on the r.h.s. of figure |4.51| - here the excitation function for total energy of all produced particles 
(diamonds) is plotted into the same frame as the (pt) excitation function. 

The origin of nuclear flow in ultrarelativistic heavy ion collisions can be seen in Fig. |4.52| , where the 
time evolution of the longitudinal and the transverse pressure in the central region of Au(llGeV)Au (b=0), 
simulated with UrQMD, is depicted. A cylindrical volume of r = 6 fm and Az = 2 fm in the center of the 
reaction is considered. The pressure is calculated according to the virial theorem: 

1 Nevents ^particles 

vp l = = — E E v\v\lvl , 

iN events n= i k=1 

where i — x,y,z denotes the three spatial directions. The contribution of this kinetic part and a part from 
quasi-potentials have been studied separately in the framework of RQMD ||372 . It takes about 5 fm/c until 



the longitudinal pressure reaches its maximum value. The longitudinal pressure has decreased considerable 
at that time. Mainly not due to scattering and momentum transfer but due to the particles that simply 
have passed the central region and are therefore not considered any more. If all particles are taken into 
account, the ratio of transverse to longitudinal pressure does not exceed 0.5. This is in line with the 
observation of larger longitudinal than transverse flow after freeze-out. 

4.3.1 Bounce— off: collective flow in the reaction plane 

Due to its direct dependence on the EoS, P(p, T), flow excitation functions can provide unique information 
about phase transitions: The formation of abnormal nuclear matter, e.g., yields a reduction of the collective 
flow fL5| . A directed flow excitation function as signature of the phase transition into the QGP has been 
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proposed by several authors |7|, A microscopic analysis showed that the existence of a first order phase 
transition can show up as a reduction in the directed transverse flow [71|. 

For first order phase transitions, the pressure remains constant in the region of the phase coexistence. 
This results in a vanishing velocity of sound c s = \jdpjde. The expansion of the system is driven by the 
pressure gradients, therefore expansion depends crucially on c\. Matter in the mixed phase expands less 
rapidly than a hadron gas at the same energy density and entropy. In case of rapid changes in the EoS 
without phase transition, the pressure gradients are finite, but still smaller than for an ideal gas EoS, and 



therefore the system expands more slowly ||374| , |375| . 

This reduction of c\ in the transition region is commonly referred to as softening of the EoS. The 
respective region of energy densities has been called the soft region [|376| , |377| , |373| , |378|| . Here the flow will 
temporarily slow down (or possibly even stall). Consequently a time delay is expected in the expansion 
of the system. This prevents the deflection of spectator matter (the bounce-off) and, therefore, causes a 
reduction of the directed transverse flow [17, 113]. The softening of the EoS should be observable in the 



excitation function of the transverse directed flow of baryons (see figure |4.53| ). 

The overall decrease of p x seen in Fig. [4.53| for Eia, > 2 AGeV both for the hadronic and the QGP 
equation of state demonstrates that faster spectators are less easily deflected (because A and t in equa- 
tion £0| are decreasing with Ei ah ) by the hot, expanding participant matter. For the QGP equation of 
state, however, these one-fluid calculations show a local minimum in the excitation function, at about 
6 GeV/nucleon. This can be related to the QGP phase transition, i.e. to the existence of the soft region 
in the EoS. 

However, one-fluid hydrodynamic calculations assume instantaneous thermalization. This becomes 
unrealistic for increasing beam energies since due to the average rapidity loss of only one unit per proton- 
proton collision, nucleons require several collisions for thermalization. A more realistic three-fluid calcu- 
lation without a first order phase-transition, in which only local thermal equilibrium within each fluid is 
assumed, yields similar flow values as the one fluid model with a phase transition (solid squares in fig. 
[4.53|) |p79| , [38011 . The position of the minimum (the magnitude of the overall effect) therefore strongly 
depends on the degree of stopping (i.e. which type of fluid-dynamical model is employed) and on the 
details of the chosen EoS and phase transition parameters. 

Moreover, for finite volumes, V < 100 fm 3 , corresponding to expected plasma volumes in heavy ion 
collisions, there is a considerable rounding in the variables e/T 4 and s/T 3 around Tc, even if there is a 
first-order phase transition for the infinite volume limit |281| . This is inferred, under simple assumptions, 
from basic thermodynamic considerations: Fluctuations of the two phases in a finite system lead to a 
smooth transition between the low temperature regime — where the hadronic phase dominates the system 
- and the high temperature regime — where the pure quark phase is most probable. Such a behavior 
has severe implications on the proposed signal as was shown in [|378j| : a smooth crossover transition within 
an assumed interval of AT = O.IT^ results in drastically reduced time delays as compared to a sharp 
transition. 

Another effect is present if explicit finite size effects and the necessary requirement of exact color- 
singletness within the quark phase ||381] , |382|| is included ||281 |. The model exhibits a barrier in the free 
energy between the two phases near the phase transition. This leads to a shift of the critical temperature 
to higher temperatures for finite volumes (see Fig. |4.54j ). The speed of sound is considerably increased in 
the mixed phase. The significance of the time delay signal for the experimental detection of a QGP phase 
in heavy ion collisions, in turn, becomes questionable. 

A second order phase transition may not exhibit this minimum in the flow excitation function: The 
existence of a minimum in p x ,dir{Ei a b) is rather a qualitative signal for a strong first order transition. If 
such a drop of p x ,dir{Ei a b) is observed, it remains to be seen which phase transition caused this behavior: a 
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hadron-quark-gluon phase transition or, e.g., a nuclear matter to resonance matter transition in confined 
hadronic matter ||249 , [383|| . 

What does the purely hadronic UrQMD model predict with respect to directed transverse flow? Fig- 
ure ^J)5| shows the averaged in plane transverse momenta for Ni+Ni and Au+Au in the 0.1-4 GeV/nucleon 
region which is accessible through experiments at SIS and AGS. Calculations employing a hard equation 
of state (full symbols) are compared to cascade simulations (open symbols). In the latter case only a 
slight mass dependence is observed. For the calculation with potentials the integrated directed transverse 
momentum push per baryon is more than twice as high for the heavier system which corroborates the 
importance of a non-trivial equation of state of hadronic matter. 

For beam energies below 5 GeV/nucleon the amount of directed transverse momentum scales in the 
very same way as the total transverse momentum in the course of the reaction. Therefore the (p x ) versus 
rapidity divided by the (pt) is identical for all beam energies in the range of scaling as can be seen in 
figure [4.56 . 



4.3.2 Squeeze— out: flow perpendicular to the reaction plane 



The most strongly stopped, compressed matter at mid-rapidity is responsible for the squeeze-out []7T 
strong dependence on the nuclear equation of state for this collective effect is seen |202| , 384 1. 



A 



Let us now show the dependence of the observed squeeze-out on increasing beam energy, transverse 



momentum and impact parameter. We define a squeeze-out ratio [[fl|] 



out /in 
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For values R ou t/in > 1 neutrons are emitted preferentially perpendicular to the reaction plane. 

The top frame of figure [4.57| shows R ou t/in for Au+Au collisions versus scaled transverse momentum 
Pt/Pproj for beam energies of 400, 600, 800 and 1000 MeV/nucleon with cuts on rapidity ( —0.15 < 
Dcm < 0.15) and impact parameter (3 < b [fm] < 9). The ratio increases monotonously with the 
transverse momentum and is independent of the beam energy if the transverse momenta are scaled with 
the projectile momenta. Therefore, higher transverse momenta must be probed with increasing beam 
energy if a clean squeeze-out signal is to be isolated. Such a scaling behavior has already been extracted 



for the in-plane bounce-off 'in the hydrodynamic model [|385|1 . It has been experimentally observed in the 
case of the neutron squeeze-out by the LAND collaboration ||386|| . The respective comparison between the 
IQMD model and LAND data (mid-rapidity, ERAT2 centrality criterion) is shown for 400 MeV/nucleon 
in the lower frame of figure |4.57| . One should note, however, that in the IQMD calculation the reaction 
plane is always known, whereas the experimental determination of the reaction plane yields fluctuations 
which might underestimate the measured ratio by 15 to 30%. The p t dependence of the squeeze-out ratio 
is even enhanced for heavier clusters, as can be seen in figure 4.59| . 

The dependence of the neutron squeeze-out on the momentum dependent interaction and the equation 
of state is shown in figure |4.58j (for Au+Au collisions at 1 GeV/nucleon): The upper frame shows the hard 
equation of state with and without momentum dependent interaction (mdi). With mdi, R ou t/in increases 
by 50% for large transverse momenta! A difference of 50% can also be observed by comparing the hard and 



soft equations of state (both including mdi) in the lower frame of figure |4.58| : At high transverse momenta 
the hard equation of state with mdi exhibits a 50% higher squeeze-out ratio than the soft equation of 
state with mdi. The hard equation of state without mdi shows a R ou t/in vs. Pt dependence similar to the 
soft equation of state with mdi. It is important to note, however, that these differences are only seen for 



105 



high transverse momenta pt/p pro j > 1- Analyzing p t -integrated data will therefore severely diminish the 
sensitivity towards the nuclear equation of state and the momentum dependent interaction. 

In fig. |4.60] , the excitation function of the coefficient v% is shown representing the so-called elliptical 



flow. t> 2 is extracted from a fit to the azimuthal distribution of nucleons according to dN/dcf) = Vq[1 + 
v% cos(0) + v 2 cos (20)]. Positive values of t>2 correspond to a preferential in-plane enhancement of the 
emitted particles while negative values describe preferred emission perpendicular to the reaction plane. 

Clearly, large differences are seen when comparing calculations with and without potentials. On the 
one hand, cascade calculations only show in-plane enhancement of the emitted particles. On the other 
hand, calculations including nucleonic potentials show a transition from in-plane to out-of-plane emission 
with increasing bombarding energy which is seen by current experiments of the EOS collaboration 
This big sensitivity to the model ingredients gives a handle to pin them down further. 
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Figure 4.51: Left: Excitation function of temperature T and average transverse expansion velocity (3f The 
figure has been taken from ||366|| . The data are from ||369| , |370| , |371| , [220j , |221|| . Right: Excitation function 
of the mean transverse momentum (squares and circles, left ordinate) in comparison with the total energy 
of produced particles (diamonds, right ordinate). Shown are the UrQMD results for Au+Au collisions 
with impact parameter 6 = 4 fm. 
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Figure 4.52: Time evolution of the longitudinal and the transverse pressure in the central region of a 
Au(llGeV)Au collision (b=0) calculated with UrQMD. 
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Figure 4.53: Excitation function of directed transverse flow, calculated in the framework of nuclear 
hydrodynamics [18, 373 , with and without deconfinement phase transition. In the case of a phase transition 
a minimum in the excitation function is clearly visible. Solid squares correspond to a 3-fluid model 
calculation with ideal nucleon gas EoS. 
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Figure 4.54: Left: Shift of the critical temperature ATq vs. the systems size V. The bag constant is 
B 1 / 4 = 200 MeV. The color singlet constraint is taken into account for the QGP equation of state. Right: 
Speed of sound (squared) c 2 s as a function of energy density e for three different cases: 

1) infinite volume of the system (full line). 

2) V = 100 fm 3 using the infinite matter EoS (dotted). 

3) V — 100 fm 3 using the EoS with color singlet constraint (dashed). 
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clS du function of the scaled rapidity for Au+Au at 
b=4 fm (UrQMD calculation). If the transverse 
flow is scaled with the mean transverse momen- 
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barding energy. 



110 



Au(XAMeV)Au 



O 400 MeV/nucl. 

□ 600 MeV/nucl. 

A 800 MeV/nucl. 

1000 MeV/nucl. 



O 400 MeV/nucl. IQMD 
■ 400 MeV/nucl. LAND 



).() 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
Pt/Pproj. 



1.6 



3.5 



3.0 



2.5 



ca 2.0 



1.5 



3.5 



3.0 



2.5 



P5 2 



1.5 



1.0 



Au(lAGeV)Au 



O hard eos with mdi 
□ hard cos without mdi 



3> 



® 



O hard cos with mdi 
□ soft eos with mdi 



3> 



CD 



% m 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 
Pt/Pproj. 



Figure 4.57: Squeeze-out ratio R ou t/in versus 
scaled transverse momentum Pt/p pro j for neutrons 
in Au+Au collisions at 400, 600, 800 and 1000 
MeV/nucleon calculated with a hard equation of 
state without momentum dependent interaction 
(top) and a comparison between the IQMD cal- 
culation and data by the LAND collaboration for 
400 MeV/nucleon (bottom). By scaling p t with 
Pproj. the ratio becomes independent of the inci- 
dent beam energy (top). The binning is identical 
for all 4 energies, the symbols were shifted to in- 
crease the readability of the figure. The calculation 
shows good agreement with the data (bottom). 



Figure 4.58: Squeeze-out ratio R ou t/in versus 
scaled transverse momentum Pt/p pro j for neutrons 
in Au+Au collisions at 1 GeV/nucl.. The upper 
frame shows a comparison between IQMD calcula- 
tions with a hard equation of state with and with- 
out momentum dependent interaction (mdi). For 
large transverse momenta pt/p pr oj > 1 the calcula- 
tion with mdi exhibits a 50% larger squeeze-out ra- 
tio. The lower frame shows a comparison between 
calculations with hard equation of state with mdi 
and soft equation of state with mdi: Again, for 
large transverse momenta Pt/ Pproj > 1 the hard 
equation of state shows a 50% higher squeeze-out 
ratio than the soft equation of state. 
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Figure 4.59: Squeeze-out ratio R ou t/in versus fragment mass Af for Au+Au collisions at 400 MeV/nucleon 
(IQMD calculation). The system has been rotated by Qfi ow into the principal axis system, which enhances 
Rout/in- The rise of R ou t/in with increasing fragment mass is clearly visible, especially for high transverse 
momenta p t /Af. The figure has been taken from ||384 1 . 
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Chapter 5 

Summary and conclusions 



This review has by no means covered the entire field of microscopic transport model applications to (ultra- 
)relativistic heavy ion collisions. To cover the entire field in one review is almost impossible due to its fast 
development pace and the width of applications which have been added in the last couple of years. Instead, 
we have selected a sample of topics which allow to elucidate different important features of microscopic 
non-equilibrium transport theory. Of course a "good" transport model must be able to come close to 
known correct data. It should predict the outcome of future experiments. However, it's usefulness must 
stretch far beyond that of a mere event-generator: it allows for the analysis of the underlying physics (the 
"input" of the transport model) and shows how the various components contribute to the final result of 
the calculation. 

In chapter ^| the concepts of microscopic transport theory have been introduced and the features and 
shortcomings of the most commonly used ansatzes are discussed. For pedagogical reasons transport theory 
is first discussed extensively for the non-relativistic case and the relativistic extension is introduced in a 
subsequent section. An introduction to the Quantum Molecular Dynamics model concludes this chapter. 

The UrQMD transport model has been described in great detail in chapter |||. Based on the same prin- 
ciples as QMD and RQMD it incorporates a vastly extended collision term with full baryon-antibaryon 
symmetry, 55 baryon and 32 meson species. Isospin is explicitly treated for all hadrons. The range of ap- 
plicability stretches from Ei ab < 100 MeV/nucleon up to Ei ab > 200 GeV/nucleon, allowing for a consistent 
calculation of excitation functions from the intermediate energy domain up to ultrarelativistic energies. 
The UrQMD model as described in chapter |3] is no static structure. Its most important purpose is to serve 
as a framework into which new transport theoretical concepts and physics ideas can be incorporated, so 
that their effects on (ultra-)relativistic heavy ion reactions may be explored. 

Among the tasks to be tackled in the near future are the implementation of a relativistic generalization 
of n-body forces, an improved treatment of the imaginary part of the relativistic hadron self-energy (off- 
shell propagation), medium-dependent cross sections, an improved treatment of resonance lifetimes and 
time-delays as well as the inclusion of partonic degrees of freedom. The latter is important for the extension 
and continued application of transport models at collider energies (RHIC and LHC). The intriguing role 
of color coherence phenomena (transparency and opacity) for fluctuations, stopping and charmonium 
production can be studied best at these energies. However, even at these high energies, hadronic degrees 
of freedom must not be neglected since hadronic interactions in the late reaction phase may considerably 
change the hadrochemical cocktail and phase space distribution of the system (e.g. via feeding). 

Applications of transport models to (ultra-)relativistic heavy ion collisions have then been reviewed in 
chapter |j. The main topics under discussion have been stopping, particle production and collective flow. 
The main emphasis in the stopping section has been on the scaling behavior of the rapidity distribution 
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with respect to the incident beam energy and collision system mass as well as the dependence of the 
stopping behavior on the (di-)quark hadron cross section (i.e. the formation time dependence). 

The section on particle production has carried the largest weight in this review: subthreshold production 
of antiprotons has been used as an example of subthreshold particle production in heavy ion collisions - 
the ideas and concepts being very similar for other particle species such as etas and kaons. Multi-step 
processes and the excitation of heavy resonances dominate this domain of particle production. Therefore 
subthreshold particle production is a sensitive probe for the investigation of in-medium cross sections 
and the nuclear equation of state. Single particle spectra have been discussed for protons, deuterons, 
tritons and pions. The proton and cluster spectra have been used to point out problems with the source 
temperature extraction from single particle spectra. Pion energy spectra show how different regions of the 
spectra may be sensitive to different sources and reaction stages. 

Particle ratios have been used to discuss problems of temperature and chemical potential extraction 
via chemical equilibrium assumptions. The rapidity dependence of certain particle ratios may allow for 
the in-depth investigation of the baryon stopping mechanism in heavy ion collisions. The UrQMD exhibits 
the properties of a free hadron gas, however with a limiting temperature of approximately 140 MeV. The 
central cell in Au+Au reactions at AGS seems to be close to local thermal and chemical equilibrium. 

The formation of resonance matter has been discussed in great detail. The fraction of excited matter in 
the collision system rises continuously from 20% at SIS up to 70% at CERN SPS energies. The properties 
and interactions of highly excited resonance matter is an important topic for further research. Creation of 
strange matter is discussed, although explicit Strang elet-formaiion has yet to be included into microscopic 
transport model calculations. Hyperclusters and ( ant i-) deuterons have served as examples for the formation 
of composite particles. The main emphasis of the discussion has been on the phase space distributions of 
anti-deuterons and the sensitivity of the respective (centrality dependent) yields on parameters such as 
the formation time and the treatment of the baryon-antibaryon annihilation cross section. 

Freeze-out distributions of mesons and hyperons at SIS and SPS energies show that global equilibration 
is not achieved in heavy ion reactions. The distributions indicate a complex non-equilibrium time evolution 
of the hadronic system. A thermal model although quite successful in fitting the data, may not be yielding 
the correct physics interpretation. 

Electromagnetic probes allow for an unhindered view into the hot and dense reaction zone. This review 
focused on dilepton production (both at BEVALAC and SPS energies) and the creation of secondary Drell- 
Yan pairs at CERN/SPS energies, two of the most current topics. The latter offer a novel explanation for 
the observed dimuon excess in the intermediate mass range of 1.5 < m < 2.5 GeV of the dimuon-spectrum. 
A brief discussion of charmonium production and suppression concludes the section on particle production. 

The final section of chapter |4] deals with collective flow. The excitation function of the directed flow 
in the reaction plane can provide a signature for the hadron gas to quark gluon plasma phase transition. 
However, finite size effects may severely diminish the signal. Transport model calculations (which in 
most cases do not contain a phase transition) serve here again as an important smooth baseline to check 
the occurrence of possible QGP signatures, namely as irregularities in the excitation function, taken in 
small steps, AE m 5 GeV or so. The sensitivity of flow perpendicular to the reaction plane to the nuclear 
equation of state is analyzed. Especially the squeeze-out excitation function seems to offer new possibilities 
for the investigation of the density dependence of the nuclear equation of state. 

Microscopic transport models like UrQMD are a unique tool to further develop the general understand- 
ing of the dynamics of heavy ion collisions and interactions over a vast energy range from the Coulomb 
barrier (several MeV per nucleon) to the highest energies currently available or planned for the future. 
Their development is by far not complete. Nevertheless they allow for important insights into the field of 
hot, dense hadronic matter created in (ultra-)relativistic heavy ion reactions. 
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